AD- A 197  152 


"V W  1 

iJHLLAi  j  i  I  ILL) 


SECURI T  V  CLASSIFICATION  OF  THIS  PAGE  (Wffen  Oafa^  rntered) 

REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

*.  REPORT  NUMBER  2.  GOVT  ACCESSION  NO. 

AF  IT/C  I/NR  88-  l  74 

3.  RE9yy|Lys  catalog  number 

sM  IL  f  ilC  LW.i 

TITLE  (and  5u6fff(e) 

KJOWL/Mi^iv  HfF£cT5  no 

PrtoT6  (U  ndfccrvvt  Ctf-'y  STA15 

5.  TYPE  OF  REPORT  ft  PERIOD  COVERED 

MS  THESIS 

6.  PERFORMING  ORG.  REPORT  NUMBER 

AUTHORf.^ 

OdmALQ  (Ld  S3  £nSSCHtO£ 

8.  CONTRACT  OR  GRANT  NUMBERfs) 

PERFORMING  ORGANI  Z  ATION  NAME  AND  ADDRESS 

AFIT  STUDENT  AT:  Urti\»fcA.&  iTy  0  f 

OV  Fort-K  -  U.K, 

10.  PROGRAM  ELEMEN  T.  PROJECT.  TASK 

AREA  ft  WORK  UNIT  NUMBERS 

.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

_  .  _ _ 

12.  REPORT  DATE 

1988 

13.  NUMBER  OF  PAGES 

u  MONITORING  AGENCY  NAME  6  AOORESSf/f  dlllrrrnt  Iron:  Controlling  Otlice) 

AFIT/NR 

Wright-Patterson  AFB  OH  45433-6583 

15.  SECURITY  CLASS,  (ol  thla  rrport) 

UNCLASSIFIED 

15a.  DECL  ASSI  FI  C  ATI  ON/ DOWN  GR  AOIN  G 

schedule 

...  1 

16.  DISTRIBUTION  STATEMENT  (ol  Ihl.i  Rrport) 

DISTRIBUTED  UNLIMITED:  APPROVED  FOR  PUBLIC  RELEASE 

17.  DISTRIBUTION  STATEMENT  (of  the  abstract  entered  In  Block  20,  If  different  from  Report) 

SAME  AS  REPORT 

.8,  supplementary  notes  Approved  for  Publ  Ic^fel ea s JAW  AFR  190-1 

LYNN  E.  WOLAVER  ^ 

Dean  for  Research  /^ndProfessional  Development  ' 

Air  Force  Institute  of  Technology 

_ Wrjqht-Patterson  AFB  OH  45433-6583 _ 

19.  KEY  WORDS  (Continue  on  reverse  side  if  necessity  and  identify  by  block  number) 

20.  ABSTRACT  (Continue  on  reverse  side  ft  necessary  end  Identify  by  block  number)  **  ___  •  • 

ATTACHED  DT!C> 

^flecte|^ 

\  AUG  1  7 1988  a  I 

v,  *  Hi 

H 

DD  1  JAN  73  1473  EDITION  OF  1  NOV  65  IS  OBSOLETE  L7 ’CLASSIFIED 

0 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  f»J.en  Data  Enlerad) 


NONLINEAR  EFFECTS  IN 
PHOTOREFRACTIVE  CRYSTALS 

By 

Donald  Ross  Erbschloe 
Merton  College 
Trinity  Term,  1988 


A  THESIS 

SUBMITTED  FOR  THE  D.  PHIL.  DEGREE 
at  the  University  of  Oxford 


The  Department  of  Engineering  Science 
Parks  Road 
Oxford 


*  88  8  16  014 


i 


ABSTRACT 


Donald  Ross  Erbschloe  Thesis  submitted  for  D.Phil.  degree 

Merton  College  Trinity  Term  1988 


NONLINEAR  EFFECTS  IN  PHOTOREFRACTIVE  CRYSTALS 

^Photorefractive  crystals  are  materials  whose  index  of  refraction  is  altered  under 
illumination  by  light.  These  crystals  are  both  photoconductive  and  electrooptic. 
When  a  nonuniform  light  intensity  pattern  is  present  in  the  material,  photocarriers 
are  generated  and  redistributed,  creating  space  charge  electric  fields  which  change 
the  refractive  index  locally.  These  crystals  are  ideal  media  for  real-time  holography, 
and  applications  include  wave  amplification,  image  processing,  phase  conjugation, 
and  laser  beam  steering  for  optical  interconnects. 

This  thesis  invest  igates  many  novel  aspects  of  the  photorefractive  effect.  A  study 
of  nonreciprocal  behaviour  identifies  a  new  important  consideration  in  the  theory 
of  two-wave  mixing  between  counterpropagating  beams — namely  the  presence  of  a 
photocurrent,  or  frequency  detuning  between  the  beams  results  in  a  spatially  vary¬ 
ing  beam  coupling.  A  numerical  treatment  of  these  important  cases  provides  the 
first  systematic  theoretical  assessment  of  the  control  of  nonreciprocal  transmission 
and  phase  shift  in  iithium  niobate,  a  representative  photorefractive  crystal.  A  com¬ 
parison  between  crystal  types  suggests  candidates  for  nonreciprocal  applications 
such  as  an  optical  diode,  l/  o  ,  ;  ’-)1  ’  <--~ 

A  study  of  bismuth  silicon  oxide,  Bi^SiC^O)  as  the  active  gain  medium  in  an 
oscillator  reveals  a  novel  feature,  the  presence  of  a  light  intensity  threshold.  For  one 
crystal  sample  no  oscillation  occurred  for  incident  intensities  less  than  0.8  mW /cm2. 
A  surprising  new  result  is  the  appearance  of  higher  diffracted  orders  in  a  crystal 
sample  with  a  small  wedge  angle  (0.036°)  due  to  wave  mixing  between  an  incident 
beam  and  its  first  codirectional  multiple  reflection.  Several  applications  for  this 
new  means  of  obtaining  beam  interaction  are  discussed — including  the  study  of 
the  photorefractive  coupling  for  very  large  grating  spacings,  the  investigation  of 
transient  phenomena,  and  laser  beam  steering. 

A  new  gain  mechanism  for  wave  mixing  in  photorefractive  crystals  is  verified 
experimentally  using  the  unidirectional  ring  resonator  and  the  linear  resonator. 
The  new  mechanism  depends  on  two,  frequency  detuned  pump  beams  and  off-Bragg 
coupling.  A  comparison  with  the  traditional  gain  mechanism  for  wave  amplification, 
using  a  single  pump  beam,  reveals  that  the  intensity  of  the  oscillating  beam  with 
the  new  mechanism  can  be  four  times  larger  than  with  a  single  pump  beam.  We 
provide  the  first  quantitative  measurement  of  oscillation  in  the  linear  resonator 
using  BinSiCho- 
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Chapter  1 
Introduction 

A  page  of  history  is  worth  a  volume  of  logic. 

Olivet  Wendell  Holmes,  Jr. 

The  driving  mechanism  for  the  optical  properties  of  a  crystal  lies  in  the  regular 
arrangement  of  its  atoms  or  molecules.  It  seems  natural  that  crystals  should  be 
used  to  investigate  the  character  of  light  and  vice  versa.  Perhaps  this  is  why  the 
evolution  of  optics  is  coincident  with  (indeed,  often  inextricably  linked  to)  the  study 
of  crystals.  The  twentieth  century  has  witnessed  tremendous  parallel  advances  in 
technology  in  the  fields  of  optics  and  solid-state  physics.  Developments  such  as  the 
laser  and  semiconductor  devices  have  infiltrated  and  revolutionized  everyday  life. 
More  than  ever  we  exploit  that  which  we  purport  to  know. 

In  1966  [5]  an  investigation  of  nonlinear  optical  properties  of  a  crystal  of  lithium 
niobate  (LiNbOa)  opened  a  new  front  of  study — the  photorefractive  effect.  This 
effect  truly  embodies  the  union  of  crystals  and  optics.  In  this  process,  as  the  name 
suggests,  the  refractive  index  is  influenced  by  light.  A  photorefractive  material  ex- 
hibitr  both  photoconductivity  and  the  linear  electrooptic  effect.  In  tact  the  best 
photorefractive  materials  such  as  bismuth  silicon  oxide  (BiuSiOao)  have  moder¬ 
ate  to  large  electrooptic  coefficients,  large  photoconductivity  parameters  and  low 
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Figure  1:  Intensity  pattern  created  by  two  interfering  sinusoidal  waves. 

dark  conductivities.  How  a  (semi-)  permanent  grating  can  be  established  in  such 
a  material  is  summarized  below.  If  two  monochromatic  plane  waves  intersect  as  in 
Figure  1,  an  interference  pattern  is  established  in  which  the  intensity,  /,  is  modu¬ 
lated  according  to: 

I(z)  -  I0  (1  +  m  cos(Kz))  (1) 

where  T0  is  the  total  intensity,  K  is  the  length  of  the  grating  vector  formed  hv  the 
two  mixing  waves  and  m  is  the  intensity  modulation  ratio.  If  this  space- varying 
sinusoidal  intensity  pattern  falls  upon  a  photorefractive  crystal  the  following  pro¬ 
cesses  occur.  These  materials  are  semiconductors  so  where  the  irradiance  is  high, 
electrons  from  donor  atoms  are  excited  into  the  conduction  band.  By  either  drift 
or  diffusion  the  free  electrons  move  to  darker  regions  where  recombination  occurs. 
The  electrons  may  recombine  and  be  re-excited  many  times  before  final  trapping 
occurs.  Motion  of  charge  carriers  will  eventually  reach  an  equilibrium  when  drift 
and  diffusion  effects  are  balanced  by  internal  fields  caused  by  space  charge  effects. 
A  modulated  charge  density  builds  up  inside  the  crystal,  mimicking  the  modulated 
intensity  pattern.  The  success  of  the  imitation  will  depend  on  a  number  of  features, 
such  as  drift  and  diffusion  lengths,  strength  of  applied  field,  etc.,  and  in  the  steady 
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state  manifests  itself  in  a  phase  mismatch  angle  between  the  intensity  fringes  and 
the  charge  concentration.  The  gradient  of  the  charge  density  expresses  the  space 
charge  electric  field  associated  with  it.  This  field  will  vary  sinusoidally  in  space  and 
through  the  linear  electrooptic  effect,  the  index  of  refraction  is  modulated  and  a 
grating  is  formed. 

This  grating  is  a  phase  hologram  [115].  A  hologram  uses  a  reference  and  an 
“object”  beam  to  create  an  interference  pattern.  This  interference  pattern,  which  is 
stored  in  the  recording  medium,  contains  phase  and  amplitude  information  about 
the  two  “writing”  beams.  If  the  hologram  is  replayed  by  the  reference  beam  alone, 
part  of  the  reference  is  diffracted  and  reconstructs  the  object  beam.  We  may  con¬ 
sider  *his  static  or  stationary  holography,  and  this  describes  the  situation  for  steady- 
state  wave-mixing  in  photorefractive  materials.  However,  to  reach  steady-state  we 
must  first  go  through  a  complicated  dynamic  process.  The  writing  beams  undergo 
diffraction  as  the  grating  is  being  formed  thus  the  hologram  is  being  written  and 
read  at  the  same  time.  These  diffracted  beams  also  contribute  to  the  interference 
field  and  so  hologram  formation  exhibits  a  wealth  of  transient  effects  before  settling 
down  to  equilibrium. 

The  interest  and  effort  into  the  study  of  photorefractivity  has  grown  steadily  in 
its  22  year  history.  Only  a  handful  of  articles  and  papers  appeared  on  the  subject  in 
the  late  1960’s.  Now  such  publications  as  Physics  Reports  (1982),  The  Journal  of  the 
Optical  Society  of  America  (February  1986),  and  The  IEEE  Journal  of  Quantum 
Electronics  (August  1986,  December  1987)  devote  issues  to  it.  Photorefractivity 
even  merits  special  dedicated  conferences  such  as  a  1986  Meeting  on  Progress  in 
Scientific  Culture  at  Erice,  Sicily  and  the  July  1987  Los  Angeles  meeting  of  the 
Optical  Society  of  America. 

There  is  more  to  the  utility  of  photorefractive  materials  than  just  as  erasable, 
real-time  media  in  which  to  store  display  holograms.  These  materials  are  optically 
nonlinear,  which  means  the  polarisation  of  the  medium  is  not  strictly  proportional 
to  an  applied  electric  field.  (It  is  an  unfortunate  consequence  of  an  overly  useful 
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Figure  2:  Photorefractive  grating  formation:  (a)  Carrier  photogeneration,  (b)  Drift, 
diffusion  and  trapping  of  carriers,  (c)  Space  charge  distribution,  (d)  Space  charge 
electric  field  (gradient  of  space  charge  distribution)  may  be  out  of  phase,  with 
intensity. (e)  Modulated  refractive  index  (via  eiectrooptic  effect.) 
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term  like  “linear”  that  the  coefficients  of  the  linear  electrooptic  effect  can  be  directly 
related  to  the  second-order  optical  nonlinear  susceptibility  coefficients.  [135])  Two- 
and  four-wave  mixing  in  photorefractive  materials  yield  such  nonlinear  phenomena 
as  energy  exchange  between  the  interacting  waves  [76],  phase  conjugation  [24],  and 
nonreciprocal  effects  [140]. 

As  we  approach  the  twenty-first  century  our  information-hungry  society  de¬ 
mands  more  and  more  from  its  technologies.  There  is  a  never  ending  quest  for 
faster  computers  and  improved  and  increased  communications  capability.  Optical 
engineering  has  become  an  important  arena  for  priority  research.  Much  of  the  early 
thrust  of  this  research  was  to  develop  optical  analogues  of  electronic  components 
such  as  switches  (bistable  optical  devices  [111]),  optical  diodes  [86],  and  logical 
gates  [80],  but  increasingly,  applications  aim  at  areas  where  optics  has  inherent 
strengths,  for  example  parallel  processing  in  associative  memories  [26]  and  neural 
networks  [81].  Photorefractive  materials  are  versatile  and  figure  prominently  in 
many  proposed  schemes  for  the  devices  mentioned. 

Before  we  discuss  studies  of  photorefractive  crystals,  it  is  worthwhile  to  review 
their  place  in  the  history  of  crystal  optics. 


1.1  History  of  Crystal  Optics 

Crystals  were  certainly  known  and  valued  to  the  ancients  in  the  forms  of  gems 
and  precious  stones.  The  first  scientific  mention  of  the  regular  shape  of  a  crystal 
(garnet)  appears  in  Theophrastus’  (372-287BC)  Concerning  Stones  [32].  To  the 
geographer  Strabo  (63BC-  24AD),  pieces  of  quartz  were  so  reminiscent  of  ice  that 
he  proposed  a  common  origin.  Thus  our  word  u  crystal”  comes  from  his  description 
“npvaraWov” .  It  means  Hto  freeze,  or  to  congeal  with  frost”  [1].  Pliny  the  Elder 
propagated  this  myth  in  his  marvelous  compendium  of  misinformation,  Natural 
History.  However,  he  did  note  the  characteristic  six  angles  of  quartz  crystals  [10]. 

The  ice-crystal  explanation  dominated  until  the  18th  century.  In  this  regard 
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Strabo’s  view  was  to  crystals  what  the  Aristotalean  view  was  to  mechanics  and 
nearly  as  long-lived.  One  of  the  more  eloquent  expressions  of  the  theory  comes  from 
John  Wycliffe  (1320-1384)  [1] 

When  the  colde  northwynde  bloweth  harde  Christall  cometh  of  the  water 

As  with  Aristotalean  dynamics,  though,  there  were  a  few  dissenting  voices  through 
the  centuries.  Agricola  (Georg  Bauer,  1490-1555)  dismissed  the  idea  of  frozen  crys¬ 
tals  as  he  noted  the  presence  of  crystalline  materials  near  Hecla,  an  active  volcano 
[101].  The  fact  that  if  crystals  were  a  form  of  ice  it  was  different  from  typical 
water-ice  was  mentioned  by  Robert  Boyle  (1627-1691).  Essentially  his  argument 
was  that  if  one  dropped  a  crystal  and  a  piece  of  ice  into  some  water  there  will  be  a 
difference.  The  crystal  will  not  float.  A  charming  refutation  of  the  then-current  hy¬ 
pothesis  of  crystals  from  frozen  dew  appears  in  Philosophical  Transactions  in  1672. 
P.  Francisco  Lana,  while  travelling  in  Italy,  was  given  some  six-sided  crystals  found 
exclusively  in  some  narrow  crevices  in  a  large  meadow.  “I  was  told,  that  they  were 
produced  from  the  dews,  because  (forsooth!)  being  gather’d  over  night,  the  next 
morning  there  would  be  found  others  at  such  a  time  only,  when  it  was  a  serene  and 
dewy  sky  .  .  .  To  Lana’s  credit,  he  argued  [82] 

It  therefore  ought  not  to  be  affirm’d,  that  a  dewy  vapour  is  of  it  self  able 
to  be  form’d  into  a  solid  gem;  because,  if  that  were  so,  such  vapours 
being  easily  carried  by  any  motion  of  the  Air  from  those  narrow  places, 
and  falling  down  in  dew  far  from  the  same,  Crystals  would  be  formed 
in  those  other  places. 

In  1669  Steno  (Niels  Stenson,  1638-1686)  set  down  the  law  of  constancy  of 
interfacial  angles  for  quartz.  This  law  was  extended  in  the  classification  hundreds 
of  types  and  shapes  of  crystals  by  such  researchers  as  Dominic  Guglielmi  and  Jean 
Baptiste  Rome  Deslisle  [32].  Strabo’s  simple  concept  was  no  longer  adequate  to 
explain  this  complicated  situation.  The  ice  theory  melted  into  oblivion. 
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If  crystals  weren’t  ice,  what  were  they?  Robert  Hooke  (1635-1703)  spectulated 
in  his  Mtcrographia  that  crystals  are  ultimately  composed  of  spherical  particles  (51). 

...  all  these  regular  features  arise  onely  from  three  or  four  several 
positions  or  postures  of  globular  particles  .  .  .  and  this  I  have  ad 

oculum  demonstrated  with  a  company  of  bullets  and  some  few  other 
simple  bodies  .  .  .  even  almost  by  shaking  them. 

It  was  a  fortuitous  accident  which  established  the  science  of  crystallography.  Rene 
Just  Haiiy  (1743-1822)  dropped  and  broke  one  of  a  group  of  calcite  crystals  belong¬ 
ing  to  a  friend.  He  noticed  that  the  fragments  had  the  same  characteristic  rhomb 
shape.  According  to  Bishop  [10], 

After  breaking  up  many  more  calcite  crystals  of  differing  external  shapes 
he  concluded  that  all  were  composed  of  a  repetition  of  a  minute  struc¬ 
tural  unit  .  .  .  called  integrant  molecules. 

(Presumably  he  first  asked  his  friend’s  permission,  else  this  is  an  example  of  exper¬ 
imental  science  at  its  most  callous.) 

In  the  19th  century  the  study  of  crystal  structure  passed  on  to  the  mathemati¬ 
cians  with  research  into  crystal  classes  and  group  theory.  Experimental  work  was 
left  to  the  study  of  interactions  of  light  with  crystals.  As  early  as  1669,  Erasmus 
Bartholinus  (1625-1698)  observed  double  refraction  in  Icelandic  feldspar  [9].  In  try¬ 
ing  to  explain  this  phenomenon,  Christian  Huyghens  (1629-1695)  discovered  the 
polarisation  of  light.  Further  advances  in  study  of  polarised  light  in  crystals  were 
made  by  David  Brewster  [13]  (1781-1868)  and  Etienne  Louis  Malus  (1776-1812), 
culminating  in  the  discovery  in  1816  that  two  rays  polarised  at  right  angles  will 
not  interfere.  This  definitive  demonstration  that  light  waves  are  transverse  was 
made  by  Augustin  Jean  Fresnel  (1788-1827)  and  Dominique  Jean  Francois  Arago 
(1786-1853)  with  the  explanation  coming  from  Thomas  Young  (1773-1829).  Fresnel 
worked  out  the  details  of  wave  propagation  in  crystals,  a  feat  which  Born  and  Wolf 
call  “one  of  the  greatest  achievements  of  natural  science* [11]. 
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Polarisation  of  light  is  easily  demonstrable  and  it  appealed  greatly  to  the  most 
eminent  popularizer  of  science  in  the  19th  century,  Micheal  Faraday  (1791-1867). 
In  the  1840’s  he  conducted  a  series  of  experiments  investigating  the  influence  of 
electric  and  magnetic  fields  on  light.  He  discovered  the  Faraday  effect  (a  non¬ 
reciprocal  rotation  of  the  polarisation  axis  in  a  crystal  subjected  to  a  magnetic 
field)  in  1845  [33]  but  failed  to  find  a  comparable  effect  using  an  electric  field.  It 
took  until  1875  and  better  laboratory  equipment  for  this  latter  observation.  In 
that  year  John  Kerr  (1824-1907)  demonstrated  the  quadratic  electrooptic  effect[68], 
a  change  in  the  refractive  index  proportional  to  the  square  of  the  applied  electric 
field.  The  linear  electrooptic  effect  followed  in  1893  from  the  studies  of  Friedrich 
Carl  Alwin  Pockels  (1865-  1913)  [50]. 

For  centuries  crystals  had  been  used  to  divulge  the  properties  of  light.  The 
history  of  light  and  crystals  comes  full  circle  in  1912  when  light  in  the  form  of  X-rays 
was  used  to  investigate  the  fundamental  nature  of  crystals.  The  classic  experiment 
exhibiting  X-ray  diffraction  was  suggested  by  Max  von  Laue  and  performed  by  W. 
Friedrich  and  P.  Knipping.  Using  Laue  photographs  produced  by  a  sodium  chloride 
crystal,  the  father-son  team  of  W.H.  and  W.L.  Bragg  ascertained  its  structure.  The 
thorough  (and  personal)  history  of  the  early  years  in  X-ray  diffraction  is  told  in 
chapter  12  of  the  Braggs’  classic  text  The  Crystalline  State  [12]. 

A  resurgence  in  interest  in  crystals  came  with  this  new  technology.  The  timing 
was  right  for  the  1920’s  saw  the  first  stirrings  of  quantum  mechanics.  An  early 
success  of  the  theory  was  the  discovery  that  a  periodic  potential  function  leads 
to  allowed  energy  states  in  the  form  of  bands  [124].  This  finding  helped  explain 
electrical  conduction,  and  led  to  the  ultimate  exploitation  of  semiconductors. 


1.2  History  of  the  Photorefractive  Effect 


In  1966,  Ashkin,  Boyd  and  their  group  of  co-workers  at  Bell  Telephone  discovered 
“optically-induced  refractive  index  inhomogeneities”  in  crystals  of  lithium  niobate 
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(LiNbOa)  and  lithium  tantalate  (LiTaOa).  Because  of  their  known  nonlinear  optical 
properties,  the  crystals  were  being  tested  in  applications  such  as  second  harmonic 
generation  and  as  electrooptic  modulators.  Exposure  of  the  crystals  to  focused  or 
unfocused  visible  laser  light  created  nonuniform  changes  primarily  in  the  extraor¬ 
dinary  refractive  index.  If  left  alone  there  was  some  relaxation  of  the  change,  but 
effects  would  still  persist  for  days  or  months.  Heating  the  crystal  (ss  170°  U)  accel¬ 
erated  the  restoration  of  the  crystal  to  the  neutral  state.  They  surmised  that  the 
electrooptic  effect  could  explain  the  index  change  but  a  more  fundamental  mech¬ 
anism  was  needed  to  explain  the  fields  required  to  induce  the  electrooptic  effect. 
Among  the  mechanisms  suggested  were  pyroelectricity  (a  spontaneous  polarization 
due  to  a  temperature  change),  microdomain  reversals,  or  a  change  in  the  structure 
of  the  crystal  cell  [5]. 

Chen  found  in  1967  that  potassium  tantalate  niobate  (KTN)  suffered  laser  in¬ 
duced  “damage”  as  well.  In  his  experiment  a  crystal  was  illuminated  with  633  nm 
He-Ne  laser  light  while  a  transverse  external  electric  field  was  being  applied.  Chen 
felt  a  long-lived  space  charge  field  could  be  established  by  the  drift  of  photoex- 
cited  electrons  from  traps  in  the  illuminated  area  and  their  subsequent  recapture  in 
traps  in  darker  regions  [18j.  In  a  later  paper,  Chen,  LaMacchia  and  Fraser  demon¬ 
strated  the  ability  of  this  “optical  damage”  in  recording  a  volume  hologram  inside 
lithium  niobate.  The  crystal  had  excellent  resolution  properties,  greater  than  1600 
lines/mm  [20).  No  applied  field  was  necessary  and  the  question  as  to  the  driving 
mechanism  of  photoexcited  carriers  in  lithium  niobate  became  important.  Johnston 
in  1970  claimed  that  a  type  of  spontaneous  polarisation  similar  to  the  pyroelectric 
effect  was  responsible  [63].  This  idea  appealed  to  Townsend  and  LaMacchia  who 
created  volume  holograms  using  barium  titanate  (BaTiOs)  [125].  Chen  likened  the 
effect  to  persistent  internal  polarisation  [19],  an  effect  well  known  in  ferroelectric 
materials  in  which  strong  internal  fields  could  last  for  days  after  space  charge  fields 
were  established  by  using  radiation  to  excite  electrons  which  can  then  drift  out  of 
the  irradiated  area  and  be  trapped  [66]. 
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In  1971,  Amodei  repoited  that  thermal  diffusion  could  account  for  the  motion 
of  excited  free  carriers  [3].  Later  that  year  he  and  co-author,  Staebler,  published 
a  classic  paper  in  which  Kogelnik’s  coupled-wave  theory  was  applied  to  hologram 
writing  and  reading.  Both  diffusion  and  drift  of  carriers  were  considered.  They  also 
predicted  a  one-way  energy  transfer  between  recording  beams  [119]. 

With  this  1971  paper,  a  coherent  picture  of  what  was  later  to  be  called  the 
photorefractive  effect  began  to  emerge.  (The  earliest  reference  I  have  found  in 
which  the  term  “photorefractive”  is  used  is  in  a  1976  article  by  Magnusson  and 
Gaylord  [87].) 

The  1970s  saw  tremendous  refinements  to  the  physical  picture  of  the  photore¬ 
fractive  effect  and  the  development  of  a  powerful  predictive  theory,  Ninomiya  was 
the  first  to  consider  dynamic  processes  in  hologram  formation.  He  considered  the 
index  modulation  as  being  linear  with  the  exposure  (the  time-average  of  the  irradi- 
ance)  [95].  Magnusson  and  Gaylord  tested  Ninomlya’s  theory  experimentally  with 
the  result  that  they  qualitatively  verified  many  predictions  of  the  theory  (i.e.,  such 
as  rapidly  varying  writing  and  reading  results  which  previous  authors  had  taken 
as  experimental  problems”).  They  suggested  coupling  experiment  with  theory  to 
diagnose  material  parameters  of  recording  media  [87].  Vahey  used  an  approach 
similar  to  Ninomiya,  but  considered  the  index  modulation  as  being  proportional 
to  the  writing  irradiance  multiplied  by  an  arbitrary  time  dependent  function.  He 
also  succeeded  in  deriving  an  expression  for  diffraction  efficiency  (how  much  of  the 
reading  beam  is  diffracted  by  the  grating)  of  the  hologram  [127].  Kim  et  al.  consid¬ 
ered  feedback  between  the  space  charge  field  and  the  free  electron  density  and  thus 
were  pioneers  in  the  incorporation  of  materials  effects  in  the  dynamic  theory.  How¬ 
ever  they  did  not  consider  how  the  feedback  would  affect  the  writing  beams  [69]. 
This  oversight  was  included  in  the  work  of  Moharam  and  Young  who  also  included 
the  photovoltaic  effect  [93],  an  internal  electric  field  which  can  drive  an  intensity- 
dependent  photocurrent  [41].  1979  saw  the  publication  of  the  two  most  successful 
theories  to  date,  theories  which  looked  closely  at  material  properties.  Feinberg  et 
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al.  proposed  the  hopping  model  in  which  charges  are  fixed  to  certain  sites  in  the 
dark,  but  under  illumination  can  move  to  neighboring  sites  with  an  irradiance  de¬ 
pendent  probability.  This  approach  had  very  good  success  in  describing  effects  seen 
using  BaTi03[36].  The  most  general  theory,  and  the  one  most  favoured  today,  is 
from  Kukhtarev  et  al.  [77,76].  Chapter  2  will  review  the  important  equations  and 
provide  an  example  of  their  theory. 


1.3  Applications  of  Photorefractive  Materials 

It  would  be  impractical  to  present  an  exhaustive  list  of  proposed  applications  for 
photorefractive  materials  here.  As  Eichler  noted  in  his  introduction  to  the  spe¬ 
cial  issue  on  “Dynamic  Gratings  and  Four- Wave  Mixing”  in  the  IEEE  Journal  of 
Quantum  Electronics  in  1986,  “. . .  more  than  500  papers  have  been  published  up 
to  1985”  [29].  Our  time  would  be  better  spent  cataloguing  areas  of  suggested  uses 
with  references  to  one  or  two  representative  papers. 

The  first  division  of  applications  is  that  of  static  (storage  or  display)  versus  dy¬ 
namic  (real-time)  holography.  Certainly  the  first  proposed  use  for  a  photorefractive 
material,  (LiNb03),  was  as  a  medium  in  which  to  store  holograms.  A  reasonable 
question  to  ask  is  how  equal  are  they  to  this  task?  The  answer,  as  it  so  often 
is  in  science,  is  “that  depends”.  It  depends  on  what  type  of  crystal  we  are  dis¬ 
cussing.  Some  crystals,  such  as  BiuSiOjo  and  BinGeOjo?  are  very  sensitive  to  light 
and  are  very  quick  (grating  writing  times  for  moderate  («  10  mW/cm2)  intensities 
are  on  the  order  of  milliseconds).  They  compare  favorably  with  high  resolution 
photographic  plates,  but  suffer  from  short  storage  times  (20-30  hours)  [54].  Other 
crystals,  like  LiNb03  and  LiTa03,  respond  in  seconds  to  minutes,  but  have  much 
longer  storage  times  (1-10  years)  [44].  The  benefit  with  photorefractives,  of  course, 
is  that  you  have  a  self- developing,  erasable  medium.  However  it  is  possible  to  “fix” 
the  holograms  in  these  crystals  by  writing  the  hologram  at  a  high  temperature  (160° 
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C)  and  reading  it  at  room  temperature.  Staebler,  et  al.  estimate  that  the  result¬ 
ing  storage  time  for  a  fixed  hologram  in  LiNbOs  is  105  years.  Because  these  are 
volume  holograms,  the  angular  selectivity  is  very  strong — one  must  be  very  near 
to  Bragg  conditions  for  high  efficiency  [115],  one  can  store  many  holograms  in  the 
same  crystal  by  slightly  altering  the  angle  with  each  exposure.  In  one  experiment  a 
whopping  500  holograms  were  recorded  in  a  single,  2  mm  thick  sample  of  LiNbOa- 
Each  hologram  had  a  diffraction  efficiency  of  2.5%  or  better  [120].  As  promising 
as  these  results  sound,  very  few  papers  appear  nowadays  on  holographic  storage  in 
photorefractives — perhaps  the  feeling  is  that  all  the  work  possible  on  this  subject 
has  been  done,  or  the  cost  of  materials  (a  1  cm3  cube  of  Bii2Si02o  sells,  typically, 
for  £500)  is  too  high. 

In  dynamic  holography  we  take  advantage  of  the  ability  to  simultaneously  write 
and  read  a  hologram  in  a  photorefractive  material.  A  useful  way  to  characterise 
dynamic  applications  is  by  the  number  of  interacting  waves. 

In  wave  mixing,  as  in  an  argument,  the  minimum  number  required  to  produce 
interference  is  two.  The  simplest  case  is  when  the  writing  beams  are  uniform  plane 
waves  and,  as  we  have  already  discussed,  this  creates  a  sinusoidal  diffraction  grat¬ 
ing.  The  process  works  just  as  well  if  one  of  the  beams  has  temporal  or  spatial 
distributions,  and  two-wave  mixing  has  been  used  to  generate  real-time  interfero- 
grams.  An  interferogram  is  a  means  of  measuring  the  relative  spatial  displacement 
between  two  surfaces.  Departures  from  parallelism  result  in  bright  and  dark  fringes: 
dark  where  the  separation  is  a  half-odd  multiple  of  the  wavelength,  corresponding 
to  destructive  interference;  and  bright  where  the  interference  is  constructive.  You 
can  use  the  same  surface  if  you  superimpose  two  images  taken  at  different  times, 
which  is  particularly  useful  if  your  object  of  interest  is  vibrating.  As  mentioned 
above  the  response  time  in  Bii2SiO20  can  be  on  the  order  of  milliseconds.  If  the 
period  of  the  vibrating  object  is  much  less  than  this  response  time,  the  mixture  of 
the  object  beam  and  reference  beam  produces  a  continuous  interferometric  pattern 
[54,67].  In  LiNb03,  where  the  response  time  is  much  longer,  it  is  more  convenient 
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to  use  the  crystal  as  a  storage  hologram  and  make  a  double  exposure  interferogram 

(88!. 

Two- wave  mixing  in  photorefractives  generally  involves  power  transfer.  In  Chap¬ 
ter  2  we  will  discuss  this  phenomenon  in  great  detail,  but  we  should  mention  here 
that  in  Bi12Si02o  this  transfer  can  be  dramatically  increased  if  there  is  a  slight  dif¬ 
ference  in  frequency  between  the  waves  (detuning  typically  of  a  few  hundred  Hertz). 
Thus  we  can  use  two-wave  mixing  to  amplify  weak  signals.  In  Bii2Si020  gains  of 
103  and  in  BaTi03  gains  of  4xl03  have  been  reported  [109,126].  This  suggests  using 
photorefractive  materials  for  laser  beam  steering.  In  this  scheme  a  crystal  is  illumi¬ 
nated  with  a  strong,  uniform  pump  beam.  A  weak  probe  beam  (frequency  shifted 
if  necessary)  is  turned  on  in  the  required  direction.  The  large  amplification  of  this 
beam  essentially  redirects  part  of  the  pump  beam  in  the  desired  orientation  [108]. 

We  have  discussed  how  to  form  a  holographic  grating  with  an  object  and  a  ref¬ 
erence  wave.  A  third  wave  can  be  diffracted  off  this  grating  and,  if  it  is  sufficiently 
intense  and  near  the  same  wavelength,  may  interact  and  contribute  to  the  dynamic 
diffraction  process.  Of  particular  interest  is  when  the  third  beam  is  counterpropa- 
gating  to  the  reference  beam.  Then  part  of  the  third  beam  will  be  diffracted  from 
the  original  grating  from  behind.  This  diffracted  beam  runs  counter  to,  and  retraces 
the  path  of,  the  object  beam.  Its  phase  is  the  reverse  of  the  object  beam.  This 
latter  property  gives  the  diffracted  beam  its  name — a  phase  conjugate.  It  has  been 
likened  to  time  reversal  [133]. 

Although  we  started  with  three  beams,  we  ended  up  with  four,  so  this  scheme 
falls  under  the  category  of  four- wave  mixing.  The  arrangement  described  above 
which  produces  a  phase  conjugate  wave  is  often  labelled  a  phase  conjugate  mirror 
(PCM).  These  mirrors  are  characterised  by  a  threshold  [70].  Any  sufficiently  intense 
beam  incident  on  the  PCM  produces  a  reflected  conjugate  beam.  (In  analogy  with 
a  classical  mirror,  the  conjugate  is  the  “reflected  wave”.  However,  considering  the 
conjugate’s  turn  tail  character,  perhaps,  the  term  “retreating  wave”  is  more  apt.) 
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It  is  possible  for  a  single  input  beam  to  generate  a  phase  conjugate.  The  pho- 
torefractive  crystal  can  act  as  the  gain  medium  in  a  linear  resonator.  Optical  noise 
from  the  input  can  be  captured  and  amplified  in  the  resonator,  creating  the  required 
counterpropagating  beams  [131].  An  even  more  intriguing  arrangement  is  the  CAT 
mirror  (named  for  the  discoverer’s  pet  feline,  whose  portrait  ignominiously  was 
aberrated  for  the  experiment)  using  BaTiOa,  in  which  internal  reflections  provide 
all  the  necessary  beams  [37]. 

Phase  conjugate  mirrors  play  important  roles  in  a  wide  variety  of  applications. 
We  can  use  their  phase  reversing  properties  to  clean  up  beams  from  unwanted  ab- 
berations  [21]  or  compress  frequency  dispersed  pulses  in  optical  fibres  [99].  Since 
the  amplitude  of  the  conjugate  wave  is  related  to  the  product  of  the  amplitudes 
of  the  other  three  mixing  waves,  PCMs  are  useful  for  image  processing  tasks  such 
as  image  division  [60],  image  subtraction  [80],  intensity  inversion  [132]  or  edge  en¬ 
hancement  [35].  The  thresholding  property  of  PCMs  can  be  used  for  correlation 
[40]  and  associative  memories  [134]. 

The  response  time  of  photorefractives  is  proportional  to  the  incident  intensity. 
Very  intense  pulses  can  be  used  to  write  gratings  on  nano-  or  even  pico-second 
time  scales  [83,114].  Thus  some  of  the  applications  discussed  in  this  section  may 
be  possible  on  ultra-fast  time  scales.  The  limitation  is  the  ability  of  the  crystal 
to  dissipate  the  power.  For  pulses  of  gigawatt  instantaneous  power,  the  repetition 
frequency  must  be  no  higher  than  a  few  Hertz  [138]. 

1.4  Outline  of  the  Thesis 

Our  exploration  of  photorefractivity  begins  in  Chapter  2  with  a  discussion  of  the 
theory  describing  two-wave  interactions.  This  comparatively  simple  case  is  the 
springboard  to  more  complicated  geometries,  and  the  recipe  we  develop  to  handle 
it  will  serve  us  well  in  later  chapters. 

Two- wave  mixing  in  photorefractive  crystals  has  been  studied  for  nearly  twenty 
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years  and  one  might  think  that  there  are  no  new  wrinkles  to  the  problem.  We  can 
classify  these  interactions  in  terms  of  the  grating  spacing  of  the  mixing  beams  (a 
small  interbeam  angle  means  a  large  grating  spacing  and  vice  versa).  Certainly 
the  range  of  grating  spacings  from  submicrons  up  to  100  microns  has  been  studied 
(see  for  example  J.M.  Heaton’s  D.Phil.  thesis  [49]).  However,  in  Chapter  3  we 
bring  up  important  considerations  for  the  case  of  counterpropagating  beams  (grating 
spacing  equals  one  half  wavelength  —  the  smallest  possible  length),  and  explore 
nonreciprocal  effects. 

With  Chapter  4  we  lean  towards  the  other  extreme  of  the  grating  spacing  scale. 
Significant  gain  is  possible  in  Bi^SiOjo  for  large  (10-50  microns)  grating  spacings 
with  moving  intensity  fringes.  We  use  this  effect  in  a  unidirectional  ring  resonator. 
At  these  large  grating  spacings,  conditions  are  favourable  for  energy  transfer  to 
higher  diffracted  orders. 

Chapter  5  describes  theory  and  experiments  involving  a  new  mechanism  for 
gain  in  photorefractive  materials.  This  amplification  depends  on  off-Bragg  coupling 
and  is  more  efficient  in  many  cases  than  our  stalwart  gain  mechanism  derived  in 
Chapter  2. 

The  production  of  higher  diffracted  orders  is  investigated  in  Chapter  6  using  a 
novel  techninque-we  take  advantage  of  a  very  slight  wedge  angle  in  a  crystal  sample. 
Two-wave  mixing  is  provided  by  a  single  incident  beam  and  its  nearly  coincident 
internal  reflection. 

We  finish  up  with  conclusions,  suggestions  for  future  work,  and  some  relevant 
appendices. 


T 


Chapter  2 

Two- Wave  Mixing  in 
Photorefractive  Crystals 


“ The  revelation  came  to  me  while  doing  homework  on  partial  differential  equations. 
I  was  a  junior  and  I  had  to  calculate  the  torsional  modes  of  vibration  in  a  metal  beam. 
That  was  when  I  realized  that  solving  a  differential  equation  gives  you  a  feeling  of 
accomplishment,  but  not  one  of  justice.  ” 

“ Then  you  were  doing  it  wrong,  son.  ” 

“The  Paradoxical  Case  of  the  Relativity  Twins” 


To  create  a  dynamic  holographic  grating  in  a  photorefractive  material  three 
processes  must  occur. 

1.  Transport  of  photoinduced  carriers. 

2.  Modulation  of  the  refractive  index  via  the  linear  electrooptic  effect. 

3.  Diffraction  of  the  writing  beams  as  the  grating  forms. 

The  first  process  can  be  described  in  terms  of  materials  equations;  the  second 
is  worked  out  by  considering  the  electrooptic  tensor;  and  the  third  is  specified  by 
field  equations.  In  this  chapter  we  adopt  the  strategy  of  divide  and  conquer — each 
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process  will  handled  separately.  This  makes  the  problem  much  more  tenable  and 
allows  the  luxury  of  pausing  occasionally  to  discuss  the  physics  involved.  By  the  end 
of  the  chapter  we  will  have  derived  an  expression  for  power  transfer  for  two-wave 
mixing  and  we  will  consider  means  to  control  and  optimize  this  effect. 

Two-wave  mixing  is  fundamental  to  understanding  the  whys,  hows  and  where¬ 
fores  of  the  photorefractive  effect.  The  interaction  of  several  waves  often  can  be 
reduced  to  a  linear  combination  of  mixing  of  all  the  possible  pairs  of  waves. 


2.1  The  Materials  Equations 

One  of  the  most  quoted  papers  dealing  with  photorefractivity  appeared  in  two 
parts  in  Ferroelectric s  in  1979.  N.  V.  Kukhtarev  and  his  coworkers  [77,76]  laid 
down  a  consistent  framework  of  materials  equations  describing  the  generation  and 
recombination  of  charge  carriers,  their  transport  in  the  form  of  a  current  density, 
and  the  establishment  of  a  space  charge  field. 

2.1.1  The  rate  equations 

In  this  band  transport  model  charge  carriers  are  liberated  from  fixed  donor  sites  in 
the  bandgap.  Thus  one  way  to  monitor  the  carrier  concentration,  n,  is  to  consider 
the  ionized  donor  density,  Np  (here  we  assume  an  n-type  semiconductor).  (There 
are  also  negatively  charged  traps,  whose  concentration  is  given  by  TV ~A — this  will 
appear  in  a  later  equation.)  These  ionized  donors  act  as  traps  for  unwary  carriers 
and  thus  both  n  and  Np  can  be  depleted. 

The  first  materials  equation  describes  the  time  rate  change  of  Np. 

=  (sI  +  /3)(Np-N+)  -7rnTV+  (1) 

The  first  term  on  the  right  hand  side  of  the  equation  is  the  creation  term.  There 
are  two  avenues  of  increasing  Np:  photo-  or  thermal  generation,  s  is  the  photoion¬ 
ization  constant,  I  is  the  driving  light  intensity  and  0  is  the  thermal  generation 
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coefficient.  Note  that  the  production  of  ionized  donors  is  proportional  to  the  den¬ 
sity  of  un-ionized  donors  (Np  —  -Nj). 

The  second  term  describes  recombination.  ~)T  is  the  recombination  coefficient, 
and  we  see  the  depletion  of  is  proportional  to  the  product  of  the  carrier  and 
trap  densities.  This  means  the  recombination  term  is  nonlinear. 

The  donor  sites  are  fixed.  Whenever  a  photon  or  a  phonon  begets  an  ionized 

donor  a  charge  carrier  is  produced.  But  the  charge  carriers  are  mobile —  there  is  an 

additional  process  by  which  the  carrier  density  can  change.  That  is  via  a  current 

density.  Our  second  materials  equations  is: 

dn  _  dN '£  _  1 
dt  dt  e 

This  equation  should  look  familiar  to  anyone  well  versed  in  electromagnetic  theory. 
It  is  the  continuity  equation,  where  e  is  the  unit  electric  charge  and  J  is  the  current 
density.  It  is  a  statement  of  the  conservation  of  charge. 


2.1.2  The  current  density 

The  flow  of  a  current  requires  some  agent  to  move  the  charges.  Our  third  materials 
equation  considers  three  possible  processes: 

J  =  penE,  4-  pi  [Np  —  Np')  -  pkgTXn  (3) 

The  first  two  terms  on  the  right  hand  side  of  the  equation  are  drift  current  den¬ 
sities.  The  first  is  due  to  a  space  charge  electric  field,  E,,  which  we  allow  to  include 
any  background  field  due  to  an  applied  voltage.  Here  p  is  the  carrier  mobility. 
The  second  term  is  a  photovoltaic  drift  current  density  which  is  proportional  to  the 
intensity  and  the  density  of  un-ionized  donors.  The  constant  of  proportionality  is 
the  photovoltaic  coeffcient,  p. 

The  final  term  is  the  diffusion  current  density,  which  depends  on  the  gradient  of 
the  carrier  concentration.  Here  kg  is  Boltzmann’s  constant  and  T  is  the  absolute 
temperature. 
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Cona  uct  io  n 


(  s  I  +  p  )  (  Nd  -  Np  ) 


donor 


trap 


valence 

Figure  1:  Band  transport  model:  Generation,  transport  and  trapping  of  n-type 
charge  carriers.  Donors  and  traps  lie  in  the  bandgap  between  the  conduction  and 
valence  bands. 

Fig.  1  summarizes  the  first  three  materials  equations.  In  photorefractive  ferro¬ 
electric  crystals,  such  as  LiNb03  and  BaTi03,  the  donor  and  trap  sites  are  impuri¬ 
ties,  mainly  transition  metals  such  as  Fe,  Cu  and  Mn  [72,71].  In  Fe  doped  LiNb03, 
for  example,  the  charge  transport  occurs  between  Fe2+  and  Fe3+  [73].  The  donor 
sites  in  Bi^SiOjo  crystals  r. re  silicon  vacancies  that  lie  2.6  eV  below  the  conduction 
band  (the  band  gap  energy  is  3.25  eV  [52]).  The  nature  of  the  traps  in  BiuSiOjo  is 
not  well  understood,  but  they  lie  0.6  eV  below  the  conduction  band.  Luminescence 
centres  (2.25  eV  below  the  band  when  occupied  by  an  electron,  1.3  eV  when  not 
occupied)  may  also  contribute  to  charge  transport  [44]. 

2.1.3  The  space  charge  field 

Our  last  equation  is  Poisson’s  equation.  The  divergence  of  the  space  charge  field  is 
related  to  the  charge  density. 


V  •  E.  =  -  (NX  +  n  -  N+) 

e.  v  ' 


(4) 
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e,  is  the  permittivity  of  the  medium  and  is  an  ionized  acceptor  density.  The 
presence  of  ionized  acceptors  is  necessary  to  insure  charge  neutrality  because  we 
will  assume,  for  most  practical  cases,  n.  That  is,  in  a  Fourier  expansion  of 

Nq  and  n,  the  coefficients  obey  3>  n,  where  i  is  an  index.  Normally  we  take 
this  expansion  out  to  first  order. 

2.1.4  Simplifying  assumptions 

In  theory  we  could  apply  any  driving  intensity,  /,  to  Eqns  1-4  and  find  the  resulting 
space  charge  electric  field,  E,.  (That  is  our  goal,  for  this  is  the  field  which  will  mod¬ 
ulate  the  refractive  index  through  the  electrooptic  effect.)  In  practice,  though,  our 
first  step  is  to  reduce  the  complexity  of  the  equations  by  making  some  simplifying 
assumptions.  Here  we  list  each  assumption  and  briefly  discuss  its  validity. 

1.  There  is  only  one  species  of  charge  carriers.  Implicit  in  the  statement  of  the 
Kukhtarev  model  is  that  only  one  type  of  charge  carrier,  electrons  or  holes,  is  in¬ 
volved  in  the  transport  process.  Orlowski  and  Kratzig  have  demonstrated  that  the 
transport  of  both  holes  and  electrons  is  deleterious  to  photorefractivity  [100].  How¬ 
ever,  simultaneous  electron/hole  transport  does  occur  in  Fe-doped  LiNbO3[130)  and 
BaTi03(71J.  The  acceptor/donor  concentration  ratio  (Fe2+/  Fe3+)  can  be  controlled 
during  the  crystal  growth  process,  and  this  ratio  determines,  to  a  large  extent,  the 
relative  contributions  of  the  electrons  and  holes  to  photoconductivity.  Thus  by  judi¬ 
cious  doping  one  can  regain  single  charge  carrier  status.  Valley  [128]  and  Strohkendl, 
et  al.  (123)  independently  have  developed  charge  transport  models  for  the  presence 
of  both  electrons  and  holes.  These  models  are  based  on  Kukhtarev’s  theory  and 
contain  similar  expressions  of  Eqns  1-3  for  holes.  Valley  considered  two  cases — one 
in  which  the  recombination/donor  centres  are  the  same  for  holes  and  electrons,  and 
one  in  which  they  are  different.  In  Bi^SiOjo  electron  transport  dominates,  while  in 
BiuGeOjo  either  holes  or  electrons  govern  photoconductivity  [102]. 
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2.  Thermal  generation  of  carriers  is  negligible.  The  justification  for  this  assump¬ 
tion  comes  from  the  low  dark  conductivities  of  most  photorefractive  crystals  i44i. 
Our  own  measurements  of  dark  current  in  Bii2Si02o  (applied  voltage  but  no  illumi¬ 
nation)  gave  values  typically  in  the  picoamp  range,  while  with  moderate  incident 
intensities  (%  10  mW/cm2)  the  current  rose  to  microamps.  Thermal  generation  in 
Eqn  1  can  be  neglected  except  for  very  low  light  intensities. 

3.  The  carrier  density,  n,  is  negligible  compared  with  the  ionized  donor  density, 
Nq.  Several  authors  have  measured  trap  densities  in  photorefractive  crystals.  In 
LiNb03[16],  BaTi03[36),  and  Bix 2Si02o[94]  the  densities  are  on  the  order  of  magni¬ 
tude  of  1022  m“3.  Carrier  densities  are  usually  much  lower  than  this.  For  example, 
a  uniform  illumination  of  100  mW/cm2  in  Bi12SiO20  produces  a  carrier  density  of 
%  2  x  1017  m~3.  Thus  N£  -  n  %  N£. 

4.  The  ionized  donor  density,  N^,  is  negligible  compared  with  the  donor  density, 
No-  Here  again  measurements  provide  the  justification.  In  Bi12SiO20,  ND  is  on  the 
order  of  1026  m~3  [109],  while  in  LiNb03  it  is  1023  m~3  [16].  Therefore  ND  -  N£  ss 

nd. 

5.  The  ionized  acceptor  density,  N is  constant.  Ionized  acceptors  must  be 
present  to  offset  the  initial  ionized  donor  density  to  insure  charge  neutrality.  We 
assume  these  acceptors  play  no  active  role  in  the  charge  transport  process. 

2.1.5  The  driving  light  intensity 

Consider  the  intensity  pattern  formed  by  the  superposition  of  two  waves  of  nearly 
the  same  wavelength  as  shown  in  Fig.  2a.  This  is  referred  to  as  nearly  degenerate 
two  wave  mixing.  The  grating  vector  (the  spatial  angular  frequency  of  the  intensity 
pattern)  is  seen  in  Fig.  2b. 

I(r)  ~  | E+ej{“+t-k+-f)  +  E- =  /„  +  0  +  c.c.  (5) 

where  /0  =  /++/_  ~  |.£'+|2  4-  \E-\2  is  the  background  intensity,  7)  =  %//+  I-  is 
the  modulated  intensity,  k+/_  =  2tt/A+/_  is  the  +/—  wave  vector,  with  A+/_  the 
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(a)  (b) 


Figure  2:  Two- wave  mixing,  (a)  The  waves  inside  the  medium,  (b)  The  grating 
vector  in  reciprocal  space. 

wavelength  in  the  medium,  \K\  —  |fc+  —  fc_|  =  2n/A  is  length  of  the  grating  vector 
with  A  the  grating  spacing,  6u  =  u/+  -u/_  is  the  frequency  detuning,  and  c.c.  means 
complex  conjugate.  For  6ui  ^  0  this  intensity  pattern  consists  of  sinusoidal  fringes 
moving  with  velocity  v  =  SwjK.  It  is  common  practice  to  set  up  a  coordinate 
system  in  which  K  is  parallel  to  one  of  the  axes,  say  z,  so  K  ■  f  =  K  z.  This  total 
light  intensity  will  drive  the  materials  equations. 

At  this  point  it  is  convenient  to  introduce  normalized  quantities  for  n,  Np, 
E,,  and  J .  Experience  shows  that  this  will  save  a  lot  of  algebra  later.  With  our 
simplifying  assumptions  our  normalized  materials  equations  may  be  written: 
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=  (l  +  mejK{vt'z>  +•  c.c.)  -  —  nn;Vn  (6) 

=  •  /„  (7) 

-  (8) 

=  1  -  tfn  (9) 

where  our  normalized  quantities  are  iVn  =  N^/N^,  rin  —n/N^ ,  l£n  =  EJEq ,  Jn  = 
J/ HeN^  Eq.  The  constants  are  h  =  / sIqNd,  m  =  A /Zo»  ^  =  /i,/rr,  rr  = 

1  /yrN^ .  The  characteristic  electric  fields  are  2T,  =  eN^/e.K,  Et  —  kgTK/e ,  2?p  = 
PlrN^/sen- 

The  next  step  is  to  assume  the  solutions  to  Eqns  6-9  are  driven  by  I  and  thus 
have  the  same  general  functional  form. 

Fn(z,t)  =  Fn 0  +  FnlejK^  +  c.c.  (10) 

Note  we  assume  steady-state  conditions  for  the  coefficients  Fn 0  and  Fnl .  These  are 
the  first  two  coefficients  in  a  spatial  Fourier  series.  We  will  take  the  directions  of 
En,  Ep  and  /„  to  be  i,  i.e.  parallel  to  K . 

There  is  one  more  obstacle  to  overcome.  Eqns  6  and  8  are  nonlinear  (they  contain 
the  terms  nnNn  and  nnE„).  To  overcome  this  difficulty  we  invoke  linearization.  We 
assume  the  amplitude  of  the  first  order  term  is  much  smaller  than  the  amplitude 
of  the  zero  order  term.  This  is  a  reasonable  assumption  if  the  intensity  modulation 
is  very  small  (m  «  1).  Then  we  may  ignore  products  of  first  order  terms.  Making 
the  substitutions  of  Eqn  10  and  applying  linearization,  we  next  equate  coefficients 
of  the  static  (d.c.)  and  exponential  (a.c.)  terms.  This  gives  us  eight  equations  and 
eight  unknowns. 


dNn 
1  dt 
dN„ 
dt 

fn 


K 


VEn 
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2.1.6  The  d.c.  equations 


^nO  —  ^na 

(11) 

2* 

o 

II 

1 

(12) 

Jnd  —  ^nO^nO  H"  ^na  *-» 

(13) 

« 1  r  5> 

II 

O 

J 

(14) 

Here  E, 0  is  the  static  electric  field  due  to  any  applied  voltage. 

2.1.7  The  a.c.  equations. 


m 

=  — «nl  +  (1  - 

^na 

jKvh)Nnl 

(15) 

Nnl 

1 

II 

(16) 

Jr* 

/  E  \  E 

”  ^na^nl  “I-  j ~jjj~  J 

(17) 

En  I 

ii 

(18) 

2.1.8  The  modulated  space  charge  field 

Our  goal  was  to  find  the  modulated  space  charge  field,  Etl,  and  we  now  have  the 
means.  The  following  method  works  well.  From  the  a.c.  equations  eliminate  Jnl, 
then  n„i,  and  lastly  Nni.  The  final  result  is: 


E,  i  =  -mEw 

^  _  _  E, o  +  Ep  —  j  Et _ 


(19) 

(20) 


where  Em  =  yrNA/nK,  b  =  Kvtj  and  =  e,yrNA  /efisIoNo  is  the  dielectric 
relaxation  time. 
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We  will  defer  a  full  discussion  of  the  physical  meaning  of  Eqn  20  until  later  in 
the  chapter  when  we  have  worked  out  the  rest  of  the  photorefractive  process.  We 
should  notice,  however,  that  E,i  is  directly  proportional  to  the  intensity  modula¬ 
tion,  m.  and  it  is  complex.  This  means  the  field  can  be  out  of  phase  with  the 
intensity  pattern.  This  is  what  we  dubbed  the  phase  mismatch  angle,  4>9,  in  Fig.  2 
in  chapter  1. 


2.2  The  Linear  Electrooptic  Effect 

The  next  question  to  ask  is  what  influence  the  space  charge  field.  E ,,  will  have  on  the 
photorefractive  medium.  The  answer  comes  from  considering  the  linear  electrooptic 
effect.  In  this  effect  geometry  and  orientation  are  crucial. 

In  many  problems  we  have  the  luxury  of  choosing  a  convenient  coordinate  sys¬ 
tem.  For  example,  in  the  previous  section  we  chose  the  z  axis  as  being  parallel  to 
the  grating  vector,  K.  We  are  much  more  restricted  when  we  consider  propagation 
of  light  in  an  anisotropic  medium,  such  as  a  crystal.  Optically  speaking  there  is 
only  one  preferred  set  of  axes — the  principal  coordinate  system. 

Simply  put,  the  principal  coordinate  system  is  the  one  in  which  the  dielectric 
tensor,  e,  is  diagonal.  When  an  electric  field,  E,  is  applied  to  a  dielectric  the  medium 
is  polarized.  The  electric  displacement,  D,  is  related  to  this  field  by: 

D  =  lE  (21) 


In  the  principal  coordinate  system,  e;9  =  0  for  i  -fc  j.  Along  the  axes,  D  is  in  the 
same  direction  as  E  [46). 

If  the  medium  is  homogeneous,  nonabsorbing  and  magnetically  isotropic  then 
the  electric  energy  density  is  given  by: 


we=12e-d 


(22) 
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Propagation  direction 


i 


Cross  section 


Figure  3:  (a)  The  index  ellipsoid  with  propagation  direction,  s',  (b)  Propagation 
eigenstates  from  the  cross-section  ellipse. 


Making  a  change  of  variable,  Xj  =  D./v/^Wgeo,  and  recalling  the  refractive  index 
=  y/to Ti  gives: 

'  <23> 

i  =  l  ni 

The  above  is  the  equation  defining  the  index  ellipsoid.  The  x*s  are  the  principal 
axes  coordinates.  Thus  in  the  principal  coordinate  system  the  index  ellipsoid  has  a 
particularly  simple  form. 


The  index  ellipsoid  is  used  to  describe  propagation  of  light  in  anisotropic  media. 
For  any  propagation  direction  there  are  generally  two  eigenstates,  each  characterised 
by  a  polarisation  and  a  phase  velocity  (refractive  index).  By  eigenstate  we  mean 
an  incident  wave  polarised  in  one  of  these  characteristic  directions  will  maintain 
its  polarisation  state  as  it  progresses  through  the  medium.  Fig  3  shows  how  to 
find  these  eigenstates.  Construct  a  plane  passing  through  the  origin  of  the  index 
ellipsoid  whose  normal  is  in  the  desired  propagation  direction.  The  intersection  of 
this  plane  and  the  ellipsoid  is  an  ellipse.  The  directions  of  the  major  and  minor  axes 
give  the  eigenpolarisations  while  the  lengths  of  the  semiaxes  are  the  corresponding 
refractive  indices. 
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The  treatment  for  the  index  ellipsoid  concerned  the  optical  electric  field.  An 
applied  field,  such  as  the  space  charge  field  developed  in  Sec  1,  may  perturb  the 
distribution  of  charges  in  the  medium  and  alter  the  index  ellipsoid.  This  is  known 
as  electrooptic  effect  and  the  perturbation  may  be  proportional  to  the  applied  field 
(linear  or  Pockel’s  effect)  or  to  the  square  of  the  applied  field  (quadratic  or  Kerr 
effect). 

In  photorefractivity  it  is  the  linear  electrooptic  effect  which  is  important.  Using 
the  contracted  index  notation  of  Yariv  and  Yeh  [135]  an  applied  field  transforms 
the  index  ellipsoid  into: 

+  r'jE])  X'  +  +  r2jEj )  X*  +  ~  7"3j'£:j)  **  + 

2r^E,x2x3  4-  2r&jEjXlx3  +  2 r%JEjX^x2  =  1  (24) 

with  repeated  summation  over  j  =  1,2,3.  The  t\;s  are  the  coefficients  of  the 
electrooptic  tensor.  The  form  of  this  tensor  depends  on  the  symmetry  of  the  crystal 
class  [30].  For  example,  several  crystal  classes  are  centrosymmetric;  that  is  they 
have  a  centre  of  symmetry  and  are  invariant  under  an  inversion  (replacing  the 
electrooptic  tensor  r  with  -r).  Performing  this  inversion  in  Eqn  24  gives  a  different 
form  for  the  index  ellipsoid  unless  all  the  r^s  are  zero.  Thus  a  necessary  (but  not 
sufficient)  condition  for  exhibition  of  the  linear  electrooptic  effect  is  that  a  crystal 
be  noncentrosymmetric.  Ref  [135]  contains  a  list  of  forms  for  the  electrooptic  tensor 
for  all  the  crystal  classes  as  well  as  the  individual  coefficients  for  several  crystals. 

Eqn  24  represents  a  transformation  of  the  original  index  ellipsoid.  In  this  trans¬ 
formation  the  lengths  of  the  ellipsoid  semi-major  and  -minor  axes  may  change  and 
their  orientation  may  suffer  a  rotation.  In  other  words  we  might  have  a  new  prin¬ 
cipal  coordinate  system  with  the  application  of  an  electric  field. 

Most  of  our  experimental  work  with  photorefractive  crystals  will  involve  Bi^SiOjo 
in  the  configuration  shown  in  Fig  4a.  The  Miller  index  notation  gives  us  the  orien¬ 
tation  of  the  principal  coordinate  system  with  no  applied  field.  Fig  4b  shows  the 
new  orientation  of  the  principal  axes  for  a  field  applied  along  the  001  or  x3  crystal 
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Figure  4:  a)  Principal  coordinate  system  (dark  axes)  for  Bii2SiO20.  b)  Reoriented 
principal  coordinate  system  for  applied  field  along  the  x3  axis. 

axis.  The  analysis  of  this  case  is  as  follows.  Bi12SiO20  belongs  to  crystal  class  23. 
It  is  a  cubic  crystal  so  with  no  applied  field  the  index  ellipsoid  is  a  sphere  of  radius 
n.  The  only  non-zero  elements  of  the  electrooptic  tensor  are  r41,r52,  and  r63,  all  of 
which  are  equal.  Assuming  an  electric  field  E  applied  along  the  principal  i3  axis, 
our  index  ellipsoid  is  transformed  into: 

— |  — I  +  -|  +  2rmExxx2  =  1  (25) 

TV  TV  TV 

In  the  new  principal  axes  system  the  index  ellipsoid  has  the  form  given  by  Eqn  23. 
The  appropriate  transformation  is: 

(26) 

=  -^  (-*,+*,)  (27) 

=  *3  (28) 


This  is  merely  a  rotation  about  the  x3  axis  through  45°.  In  the  new  primed  principal 
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coordinate  system  the  index  ellipsoid  is: 

(1  -  r63£)  x?  4-  (1  -  r63JE;)  x'22  -  ~  =  1  (29) 

The  quantity  r63£  1/n2  so  the  new  ellipsoid  axes  lengths  are 

rii  =  n+-n3r63E  (30) 

1  , 

n2  =  n  —  -n3r63E  (31) 

n3  -  n  (32) 


In  the  primed  principal  coordinate  system  light  polarised  in  the  x[  or  x2  direction 
will  experience  a  refractive  index  change,  but  light  polarised  along  the  i3  axis  will 
not. 

In  general  we  can  describe  the  change  in  refractive  index  in  terms  of  a  change 
in  the  dielectric  constant  (for  optical  frequencies)  given  below: 

Aer  =  e,lreff£  (33) 

where  reff  is  the  appropriate  electrooptic  tensor  component  or  combination  of  com¬ 
ponents. 


2.3  The  Field  Equations 

Consider  the  mixing  of  two  waves  in  a  photorefractive  medium  as  shown  in  Fig  5. 
How  will  the  modulated  refractive  index  found  in  Sec  2.2  affect  the  waves?  We 
need  to  find  the  expression  of  the  governing  wave  equation.  Here  we  will  neglect 
any  vectorial  complications  such  as  optical  activity  or  birefringence.  These  effects 
will  be  discussed  briefly  later  in  the  chapter. 

The  superposition  of  the  two  waves  is  given  by: 


E  =  E+{z)e~A^  *—+*)  + 


(34) 
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Figure  5:  Orientation  for  two-wave  mixing 
The  modulated  dielectric  constant  is: 

er  =  er0  +  +  Acr'e-X*  (35) 

Here  means  complex  conjugate. 

The  derivation  of  the  scalar  wave  equation  is  a  straightforward  exercise  and 
can  be  found,  for  example,  in  Ref  (115).  We  take  advantage  ol  the  fact  that  the 
angular  frequency  difference,  Su,  between  the  waves  is  very  small  compared  to  either 
angular  frequency.  Thus  u>+  %  u>_.  This  also  means  the  time  derivatives  of  e  are 
much  smaller  than  the  time  derivatives  of  E.  The  scalar  wave  equation  is  given  by: 

V2E  A^fiotrtoE  =  0  (36) 

where  uj  is  the  angular  frequency  of  the  two  waves,  /r0  is  the  permeability  of  free 
space,  eT  is  the  relative  dielectric  constant  of  the  medium,  e0  is  the  permittivity  of 
free  space,  and  E  is  the  total  complex  amplitude. 

Putting  Eqns  34  and  35  into  36  yields: 

(V2£+  -  2 jk+  ■  VE+  -  k2E+)  e *0  + 

(V2£_  -  2 jk_  •  V£_  -  k2E _)  + 

u>Wo  (er0  +  Aere^,f  r'-<“‘)  + 

(E+e-}{ii+f'-w+t)  4-  =  0 


(37) 
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K 

K 

K 

Figure  6:  Bragg  closure.  The  propagation  directions  corresponding  to  -  K  and 
fc_—  K  are  neglected. 

where  jfc+|  ~  |^-|  =  k. 

As  is  stands  the  above  equation  is  difficult  to  work  with-  It  is  time  for  some 
simplifying  assumptions. 

2.3.1  More  simplifying  assumptions 

1)  Neglect  terms  involving  the  second  derivative.  This  is  reasonable  because  E+ 
and  E-  vary  slowly  in  space.  It  is  a  common  assumption  and  is  known  as  the 
“slowly  varying  amplitude”  approximation  [44].  An  alternative  way  of  reaching  this 
conclusion  is  to  look  at  the  first  derivative  terms.  These  are  multiplied  by  fc  which 
is  a  very  large  number.  Thus  these  terms  dominate  the  second  derivatives. 

2)  Neglect  higher  diffracted  orders.  The  vectors  k+,  and  K  form  a  Bragg  clo¬ 
sure.  In  Eqn  37,  an-:  - 

c-i(*+ 0.  These  products  correspond  to  valid  propagation  directions  for  the  solu¬ 
tion  given  by  Eqn  34.  The  products  arlc[  e-j(H-  r-^-t)ej[R  r-iut) 

suggest  higher  diffracted  orders  which  do  not  fit  the  Bragg  closure.  In  a  volume 
grating  angular  selectivity  is  very  strong  and  any  departure  from  the  Bragg  condi¬ 
tion  greatly  reduces  the  diffraction  efficiency  [115],  This  assumption  will  be  valid 
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except  for  very  small  interbeam  angle?. 

3)  1  se  one  dimensional  theory.  Assume  the  incident  beams  make  very  smaii 
angles  with  the  x  axis  so  A’_  _  rr  ~  kx. 

With  these  three  assumptions  and  the  relationship 

k2  =  Vi,'_/Xot0tr0  (38) 


we  obtain 


dEj.  .Ac,*  k  ^ 

dx  J  2er0 
dE _  A  erk 


-7 -  =  (40) 

dx  4trQ 

(In  reflection  geometry  (see  Fig  2.7c)  the  right  hand  side  of  Eqn  40  is  positive  rather 
than  negative.) 

2.3.2  Two-wave  mixing  coupled  wave  equations 

Using  Eqn  33,  Eqns  39  and  40  can  be  written: 


dE+  I-  r 

IT  = 

17  -  jr‘”uTLE- 


where  T  =  jfc0n3rc^i£)a.|  with  k0 ,  the  free  space  wave  number,  and  <ES 
In  terms  of  the  intensity  and  phase  of  each  wave  we  may  write: 


=  arg(£’w). 


=  2r  sin 


-  -a2r  sin  <£ 


9/+  +  /_ 

.  ui- 


—  r  cos  $ 


9/+  +/. 
/_ 


—  af  cos  $ 


9  Em  E 
.  E 


E  +  E 


Here  a  =  -t-1  for  transmission  geometry  and  —1  for  reflection  geometry. 


CHAPTER  2:  TWO-  WAVE  MIXIXG 


33 


a)  Transmission  b)  Reflection  c|  Ccunterprccaaating 

geometry  geometry  geometry 

Figure  7:  Geometries  for  Two-Wave  Mixing,  (a)  Transmission,  (b)  Reflection,  (c) 
Counterpropagating. 

2.3.3  Consequences  of  the  coupled  wave  equations 

The  coupling  in  Eqns  43  -  46  is  highly  dependent  on  the  phase  mismatch  angle,  <I>9. 
If  ~  0  there  is  no  coupling  in  intensity  between  the  beams.  When  <!>,,=:  90°  the 
coupling  is  strongest  and  we  may  expect  maximum  power  transfer.  The  reverse  is 
true  for  the  phase  equations.  But  what  does  coupling  of  the  phase  mean  physically? 
It  means  the  phase  difference  between  the  beams  can  change  as  one  moves  across 
the  crystal.  This  means  the  fringes  of  the  steady-state  grating  are  bent  93] . 

Consider  the  symmetric  geometries  shown  in  Fig  7.  We  should  point  out  that  in 
general  the  crystal  orientations  differ  between  transmission  and  reflection  geometries 
to  take  best  advantage  of  the  electrooptic  effect.  The  counterpropagating  beam 
configuration  is  a  special  case  of  the  reflection  geometry  and  will  be  discussed  in 
detail  in  Chapter  3.  In  this  antiparallel  case  the  steady-state  grating  fringes  are  not 
bent,  but  chirped  as  a  result  of  the  phase  difference  in  the  beams. 

Power  conservation  is  different  for  transmission  and  reflection  geometries.  For 
the  transmission  case  the  sum  /+  +  /_  is  constant  while  it  is  easily  shown  that  the 
reflection  case  gives  constant  1+  —  /_.  We  can  use  power  conservation  to  obtain 
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solutions  for  the  coupled  wave  equations.  For  the  transmission  case,  substitute 
/j.  =  /0  —  I-  into  the  differential  equations  and  solve. 


M*) 

/_(*) 

4>+{x) 

4>-{x) 


3r  exp(2T  sin  $gx) 

1  4-  (3t  exp(2Fsin  <I>s.r)  ° 

1  r 

1  + /3r  exp(2r  sin  ° 

<£+(0)  +  T  cos  $gx  t -  -  cot  $g  In 


_ 1  +  (3r _ 

1  +  0T  exp(2T  sin  $gx) 


<M°)  -  -  cot  $g  In 


_ 1  +  3r _ 

\  +  0T  exp(2 r  sin  $9x) 


Here  /3r  =  /*(0)//_(0),  the  input  beam  ratio. 


(47) 

(48) 

(49) 

(50) 


2.3.4  Crystal  parameters 

Different  crystal  types  generally  have  different  photorefractive  responses.  There  may 
even  be  variation  between  crystals  of  the  same  type  but  from  different  manufacturers 
[31].  The  presence  of  impurities  or  dopants  may  lead  to  variation  in  effects  such  as 
absorption,  relative  concentrations  of  acceptors  and  donors,  or  the  strength  of  the 
photovoltaic  field.  Ultimately  the  individual  behaviour  of  a  crystal  sample  is  due 
to  its  distinct  combination  of  crystal  parameters. 

Table  1  lists  typical  values  from  the  literature  for  the  parameters  which  make 
up  E,i  for  Bii2Si02oi  LiNb03,  and  BaTi03.  These  are  the  crystals  we  will  consider 
theoretically  and  experimentally  in  the  rest  of  the  thesis.  LiNb03  and  BaTi03  are 
uniaxial  crystals  and  have  ordinary  (n„)  and  extraordinary  (ne)  refractive  indicies. 
There  are  two  important  electrooptic  tensor  components  for  the  applied  electric  field 
and  grating  vector  direction  along  the  001  axis.  In  LiNb03  r33  is  the  appropriate 
value  for  transmission  geometry  and  ri3  for  reflection  geometry.  We  have  taken  p 
as  l.OxlO-30  m3/V  to  give  a  photovoltaic  field  of  2.5  kV/cm,  a  value  we  measured 
for  one  sample  of  LiNb03.  The  photovoltaic  effect  is  negligible  in  BaTi03  and 
Bij2Si02o.  BaTi03  has  a  very  large  electrooptic  tensor  component,  r51,  which  en¬ 
hances  two- wave  mixing  in  transmission  geometry.  The  r]3  component  is  important 
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Parameter 

BiijSiOjo 

LiNb03 

BaTi03 

ND(m~3) 

l.OxlO25 

l.OxlO23 

l.OxlO25  1 

> 

j 

? 

l.OxlO22 

l.OxlO22 

2.0xl022 

e. 

56 

fn  =  78 
£33  —  32 

e„  =  4300  j 
£33  =  168 

n  (refractive  index) 

2.62 

na  =  2.29 
nt  =  2.20 

n0  =  2.44 
ne  -  2.37 

reff  (m/y) 

3.4xl0-12 

r13  =  9.6xl0-12 
r33  =  3.1xlO-!I 

r13  —  2.6xl0-n 
r5i  =  1.6X10-9 

H  (m2/'Vs) 

3.0xl0~6 

l.OxlO'5 

5x1 6 -B 

s  (m2/J) 

2.0xl0-5 

2.5xl0-5 

3xl0-6 

p  (m3/ V) 

0 

l.OxlO-30 

0 

7,(m3/s) 

1.6xl0-17 

l.OxlO-15 

5xl0-14 

Table  1:  Crystal  parameters  for  Bit2Si02o,  LiNb03,  and  BaTi03  for  A  =  514  nm. 
Refs  [49,135,130,16,86] 

in  reflection  geometry.  The  values  listed  are  for  a  wavelength  of  514  nm. 


2.3.5  Optimization  of  intensity  coupling 

The  characteristic  fields  are  Et,  E ,,  Em,  and  Ep.  Et  is  a  diffusion  field.  E ,  is 
proportional  to  the  maximum  value  the  space  charge  field  can  have.  It  is  1  /27T 
times  the  field  due  to  the  separation  of  ionized  donors  and  acceptors  by  one  grating 
spacing  [117].  Em  is  a  consideration  during  the  transient  period  while  the  grating 
forms  [49]  or  when  moving  fringes  are  present,  i.e.  when  the  two  mixing  beams  have 
different  frequencies  [121].  The  photovoltaic  field  is  Ep.  It  is  a  constant  field  and 
can  be  present  in  ferroelectric  materials  such  as  LiNb03  and  has  been  measured  in 
some  samples  to  be  as  strong  as  100  kV/cm  [72]. 

Et ,  E, v  and  Em  depend  on  K  and  therefore  on  the  interbeam  angle.  Fig  8 
depicts  the  relative  strengths  of  these  fields  as  functions  of  the  interbeam  angle  for 
Bi12Si02o  illuminated  by  light  of  wavelength  514  nm  at  a  temperature  of  293  K. 
Here  the  interbeam  angle  is  measured  in  the  medium.  As  the  angle  between  the 
beams  increases  the  grating  spacing  decreases.  This  has  the  expected  effects  of 
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increasing  Et,  but  decreasing  Eq.  For  small  interbeam  angles  (8  <  5°)  we  may  take 
£t<C  Eq.  In  all  cases  Em Eq.  Et<£.  Em only  for  very  small  (<  1°)  angles  but  for 
angles  >  10°  the  reverse  is  true. 

The  above  approximations  may  be  used  to  simplify  the  expression  for  Eal.  The 
interbeam  angle  ranges  (small/large)  where  the  approximations  are  valid  are  useful 
to  delineate  photorefractive  behaviour.  Consider  the  case  with  no  moving  fringes 
( b  =  0).  For  small  interbeam  angles  (large  grating  spacings)  we  take  Eq  E,0  'J> > 
Et.  Then 

E,i  ~  -mE, o  (51) 

Since  m  is  real,  E,i  is  in  phase  with  the  intensity  pattern  and  =  0. 

There  are  two  cases  in  which  E,i  becomes  pure  imaginary.  If  there  is  no  applied 
field,  E, o  =  0,  then  for  Eq  »  Et 

E,i  ~  jmEt  (52) 

For  large  interbeam  angles,  Eq  <  Et  and 

E. i  %  jmEq  (53) 

Here  the  modulated  space  charge  field  is  90°  out  of  phase  with  the  intensity  pattern. 

The  small  interbeam  angle  case  is  called  drift  dominated  and  is  characterized  by 
little  or  no  intensity  coupling,  but  pronounced  bending  of  the  hologram  fringes  at 
steady-state  due  to  strong  phase  coupling.  When  the  interbeam  angle  is  large  (small 
grating  spacing)  or  when  no  external  field  is  applied,  the  diffusion  process  governs 
charge  transport.  This  is  called  the  diffusion  dominated  regime.  Now  significant 
power  transfer  between  the  beams  is  possible.  Staebler  and  Amodei  [119]  were  the 
first  to  describe  these  two  extremes. 

Figs  9a-e  demonstrate  these  effects  for  Bi12Si02c  using  the  crystal  parameter 
values  of  Table  1.  In  the  first  the  phase  mismatch  angle  is  plotted  as  a  function  of 
applied  field,  E,o,  for  five  interbeam  angles  (1°,  2°,  5°,  10°,  and  20°).  At  E,0  =  0 
the  only  possible  process  for  charge  transport  is  diffusion,  thus  we  expect  =  90°. 
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Figure  8.  Characteristic  fields  versus  interbeam  angle  (in  the  medium),  (a)  Trans 
mission  geometry,  (b)  Expanded  view  of  (a). 
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Oratmg  spacing  (pm) 

0-102  0-1  0099  0-0985  0-0981 


Interbeam  angle  (degrees) 


Figure  8:  Characteristic  fields  versus  interbeam  angle  (in  the  medium),  (c)  Reflec¬ 
tion  geometry. 

For  small  interbeam  angles  even  small  electric  fields  push  towards  0  or  180°, 
while  for  large  angles,  diffusion  continues  to  dominate  over  drift.  Fig  9b  shows 
the  symmetric  dependence  of  the  magnitude  of  the  coupling  term  T  on  E,0.  It  is 
symmetric  because  we  have  no  moving  fringes  and  no  photovoltaic  field.  This  plot 
illustrates  an  important  trade-off,  the  phase  mismatch  is  more  favourable  for  power 
transfer  in  the  diffusion  dominated  case,  but  in  the  presence  of  an  applied  field  the 
coupling  coefficient,  T,  has  a  smaller  magnitude  than  for  the  drift  dominated  case. 
The  intensity  coupling  depends  on  2 T  sin  and  we  might  expect  the  magnitude 
of  this  term  to  be  optimized  for  a  particular  interbeam  angle  (or  grating  spacing). 
Fig  9c  shows  that  coupling  gets  stronger  with  increasing  interbeam  angle  for  small 
applied  fields,  but  not  necessarily  for  greater  fields.  Fig  9d  plots  2Tsin'l,9  versus 
interbeam  angle  for  E,0=  0,  2,  5,  8,  and  10  kV/cm.  There  is  a  clear  optimum  angle 
for  the  stronger  fields.  In  Fig  9e  we  plot  the  gain  (/_(d)//_( 0)),  given  by  Eqn  48,  for 
an  interaction  length  of  1  cm,  no  moving  fringes  and  an  applied  field  of  10  kV / cm 
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Interbeam  angle  (degrees) 

Figure  9:  Coupling  in  degenerate  two-wave  mixing,  (e)  Gain  (/+.  (d)/  J+(0)  )  vs 
interbeam  angle. 

for  various  incident  beam  ratios,  l3r.  The  beam  ratio  has  no  effect  on  the  optimum 
interbeam  angle.  The  smaller  /3r  is,  the  greater  the  gain.  Thus  for  flr  =  10~3  the 
gain  can  be  as  high  as  40. 

Another  key  feature  borne  out  by  Figs  9c  and  9d  is  that  the  energy  transfer  is 
just  in  one  direction.  That  is  the  coupling  term  has  the  same  sign  regardless  of  the 
combination  of  beam  ratio,  interbeam  angle,  and  applied  field.  In  our  example,  I± 
always  gains  at  the  expense  of 

For  the  applications  oriented  reader,  it  is  evident  from  Fig  9c  that  the  applied 
field  is  a  useful  control  parameter  for  beam  coupling  in  BiuSiOjo,  particularly  for 
small  interbeam  angles.  Changing  the  interbeam  angle  also  influences  the  beam 
coupling,  but  shifting  beam  directions  is  undesirable  in  applications  requiring  fixed 
optics. 
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2.3.6  Optimization  with  moving  fringes 


The  preceding  discussion  highlights  an  unfortunate  aspect  of  photorefractivitv  with 
Bii^SiOio-  An  applied  field  is  necessary  to  obtain  a  large  coupling  coefficient.  T. 
but  this  field  takes  us  away  from  the  high  gain,  diffusion-dominated  regime. 

Fortunately  there  is  a  way  to  have  the  best  of  both  worlds,  to  maximize  F  and 
have  <1>9  =  90°.  The  trick  is  to  use  moving  fringes.  This  was  first  demonstrated  ex¬ 
perimentally  in  1981  by  Huignard  and  Marrakchi  [56j.  The  following  year  Stepanov, 
et  al.  explained  this  effect  in  term.  of  “running"’  holograms  11211.  The  process  of 
reading  a  hologram  in  a  photorefractive  material  in  the  presence  of  an  electric  field 
results  in  a  moving  grating.  A  resonant  condition  is  achieved  if  the  recording  inten¬ 
sity  pattern  runs  synchronously  with  this  moving  grating.  Treatments  of  moving 
gratings  have  been  provided  theoretically  by  Valley  [129]  and  both  experimentally 
and  theoretically  by  Refregier.  et  al.  [109],  Gunter  [45],  and  Huignard  and  Mar¬ 
rakchi  [57]. 

Let  us  reconsider  Eqn  20  for  the  case  b  yt  0.  Using  the  approximations  of  Sec  2. 
we  obtain  for  large  interbeam  angles 

ft.-"?*-  <«4> 

For  negligible  6,  we  are  still  in  the  diffusion  dominated  regime  with  =  x90°.  As 
the  fringe  velocity  grows  in  magnitude  the  real  part  of  E,x  becomes  increasingly 
important.  For  \bE„/Em\  >  1,  E,x  %  -mEm/bEq.  Here  =  0  and  we  have 

no  power  transfer  between  the  beams.  Clearly  moving  fringes  will  not  facilitate 
intensity  coupling  for  large  interbeam  angles. 

For  small  interbeam  angles 


E,  i  =  —  m 


JiO 


(55) 


With  a  bit  of  algebra  we  can  show  that  the  imaginary  part  of  E,\  has  a  maximum 
for  b  =  —Em/E,o  in  the  limit  E ^  Em.  Figs  10a-d  give  a  comparative  study  of 
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Figure  10:  Coupling  in  nearly  degenerate  two-wave  mixing,  (c)  T  sin$9  vs  E,o ■  (d) 
T  sin  vs  fringe  velocity.  In  (d)  the  different  curves  refer  to  grating  spacings  (/im). 
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$3,  T  and  coupling  in  Bi12SiO20  to  that  for  zero  fringe  velocity.  We  have  chosen 
one  particular  interbeam  angle,  0.75°,  corresponding  to  a  grating  spacing  of  15  fim. 
Fig  10a  shows  as  a  function  of  E, 0  for  various  fringe  velocities.  We  see  that 
=  90°  is  realizable  for  E, 0  ^  0  with  an  appropriate  fringe  velocity.  Happily  this 
combination  of  field  and  velocity  is  also  where  T  has  its  peak,  as  seen  in  Fig  10b. 
The  resonant  condition  is  apparent  with  applied  field  as  the  control  parameter  in 
Fig  10c  and  with  fringe  velocity  as  the  control  parameter  in  Fig  lOd. 

rsin$s  for  moving  fringes  is  nearly  an  order  of  magnitude  larger  than  that 
obtained  with  no  moving  fringes.  If  we  tried  to  reproduce  Fig  9e  for  gain  with  a 
1  cm  long  crystal,  we  would  find  essentially  complete  power  transfer  between  the 
beams. 

Highly  efficient  gain  in  Bi12Si02o  and  related  crystals  may  be  achieved  using 
an  applied  a.c.  field.  First  proposed  by  Stepanov  and  Petrov  [122)  this  method 
employs  an  a.c.  field  of  frequency  u)ac  »  rd.  The  photoinduced  carriers  do  not  have 
time  to  move  far  before  the  direction  of  the  applied  field  changes.  Thus  the  time 
average  for  is  90°  and  optimum  gain  is  realized.  This  technique  has  recently 
been  used  in  photorefractive  gallium  arsenide  [78]. 

2.3.7  Other  effects 

In  the  course  of  our  treatment  of  the  photorefractive  effect  we  have  made  several 
simplifying  assumptions.  Some,  such  as  neglecting  thermal  generation  of  carriers, 
have  a  sound  basis  in  experimental  evidence.  Other  assumptions  were  made  for  the 
ease  of  the  theory.  Let  us  discuss  briefly  the  modifications  necessary  if  we  were  to 
include  these  effects. 

1.  Large  modulations.  To  linearize  the  materials  equations  we  expanded  n,  Np, 
J ,  and  E,  out  to  first  order  terms  in  a  Fourier  series  and  then  neglected  products 
of  first  order  terms.  In  other  words  we  assumed  small  modulations.  It  is  easy  to 
show  that  each  first  order  coefficient  is  proportional  to  the  intensity  modulation, 
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m.  If  the  gain  is  large  (strong  coupling  or  large  interaction  lengths)  then  m  might 
not  be  small.  This  problem  has  been  addressed  by  Refregier,  et  al.  [109],  Ochoa,  et 
al.  [97],  and  most  recently  by  Au  and  Solymar  [6] .  Higher  order  Fourier  terms  are 
included  in  a  pertubation  treatment.  Ref  [6]  modeled  Bii2SiO20  and  found  that,  in 
general,  the  imaginary  part  of  Etl  decreased  for  larger  modulation  depths.  We  also 
made  the  claim  that  Nq  ^  n,  or  specifically,  Nq1  n^.  If  we  look  at  the  a.c.  Eqns 
15-18,  then  this  latter  condition  is  no  longer  true  for  very  small  interbeam  angles 
and  very  weak  applied  fields.  For  example,  using  Bi12SiO20  with  no  applied  field,  no 
moving  fringes,  an  interbeam  angle  of  0.1°,  a  temperature  of  300°,  and  a  wavelength 
of  514  nm,  we  find  that  ~  [«j  |.  Fortunately,  most  of  our  applications  will 

involve  larger  interbeam  angles,  stronger  applied  fields,  or  both. 

2.  Absorption.  The  photorefractive  process  would  not  work  without  attenuation 
of  the  mixing  waves.  If  none  of  the  incident  light  was  absorbed,  no  photocarriers 
would  be  produced,  and  there  would  be  no  space  charge  field  or  refractive  index 
modulation.  With  moving  fringes  there  is  another  way  in  which  absorption  reduces 
the  gain.  The  total  intensity  will  not  be  constant  across  the  medium.  The  dielectric 
relaxation  time  is  inversely  proportional  to  the  total  intensity.  A  spatially  varying 
intensity  means  the  optimum  fringe  velocity  changes  across  the  crystal.  Only  one 
“slice”  of  the  crystal  will  be  at  the  resonance  condition. 

For  the  transmission  geometry  there  are  analytic  solutions  to  the  coupled  wave 
equations  in  the  presence  of  absorption  (see  Chap  4)  but  a  numerical  treatment  is 
necessary  for  the  reflection  case. 

3.  Vectorial  Effects.  We  used  a  scalar  treatment  for  the  field  equations.  But 
the  vectorial  nature  of  the  problem  is  important  when  we  consider  birefringence  or 
optical  activity.  For  example,  the  optical  activity  of  Bi12SiO20  (the  rotation  of  the 
polarisation  axis  is  38°/mm  for  514  nm)  will  reduce  photorefractive  gain.  To  see 
why,  consider  the  transmission  geometry  for  the  crystal  orientation  in  Fig  5.  Take 
the  two  incident  beams  polarised  in  the  y  direction  to  maximize  the  electrooptic 
effect.  If  the  crystal  is  10  mm  long  then  the  polarisation  of  the  waves  will  undergo  a 
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little  more  than  one  complete  revolution  during  transit.  Half  the  interaction  length 
will  be  spent  with  polarisations  close  to  the  z  direction.  But  the  electrooptic  effect 
is  null  for  beams  polarised  in  z.  There  will  be  little  or  no  coupling  in  these  regions 
and  therefore  the  overall  gain  is  reduced. 

John  Heaton’s  D.Phil.  thesis  gives  a  thorough  treatment  of  absorption,  optical 
activity,  and  birefringence  [49] .  One  must  resort  to  numerical  methods  to  include 
these  effects. 


2.4  The  Photorefractive  Effect  as  a  Nonlinear  Prob¬ 
lem 

The  term  nonlinear  has  cropped  up  many  times  in  the  course  of  this  chapter.  We 
finish  our  discussion  of  two-wave  mixing  by  reviewing  the  inherent  nonlinearity  of 
the  photorefractive  process. 

The  most  important  nonlinear  character  was  mentioned  in  Section  2.2.  The 
polarisation  of  the  medium  is  described  by  D  =  lE.  But  the  dielectric  tensor,  e, 
itself  depends  on  E,  thus  we  have  a  nonlinear  relationship.  In  particular,  e,  via  the 
electrooptic  effect  due  to  the  modulated  space  charge  field,  depends  on  the  prod¬ 
uct  of  the  mixing  optical  field  amplitudes.  When  writers  describe  photorefractive 
materials  as  being  nonlinear,  this  is  the  property  to  which  they  refer. 

The  material  equations  are  nonlinear  in  the  sense  that,  they  contain  products  of 
intensity  driven  terms  such  as  the  space  charge  electric  field  and  carrier  concentra¬ 
tion.  We  resorted  to  a  linearization  process  to  arrive  at  analytic  expressions  for  the 
modulated  space  charge  field. 

Finally  the  field  (coupled  wave)  equations  were  expressed  in  terms  of  the  prod¬ 
uct  of  complex  amplitudes.  The  lossless  case  given  by  Eqns  43  and  44  is  a  pair 
of  nonlinear  coupled  differential  equations.  Fortunately  these  equations  have  an 
analytic  solution. 


Chapter  3 

Nonreciprocal  Effects  in 
Photorefractive  Materials 


So  flew  the  god  and  the  virgin — he  on  wings  of  love,  and  she  on  wings  of  fear. 

The  Legend  of  Apollo  and  Daphne 
Bullfinch’s  Mythology 


In  Chapter  2  we  learned  that  the  phases  of  two  mixing  waves  could  be  coupled 
in  a  nonlinear  photorefractive  medium.  The  particular  case  of  counterpropagating 
beams  was  investigated  first  by  Yeh  [136,139].  A  1985  Letter  by  Solymar  and  Wilson 
proposed  a  novel  application  for  this  phase  coupling  [118].  They  suggested  using 
photorefractive  materials  in  a  controllable  nonreciprocal  (NR)  optical  phase  shifter. 
This  is  a  device  through  which  a  beam  transitting  in  one  direction  will  suffer  a 
certain  phase  shift,  but  a  beam  travelling  in  the  opposite  direction  may  undergo  a 
different  phase  shift. 

NR  devices,  such  as  phase  shifters,  are  well  known  in  microwave  applications  [23] 
and  more  recently  have  become  important  in  optical  waveguides  [98,17,92].  Many 
applications  involve  a  180°phase  shift  (gyrator)  and  the  devices  invariably  depend  on 
a  gyrotropic  or  Faraday  rotation  effect.  NR  phase  shifts  of  up  to  360°  are  available 
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using  these  devices  so,  to  be  competitive,  a  photorefracti ve  NR  phase  shifter  should 
be  able  to  produce  similar  sized  effects.  Solymar  and  Wilson  suggested  using  the 
NR  phase  shifter  in  a  fibre-optic  gyroscope  to  compensate  for  the  Sagnac  phase 
shift.  This  operation  would  be  analogous  to  slewing  a  mechanical  gyroscope  to  a 
new  orientation.  Yeh  has  proposed  using  the  effect  to  decouple  counterpropagating 
modes  in  a  bidirectional  ring  resonator  [139].  This  latter  use  has  recently  been 
demonstrated  in  a  resonator  using  BaTi03  as  the  gain  medium  [62]. 

In  this  chapter  we  will  explore  the  feasibility  of  a  NR  phase  shifter  using  pho- 
torefractive  materials.  First  we  look  at  an  analytic  treatment  of  two-wave  mixing 
for  counterpropagating  waves  [59,136].  We  then  argue  that  this  approach  requires  a 
numerical  treatment  for  several  important  cases,  and  we  will  compare  results  using 
analytical  and  numerical  methods.  Once  we  have  an  effective  means  to  solve  the 
coupled  wave  equations,  we  can  look  in  detail  at  NR  transmission  and  phase  shift, 
the  expected  magnitude  of  these  effects,  and  the  means  to  control  them.  The  chap¬ 
ter  finishes  with  an  overall  appraisal  of  using  photorefractives  for  NR  applications 
and  a  discussion  of  the  limitations  of  our  treatment. 

3.1  Reflection  Holograms 

3.1.1  Background 

Most  applications  of  two- wave  mixing  in  photorefractive  materials  use  transmission 
geometry.  There  are  several  reasons  why  reflection  geometry  is  less  popular. 

1.  Large  interbeam  angles,  such  as  would  be  used  in  reflection  geometry,  are 
associated  with  the  diffusion  dominated  region  with  the  phase  mismatch  angle, 
$ff,  close  to  90°.  Unfortunately,  overall  gain  is  generally  reduced.  For  example, 
in  Chapter  2  we  found  in  Bi12Si02o  the  interbeam  angle  for  optimum  gain  was 
usually  less  than  10°. 
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2.  Reflection  geometry  often  involves  different  crystal  orientations  than  those 
for  the  transmission  case  to  exploit  fully  the  electrooptic  effect.  This  requires 
different  polished  faces,  or  perhaps,  an  alternative  crystal  cut.  There  is  a 
greater  demand  for  transmission  geometry  orientations,  and  these  are  what 
the  crystal  suppliers  provide. 

3.  Experimental  layout  for  reflection  geometry  is  generally  more  complicated.  In 
general,  a  greater  number  of  optical  elements  are  required  to  get  a  split  beam 
around  the  crystal  to  send  it  through  the  rear  face. 

4  In  the  reflection  case  the  coupled  wave  equations  have  a  two-point  boundary 
condition.  This  increases  the  difficulty  of  a  mathematical  solution. 

There  are  certainly  fewer  papers  dealing  with  reflection  geometry  as  opposed 
to  transmission  geometry  for  photorefractive  materials.  Counterpropagating  beams 
impinging  on  a  crystal  exhibiting  optical  activity  was  studied  by  Kukhtarev,  Dov- 
galenko  and  Starkov  [74],  Diffraction  efficiency  in  both  reflection  and  transmission 
geometries  was  the  subject  of  a  paper  by  Heaton,  et  al.  [47]  A  reflection  grat- 
ing  can  form  when  a  single  beam  interacts  with  its  reflection  off  the  back  face 
of  a  crystal.  The  treatment  of  this  case  was  studied  theoretically  by  Yeh  [137] 
and  experimentally  by  MacDonald,  et  al.,  in  BaTi03[86],  and  by  Zha  and  Gunter 
in  KNbO3[140].  Ja  generalized  two-wave  mixing  with  multiple  reflections  [58]. 
Dynamic  self-enhancement  in  lithium  niobate  using  reflection  geometry  has  been 
demonstrated  by  Shvarts,  et  al.  [113] 

3.1.2  An  analytic  treatment 

Two  of  the  most  detailed  theoretical  studies  of  mixing  of  counterpropagating  beams 
come  from  Yeh  [136]  and  Ja  [59].  This  latter  paper  includes  a  numerical  treatment 
for  the  lossy  case.  The  method  for  the  lossless  case  follows. 

We  start  with  the  coupled  wave  equations: 
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dl+ 

dz 

dl_ 

dz 

d<t>+ 

dz 
d<f> _ 
dz 


To  solve  the  intensity  equations,  introduce  A  =  1+  —  /_  and  £  =  1+  +  Then 
we  can  rewrite  Eqns  1  and  2  as 


2T  sin  $9 

+  /_ 

(1) 

2r  sin  y~~j~ 

1+  T  1- 

(2) 

r  cos  $9 

i+  -r  /_ 

(3) 

T  cos  $g 

9 1+  +  I- 

(4) 

dA 

dz 

dZ 

dz 

These  two  equations  lead  us  to 


0  ,  A  =  constant 

r  •  *  SJ  ~ 

1  sin  <P„  — — - 


(5) 

(6) 


I+(z)  =  i/+0  (7l  -  M2  +  40roe2r,in*,z  +  (1  -  &„))  (7) 

*-(*)  =  \I+a  (v/l  -  &„)2  +  40roe’™***  -  (1  -  &.))  (8) 

where  /+o  =  /+(0)  and  /?ro  =  /_(0)//+(0).  The  above  two  expressions  are  in¬ 
convenient  in  that  we  rarely  know  /_(0)  in  advance.  It  is  useful  to  introduce 

0r  =  /_(d)//+(0),  the  input  beam  ratio,  where  d  is  the  crystal  length.  It  is  a 
straightforward  calculation  to  show 


0ro 


0r(l  +  dr) 

(3r  +  e2r,itt*'d 


(9) 


/+(*)  =  -C  +  \JC 1  +  Bejr,in*'*  (10) 

!_(*)  =  C  +  \jc2  +  JSe,r“n*»*  (11) 


Using  this  we  obtain 
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B  =  3roI+02 

_  a  t  2  __L±A _ 

“  ft  +  e2r,i“**d 

C  =  i/+o(/3r„-l) 

ft2  -  e2r,in*«d 
-  2/+0  f3r  +  e2rin i'd 

Eqns  10  and  11  can  be  put  into  the  phase  equations  to  obtain 

<£+(2)  =  <j>+(0)  +  jTcos^gZ  +  *  cot  $sln(TVr_) 
<P-{z)  =  -  ^rcos$ff2  +  j  cot  ^sin(r+r_) 


(12) 


(13) 


(14) 

(15) 


where  T+  =  /+(d)//+(Q)  and  T_  =  /_(0)//_(d)  are  the  transmissions  of  the  two 
waves.  A  useful  relationship  is 


/+(*)/_(*)  =  Se2r,m*»z 


(16) 


We  can  use  this  to  show  the  transmissions  are 

_  (l+ft)e2r»in*»d 
+  ft  +  e2r  "n**d 

r  =  i  +  ft 

ft  +  e2r.in*fd 

We  next  obtain  the  expression  for  the  NR  phase  shift,  Aft 


(17) 

(18) 


A  <(> 


(<t>+(d)  -  M 0))  -  (0_(o)  -  4>.{d)) 


In  the  limit  — *  0,  Eqn  20  becomes 


(19) 

(20) 


a  <t>  =  |r«f 


ft  -  1 
ft  +  1 


(21) 


The  above  expression  is  what  prompted  Solymar  and  Wilson  to  suggest  a  NR 
phase  shifter.  A  large  product  I'd  for  an  input  beam  ratio  much  different  than 
1  can  lead  to  a  substantial  NR  phase  shift.  Numerical  results  from  the  previous 
chapter  suggest  values  for  T  close  to  15  cm-1  for  BijjSiOao  are  possible.  This  was 
for  transmission  geometry  and  moving  fringes.  The  question  we  need  to  ask  is  how 
large  can  the  coupling  term  be  for  reflection  geometry? 
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Figure  1:  Power  conservation  in  two- wave  mixing,  (a)  Transmission  geometry,  (b) 
Reflection  geometry. 

3.1.3  A  numerical  treatment 

To  answer  the  question  of  the  size  of  the  coupling  term  we  need  to  obtain  an 
expression  for  the  modulated  space  charge  field,  E,\ .  Before  using  the  procedure 
of  Chapter  2,  we  should  consider  power  conservation  in  the  reflection  geometry. 
Fig  1  gives  a  comparison  of  the  two  geometries.  In  the  transmission  case  the  total 
intensity  is  constant.  For  wave  mixing  in  the  reflection  case,  A  =  1+  —  In  our 
analytic  treatment  this  constant  difference,  A  =  —2 C.  The  important  feature  is 
that  the  total  intensity  varies  across  the  crystal. 

How  will  a  background  intensity  gradient  affect  our  previous  treatment?  The 
d.c.  equations  in  Section  2.1.6  will  not  be  changed.  This  means  the  zero  order  (or 
background)  current  density,  J0,  is  constant.  This  is  a  flow,  and  is  the  product 
of  a  density  and  a  velocity.  The  density  is  given  by  the  carrier  density,  n0,  which 
is  proportional  to  the  total  intensity  in  a  given  region.  The  velocity  is  the  drift 
velocity,  and  this  is  proportional  to  the  field  E,q.  Since  J0  is  constant  but  n0  is  not, 
this  means  E&  must  vary  in  space.  This  is  an  effect  not  considered  by  previous 
researchers.  There  is  an  additional  consideration  when  moving  intensity  fringes  are 
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z=0  z  =  d 


Figure  2:  Nearly  degenerate  two-wave  mixing  of  counterpropagating  waves. 

present.  This  is  essentially  the  same  argument  given  in  Chapter  2  for  gain  reduction 
when  absorption  is  present.  Coupling  is  optimized,  for  a  given  total  intensity,  for 
one  particular  fringe  velocity.  If  the  total  intensity  varies  across  the  crystal,  then 
one  fringe  velocity  cannot  satisfy  the  optimum  condition  everywhere — the  gain  is 
less. 

The  treatment  of  Sec  2.1  can  be  modified  easily  to  include  a  spatially  varying 
intensity.  Here  the  input  waves  are  as  shown  in  Fig  2.  Analogous  to  Eqns  2.5  and 
2.10  the  driving  intensity  and  normalized  quantities  have  the  forms: 

I(z,t)  =  /«(*)+ /,(*)e*'jr^-)+c.e.  (22) 

Fn(z,t)  =  Fn0(z)  A  Fnl(z)ejK^  +  c.c.  (23) 

We  assume  steady-state  conditions  for  the  coefficients.  The  procedure  to  find  Efl  is 
exactly  as  outlined  in  Chapter  2,  except  all  the  terms  involving  the  total  intensity, 
such  as  b  and  h  can  have  z  dependence.  One  final  assumption  is  needed,  namely, 
the  zero  order  terms  vary  slowly  in  space.  On  the  scale  of  a  grating  spacing  their 
derivatives  with  respect  to  z  are  negligible  compared  with  those  of  the  first  order 


terms. 
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With  these  assumptions,  we  obtain  the  following  results 


E,o{z) 


E,\(z) 


fien0 


m _ E,0(z)  4-  Ev  -  jE, _ 

m  1  +  +  6(z)^  -  H*)(l  +  fM) 


(24) 

(25) 


We  should  mention  that  in  the  derivation  of  the  above,  Jo  was  found  to  be  constant. 
As  long  as  a  current  is  present,  E, 0  varies  spatially,  since  the  background  carrier 
density,  n0,  depends  on  the  total  intensity.  In  this  case  a  further  condition  on  E, 0 
is 


V  =  [d  E.0(z)dz  (26) 

Jo 

where  V  is  any  applied  voltage  (possibly  equal  1 3  zero  for  the  case  where  only  the 
photovoltaic  field  drives  the  current).  If  EM  varies  spatially,  so  does  2?fl.  If  there  is 
no  current,  then  E.0  is  constant,  but  E,\  will  still  be  spatially  dependent  if  b  ^  0. 

Previous  experiments  demonstrating  NR  transmission  used  crystals  in  open  cir¬ 
cuit  configurations  [86,140],  that  is  no  electrodes  or  wires  and  no  current  flow.  If 
no  voltage  is  applied  and  J  =  0,  then  E&  is  constant.  This  leads  to 


£.iopen(*)  =  m- 


jEt 


(27) 


This  is  identical  to  Eqn  25  for  E,0  —  ~EP,  except  the  last  term  in  the  denominator 
was  6(1  +  Et/Em).  For  reflection  geometry,  Et/Em  3>  1,  so  the  two  expressions 
reduce  to  the  same  thing.  In  other  words,  the  open  circuit  configuration  can  be 
included  as  a  special  case  of  the  closed  circuit  treatment,  for  E,o  —  -Ep. 

The  coupled  wave  equations  are  given  by  Eqns  1-4,  except  now  T  and  may 
vary  spatially. 


3.1.4  A  comparison  of  the  treatments 

It  is  important  to  compare  results  using  variable  coupling  with  those  using  constant 
coupling.  The  two  crystals  we  will  consider  are  BiijSiOjo  and  LiNb03.  LiNb03 
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Figure  3:  Crystal  orientations  for  counterpropagating  waves.  (a)  LiNb03. 
(b)  Bii2Si02o- 

has  the  advantages  of  a  larger  electrooptic  coefficient  and  a  photovoltaic  effect.  The 
crystal  orientations  for  counterpropagating  beams  are  shown  in  Fig  3.  In  both  cases 
any  applied  voltage  will  be  in  the  z  direction. 

Consideration  of  the  electrooptic  effect 

The  photovoltaic  field  in  LiNb03  acts  along  the  001  axis.  To  take  full  advantage  of 
this  field,  the  two-wave  mixing  grating  vector  needs  to  be  in  this  direction  as  well. 

LiNb03  belongs  to  crystal  class  3m.  In  the  contracted  notation  of  Yariv  and 
Yeh  (135j,  the  electrooptic  tensor  is 
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An  electric  field  applied  along  the  001  (or  z  or  crystal  c)  axis  does  not  rotate  the 
principal  axes,  but  changes  their  lengths. 

AnIlV  =  -\n3orlzE  (28) 

An*  =  -\nlr33E  (29) 

For  light  polarised  in  either  the  x  or  y  direction  the  relevant  electrooptic  coefficient 
is  r13.  This  component  is  roughly  one  third  that  of  the  largest  coefficient  (r33),  but 
this  orientation  was  chosen  specifically  to  take  full  advantage  of  the  photovoltaic 
field. 

For  Bi12SiO20  the  field  will  be  along  the  001  axis  to  optimize  the  electrooptic 
effect.  This  case  was  worked  out  in  Chapter  2  and  we  will  use  the  same  input 
polarisations,  either  in  the  x  or  y  directions. 

The  numerical  method 

The  variable  coupling  terms  and  the  mixed  boundary  conditions  require  a  numerical 
treatment  for  the  solution  of  the  coupled  wave  equations.  A  common  technique  to 
handle  mixed  boundary  conditions  is  the  shooting  method  [104,61].  The  basic 
approach  is  to  make  reasonable  guesses  for  all  unspecified  function  values  at  one 
boundary,  integrate  to  the  other  boundary  and  compare  the  computed  function 
values  with  the  known  boundary  conditions.  Revised  guesses  are  made  and  the 
cycle  repeats  until  the  final  function  valuer  arc  within  a  cct  tolerance. 

In  our  case  we  do  not  know  /_(0)  in  advance.  We  use  the  analytic  solution  in 
Eqn  11  to  make  our  first  guess.  (Here  we  take  the  total  intensity  as  /+(0)  +  I-{d) 
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and  E, o  =  V!d  in  calculating  *J>,  and  r.j  /_ (0 )  is  known,  and.  for  convenience,  we 
take  0  4.(0)  =  o_(0)  =  0.  This  is  justified  since,  as  revealed  in  Eqns  14  and  15.  there 
is  an  undetermined  constant  in  each  solution. 

The  set  of  four  coupled  differential  equations  is  integrated  using  a  fourth-order 
Runge-Kutta  routine.  If  the  calculated  value  for  /-(d)  is  not  within  a  tolerance  of 
0.1  %  then  a  new  guess  for  /_(0)  is  found  via  a  Newton-Raphson  method. 

A  further  integration  is  required  after  each  iterative  cycle  to  choose  J0  such  that 
the  condition  in  Eqn  26  is  satisfied.  Combining  Eqns  25  and  26  gives 

.  V  +  E„d 

J  o  —  11C — j - 

Jo  a’-z /n0(/) 

where  the  total  intensity,  /(z),  is  updated  on  each  iteration. 

Figs  4-6  show  comparisons  between  the  numerical  method  outlined  here  and  the 
analytic  and  numerical  methods  of  Refs  [59,136]  which  assume  constant  coupling. 
In  all  plots  the  crystal  length  is  taken  to  be  1  cm,  the  wavelength  of  the  incident 
light  is  514  nm,  and  the  temperature  is  300  K.  For  LiNb03  the  photovoltaic  field  is 
chosen  to  be  15  kV/cm,  a  typical  value. 

The  coupling  in  LiNb03  is  strong  and  there  is  significant  variation  of  intensity 
across  the  crystal.  This  is  seen  in  both  the  intensity  and  phase  solutions  to  the 
coupled  wave  equations  (Figs  4a, b  and  5a, b).  Figs  4c  and  5c  depict  the  variation  of 
Erf.  Figs  4d  and  5d  illustrate  that  the  coupling  terms  do  vary  across  the  crystal. 

The  conditions  for  Fig  4  were  chosen  to  illustrate  the  case  of  no  moving  fringes. 
The  input  beam  ratio  is  0.1  and  the  variable  coupling  solution  for  intensity  shows 
less  power  transferred  between  the  beams.  The  difference  between  the  analytic  and 
numerical  treatments  shown  in  Fig  5a,  where  only  the  photovoltaic  field  drives  the 
current,  is  not  as  disparate.  That  is  because  the  intensity  coupling  term,  F  sin^9, 
does  not  have  a  large  spatial  variation.  What  is  interesting,  though,  is  that  the 
phase  coupling  term,  rcos$B,  changes  sign  as  we  move  across  the  crystal.  The 
difference  in  behaviour  for  the  phase  solutions  between  the  two  methods,  as  seen  in 
Fig  5b,  is  quite  dramatic. 
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(a)  (b) 


(c)  (d) 


Figure  5:  LiNb03.  V  =  0,  /+(0)  =  100  mW/cm2,  J0=  5.52  x  lO’6  A/m2,  fringe 

velocity  -  -5.0  x  10~9  m/s,  no  absorption.  Intensity  is  normalised  to  /+(0),  phase 
is  normalised  to  n  v 
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Fipire  6:  Bi12Si02o  V  -  10  kV,  /+(0)  =  10  mW/cm2,  J0=  0.0456  A/m2,  fringe 

velocity  -  1.5x  10~  m/s,  a  =  2  cm'1.  Intensity  is  normalised  to  U{ 0),  phase  is 
normalised  to  n . 
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Coupling  is  not  strong  in  BijiSiO^o  and  the  two  treatments  yield  essentially  the 
same  solutions.  If  absorption  is  present,  though,  the  intensity  variation  across  the 
crystal  is  greater,  and  some  divergence  between  the  methods  might  be  expected. 
Since  we  are  using  a  numerical  method  to  treat  the  variable  coupling  case,  it  is 
relatively  simple  to  include  absorption  in  our  treatment.  The  intensity  coupled 
equations  become 


dU 

dz 

IJ_ 

_  2r  sin  <1>S  a/j. 

1+  -t- 

(31) 

dl. 

dz 

=  2r  sin  4>3  — - - a/_ 

1+  -r  I- 

(32) 

We  can  compare  this  with  the  numerical  treatment  with  absorption  for  constant 
coupling  in  Ref  [59] .  Figs  6a-d  consider  absorption  in  the  case  of  Bi^SiOjg.  Here  ot 
is  taken  as  2  cm'1,  a  value  found  experimentally  for  some  Bi^SiOzo  samples  [31]. 
The  difference  in  the  two  treatments  is  more  apparent  in  the  phase  solutions. 

We  argued  at  the  beginning  of  this  subsection,  from  physical  principles,  that 
moving  fringes,  absorption,  or  a  photocurrent  leads  to  variable  coupling  in  reflection 
geometry.  That  was  confirmed  in  the  derivation  of  E, j.  That  this  variable  coupling 
is  a  significant  effect  is  seen  in  Figs  4-6. 

We  have  a  suitable  numerical  treatment  with  which  to  solve  the  coupled  wave 
equations  for  counterpropagating  beams.  Now  we  have  the  means  to  calculate  NR 
effects  and  consider  various  control  parameters.  Two  potential  control  parameters 
are  applied  voltage  and  fringe  velocity.  It  is  unfortunate  that  these  represent  the 
cases  where  we  most  need  the  numerical  treatment. 


3.2  Nonreciprocal  Effects 

The  origin  of  the  laurel  tree,  according  to  ancient  Greeks,  began  with  an  unkind 
remark.  One  day  Apollo  chided  Cupid  for  playing  with  arrows — such  devices  were 
surely  meant  for  more  mature  hands.  Apollo’s  assessment  of  Cupid’s  character  was 
spot  on,  for  no  sooner  had  he  turned  his  back,  but  Cupid  fired  a  gold-tipped  arrow 
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Figure  7:  Principle  of  reciprocity.  Source  at  P\  and  detector  at  may  be  inter¬ 
changed  and  still  produce  the  same  effect. 

at  him.  This  caused  Apollo  to  fall  hopelessly  in  love  with  Daphne,  a  wood  nymph. 
The  god  of  love  then  struck  Daphne  with  a  lead-tipped  arrow.  Her  heart  hardened 
and  she  could  love  no  man,  particularly  Apollo.  She  fled  and  Apollo  pursued.  He 
had  nearly  overtaken  the  nymph  when  her  father,  Peneus,  interceded  and  changed 
her  into  a  laurel  tree.  This  is  an  interesting  and  early  example  of  the  dramatic 
consequences  of  nonreciprocity. 

To  understand  NR  effects,  we  should  start  with  the  principle  of  reciprocity.  Fig 
7  illustrates  this  principle,  most  commonly  attributed  to  Helmholtz  [42].  A  source 
and  a  detector  may  be  placed,  one  at  P\  and  the  other  at  P2.  If  the  system  through 
which  the  light  signal  passes  is  linear,  then  reciprocity  holds:  the  measured  signal 
(amplitude  and  phase)  will  be  the  same  whether  the  source  is  at  P\  and  the  detector 
is  at  Pi  or  vice  versa.  If  the  intervening  system  is  nonlinear,  then  NR  behaviour  is 
possible. 

For  two-wave  mixing  in  photorefractive  materials  there  are  two  types  of  NR 
effects — transmission  and  phase  shift. 
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3.2.1  Nonreciprocal  transmission 

NR  transmission  is  merely  a  fancy  name  for  power  transfer.  If  one  mixing  wave 
gains  at  the  expense  of  the  other,  then  T_  =  Tl.  This  effect  was  demonstrated 
experimentally  in  KNb03:Mn  by  Zha  and  Gunter  1401.  For  an  open  circuit  crystal 
with  beams  of  equal  input  intensity,  they  measured  T 1  =  113%  and  XL  =  57%. 
MacDonald,  et  al..  used  BaTi03  to  investigate  "assymmetric  transmission’’  [85).  A 
single  input  beam  was  used  and  the  counterpropagating  beam  was  the  back  surface 
reflection  of  the  input.  They  found  the  transmission  through  the  crystal  depended 
on  whether  the  input  was  parallel  to  the  crystal  c  axis  (T  —  52%)  or  antiparallel 
(T  =  39%).  The  asymmetry  is  a  consequence  of  the  one-way  power  transfer  in 
photorefractive  crystals.  For  one  orientation  of  the  c  axis  the  reflected  beam  gained 
at  the  expense  of  the  input,  for  the  opposite  orientation  it  lost  power.  A  complete 
treatment  of  the  problem  involves  multiple  reflections.  The  reflection  at  the  rear 
surface  depends  on  the  steady-state  intensity  solution.  This  is  worked  out  in  detail 
in  Ref  [137].  If  this  effect  could  be  enhanced,  one  could  use  a  photorefractive  crystal 
as  an  optical  diode. 

The  potential  control  parameters  for  NR  effects  are  applied  voltage,  beam  ratio 
and  fringe  velocity  We  investigate  these  controls  using  LiNb03  as  the  medium, 
and,  in  most  cases,  consider  the  effect  of  different  values  for  the  photovoltaic  field. 
In  ail  plots  the  crystal  length  was  1  cm,  the  temperature  300  K,  and  the  wavelength 
514  nm.  The  crystal  parameters  are  taken  from  Table  2.1. 

Figs  8  a  and  b  show  the  transmissions  as  functions  of  applied  voltage  for  input 
beam  ratios,  /?,,  of  1  and  0.1  with  no  moving  fringes.  Substantial  power  is  transferred 
from  the  plus  beam  to  the  minus  beam.  Physically  this  is  due  to  being  in  the 
diffusion  dominated  regime  and  having  a  long  interaction  length.  The  photovoltaic 
field  is  taken  to  act  in  the  +z  direction,  so  a  positive  polarity  voltage  enhances  the 
NR  transmission,  while  a  negative  polarity  detracts. 

Transmission  (T_)  as  a  function  of  beam  ratio  is  shown  in  Fig  9  for  no 


CHAPTER  3:  SONRECIFROCAL  EFFECTS 


Figure  8:  Transmission  vs  applied  voltage,  (a)  (3r=  1.  (b)  0T—  0.1 
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Figure  9:  Transmission  vs  beam  ratio. 

applied  voltage  and  no  moving  fringes.  A  photovoltaic  field  increases  the  trans¬ 
mission,  and  the  gain  is  larger  for  smaller  beam  ratios. 

Some  interesting  effects  are  seen  when  moving  intensity  fringes  are  present  (Fig 
10a).  Here  no  voltage  is  applied,  the  beam  ratio  is  1,  and  the  total  input  intensity 
is  200  mW/cm2.  With  no  photovoltaic  field  the  coupling  is  still  strong  and  the 
behaviour  is  symmetric.  As  this  field  increases  the  asymmetry  with  fringe  velocity 
is  more  apparent.  The  maximum  NR  transmission  is  no  longer  at  v  =  0,  and 
for  positive  fringe  velocity  it  is  possible  to  change  the  direction  of  power  transfer. 
The  key  feature  here  is  that  moving  fringes  can  be  used  to  reduce  the  NR  effect. 
The  magnitude  of  fringe  velocity  (~  10-8  m/s)  at  this  relatively  strong  intensity  is 
very  small.  This  corresponds  to  a  detuning  frequency  of  0.7  Hz.  If  the  intensity  is 
<  reduced,  the  detuning  frequency  required  to  achieve  the  same  effect  goes  down. 

In  Fig  10b,  Ev  =  15  kV/cm,  V  =  10  kV,  0r=  0.1,  and  the  total  intensity  is  200 
mW/cm2.  Here  we  plot  J0  as  well.  Notice  that  J0  has  a  minimum  close  to  the  fringe 
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velocity  at  which  the  NR  transmission  is  optimized.  Physically  we  can  understand 
this  by  looking  at  Fig  4a  for  the  solution  to  the  intensity  coupled  wave  equations. 
The  steady-state  background  current  density  depends  on  the  background  carrier 
density,  which  is  driven  by  the  total  density.  For  very  strong  coupling,  the  space 
average  total  intensity  is  less  than  the  total  input  intensity.  This  means  the  current 
density  is  reduced.  This  is  a  useful  effect  for  device  applications.  We  can  monitor 
beam  coupling  by  measuring  the  photocurrent — there  is  no  need  to  sample  either 
of  the  beams. 

3.2.2  Nonreciprocal  phase  shift 

Phase  coupling  between  counterpropagating  beams  was  evident  in  Figs  4-6.  Thus 
we  expect  a  difference  in  phase  shifts  encountered  by  the  two  mixing  waves.  We 
need  to  investigate  the  influence  of  various  control  parameters  on  this  effect. 

We  define  NR  phase  shift,  &<j>,  as  in  Eqn  19.  We  use  LiNb03,  a  crystal  length 
of  1  cm,  temperature  of  300  K,  and  a  wavelength  of  514  nm.  Figs  11a  and  b  show 
A <f>  as  a  function  of  applied  voltage  for  no  moving  fringes  and  beam  ratios  of  1  and 
0.1.  Here  the  NR  phase  shift  is  normalised  to  ir.  A  photovoltaic  field  is  necessary  to 
produce  an  appreciable  effect.  It  is  evident  the  sign  of  A</>  depends  on  beam  ratio. 
This  is  shown  in  Fig  12  for  an  applied  voltage  of  10  kV. 

The  effect  of  fringe  velocity  is  seen  in  Figs  13a-d.  Here  V  =  10  kV  and  /+(0)  = 
100  mW/cmz.  Three  values  of  (3r,  1,  0.4,  and  0.01,  are  shown.  Note  the  difference 
in  behaviour  for  the  plots  0T  =  1  and  f3r  =  0.01.  The  sign  of  A<£  for  a  given  fringe 
velocity  is  essentially  reversed  between  the  two  cases.  3r  =  0.4  lies  in  the  transition 
range. 

The  largest  NR  phase  shift  is  roughly  — tt/2  for  Ep  =  15  kV/cm  and  0T  =  0.01. 
Fig  13d  shows  the  close  relationship  between  input  intensity  and  optimum  fringe 
velocity.  The  beam  ratio  in  this  figure  is  kept  constant  while  the  intensity,  as  given 
by  /+(0),  varies.  Each  curve  represents  a  different  fringe  velocity  and  each  one  has 
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Figure  12:  NR  phase  shift  vs  beam  ratio. 

an  optimum  for  a  different  value  of  i+(0),  (The  guiding  principle  for  photorefractives 
is  the  greater  the  intensity,  the  faster  the  reaction.)  The  optimum  value  of  A<p  is 
the  same  for  each  curve.  This  suggests  that  a  device  capable  of  seeking  out  this 
optimum  could  be  used  to  measure  intensity  by  varying  fringe  velocity,  or  to  detect 
movement  by  varying  intensity. 


3.3  Assessment  and  Design  Considerations 

3.3.1  Control  parameters 

In  the  previous  section  three  control  parameters  were  studied — applied  voltage, 
beam  ratio,  and  fringe  velocity.  (The  effect  of  a  fourth  parameter,  the  input  intensity 
was  found  to  be  dependent  upon  the  fringe  velocity.)  Let  us  discuss  each  control 
separately. 

1.  Applied  voltage.  In  the  absence  of  a  photovoltaic  field,  applied  voltage  is 
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Figure  13:  NR  phase  shift  vs  fringe  velocity,  (a)  0r=  1.  (b)  (3, 
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not  a  useful  control  parameter  for  NR  behaviour.  Its  effectiveness  is  increased  for 
smaller  beam  ratios.  To  apply  a  voltage,  transparent  electrodes  are  required. 

2.  Beam  ratio.  In  most  practical  applications,  the  input  beam  ratio  will  be  fixed 
in  advance,  or  allowed  to  reach  some  steady-state  value  {such  as  in  a  feedback  circuit 
like  a  bidirectional  ring  resonator).  Thus  this  parameter  has  limited  potential  as  a 
control.  In  general,  the  lower  the  input  beam  ratio,  the  greater  the  NR  effects.  For 
a  given  fringe  velocity,  the  beam  ratio  determines  the  sign  of  the  NR  phase  shift. 

3.  Fringe  velocity.  This  is  the  best  candidate  for  a  control  parameter  as  it 
offers  the  widest  range  for  NR  effects.  The  reason  it  is  so  effective  is  that  coupling 
between  the  beams  is  only  strong  for  a  small  range  of  fringe  velocities  close  to  0.  If 
a  strong  photovoltaic  field  exists,  one  can  use  fringe  velocity  to  change  the  direction 
of  power  transfer.  There  are  two  drawbacks  to  using  fringe  velocity  as  a  control. 
In  many  applications  it  will  be  undesirable  to  frequency  shift  one  beam  relative  to 
the  other.  In  such  a  case  moving  fringes  can  be  produced  by  moving  the  crystal  at 
the  desired  speed  by  mounting  the  crystal  on  a  piezoelectric  translator  driven  by  an 
applied  sawtooth  voltage.  The  second  drawback  is  the  very  slow  velocities  required 
to  achieve  strong  NR  effects,  even  for  moderate  intensities.  Fringe  velocity  may  be 
impractical  to  use  for  low  intensity  two-wave  mixing. 


3.3.2  Comparison  of  crystal  types 

The  previous  section  concentrated  on  LiNbOa  as  the  nonlinear  medium.  Table  1 
compares  NR  effects  for  three  crystal  types  (Bi12SiOjo:  BaTiC>3,  and  LiNb03),  using 
fringe  velocity  as  the  control.  These  are  representative  values  for  similar  conditions. 
In  each  case  the  crystal  length  was  1  cm,  the  wavelength  514  nm,  the  temperature 
300  K,  the  input  beam  ratio  was  1,  the  total  input  intensity  was  200  mW/cm2,  and 
the  beams  propagated  along  the  crystal  c  axis.  There  will  be  some  minor  differences 
between  the  crystals,  though,  to  increase  NR  effects — namely,  the  LiNbC>3  sample 
should  have  a  photovoltaic  field  (we  took  the  value  10  kV/cm  in  our  comparison), 
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Maximum 

Fringe  Velocity 

Crystal 

Transmission 

NR  Phase  Shift 

Scale  (m/s) 

* '  - — 

Bil2Si02o 

1.33  i 

6.1  x  KT3- 

10~6 

BaTi03 

wmmmm 

1CT7 

LiNb03 

HBBH 

jj 

1(T8 

Table  1:  Comparison  of  NR  effects  for  different  crystals. 

while  in  Bi^SiOjo  we  should  apply  a  voltage  (here  we  used  a  value  of  10  kV  ). 
LiNb03  seems  to  give  the  best  combination  of  transmission  and  NR  phase  shift, 
but  requires  the  slowest  fringe  movement.  If  a  strong  photovoltaic  field  is  present, 
one  can  use  fringe  velocity  to  change  the  direction  of  power  transfer. 

3.3.3  Magnitude  of  nonreciprocal  effects 

This  study  revealed  that  the  largest  calculated  NR  phase  shift  (<  7r/2)  is  less  than 
expected  from  Eqn  21.  In  that  expression,  when  — *  0  then,  for  small  beam 
ratios,  A<f>  «  Vd.  The  difficulty  with  this  simple  analysis  is  that  V  and  do 
not  have  constant  values  for  the  combinations  of  control  parameters  necessary  for 
drift  domination.  But  we  can  consider  what  happens  when  intensity  coupling  is 
reduced — this  should  correspond  to  conditions  where  phase  coupling  is  strong.  For 
example,  in  Fig  10a  there  are  specific  values  of  fringe  velocity  where  T+  =  T_.  For 
these  velocities,  the  spatial  average  of  should  be  close  to  0.  When  we  compare 
&4>  in  Fig  13a  we  see  these  velocities  are  not  where  the  NR  phase  shift  is  optimized. 
In  other  words,  conditions  which  take  the  average  value  of  to  0  also  tend  to 
reduce  the  magnitude  of  T  This  makes  sense,  T  and  $9  are  not  independent,  but 
are  linked  through  E,\. 


CHAPTER  3:  SOS  RECIPROCAL  EFFECTS 


Figure  14:  Measurement  of  NR  phase  shift  for  nondegenerate  two-wave  mixing. 

3.3.4  Measurement  of  nonreciprocal  phase  shift 

It  is  difficult  to  measure  the  phase  of  a  wave  directly.  Instead  it  is  common  practice 
to  combine  the  wave  with  a  reference  wave  (of  unspecified  phase)  and  measure  the 
phase  difference.  This,  in  essence,  is  what  is  done  in  an  interferometer.  In  a  similar 
way  we  do  not  measure  the  NR  phase  shift,  but  its  change  as  we  adjust  control 
parameters.  A  scheme  involving  frequency  shifted  input  beams  is  illustrated  in  Fig 
14.  After  mixing  in  the  crystal  each  transmitted  beam  is  combined  with  a  reference 
beam  in  a  Mach-Zehnder  interferometer  arrangement  to  produce  a  moving  inter¬ 
ference  fringe  pattern  which  can  be  measured  bv  photodetectors  1  and  2.  Both 
intensity  patterns  have  the  same  angular  frequency  but  independent  phases.  Part 
of  the  difference  in  phases  is  due  to  the  NR  phase  shift,  At}>NR,  the  rest  a  conse¬ 
quence  of  optical  path  differences  between  the  various  combining  waves.  A<j>OPL. 
For  convenience  we  may  separate  these  phases  and  consider  the  intensitv  patterns 
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to  be  described  by  .4,  s'\n(6u:t  -  A(t>OPL;  and  U  —  A  O  R;  The  combination 

of  these  two  signals  gives  another  sinusoidal  function  with  amplitude  .47-.  The  total 
phase  difference  is  related  to  the  amplitudes  of  these  three  signals  bv 


~  ^doPL  —  cos 


-1  -4  r  ~  1  -^1  ~  ••*?) 
2.4,. 4  2 


(33) 


We  can  eliminate  Ad>OP[j  by  changing  control  parameters  and  making  a  second 
measurement.  We  are  left  with  the  change  in  AtpNR.  A  bit  more  algebra  is  required 
to  extract  the  transmissions  from  the  amplitudes. 


3.3.5  Experimental  attempts 

Our  LiNb03  samples  proved  very  difficult  to  work  with.  Previous  research  with 
crystal  samples  in  the  group  was  plagued  by  noise  and  non-repeatability  [91].  The 
noisiness  of  signals  hampered  our  effects  as  well.  Of  the  dozen  or  so  crystals  of 
LiNb03  available,  only  one  was  found  with  the  proper  cut — polished  faces  normal  to 
the  c  axis.  No  BinSi02o  crystals  had  the  proper  orientation.  The  cample  had  a  large 
interaction  length,  1.3  cm,  and  the  crystal  was  used  in  an  open  circuit  configuration. 
A  set-up  was  similar  to  Fig  14  except  the  detectors  just  measured  the  transmitted 
signal — the  mirrors  were  blocked  off  and  no  moving  fringes  were  produced  at  the 
detectors.  Moderate  intensity  (10-30  mW/cm2)  input  beams  produced  independent, 
noisy,  transmitted  signals.  The  fringe  velocity  was  limited  to  no  slower  than  5  x  10~9 
m/s  and  for  this  value,  little  power  transfer  was  expected.  A  new  attempt  was  made 
mixing  unexpanded  beams.  Fig  15  shows  typical  simultaneous  readings.  Coupling 
between  the  beams  is  evident  in  that  one  trace  has  peaks  where  the  other  trace  has 
depressions  and  vice  versa.  The  noise  was  not  due  to  any  fluctuations  in  the  laser 
output.  Removing  the  crystal  produced  flat,  steady  signals  at  the  detectors. 

We  borrowed  a  second  sample  of  LiNb03  with  polished  c  faces  from  the  Univer¬ 
sity  of  Osnabruk.  This  had  a  large  photovoltaic  field  (35  kV/cm)  but  a  very  small 
interaction  length  (2  mm),  and  produced  similar  noisy  traces  with  hints  of  beam 
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Figure  15:  Typical  traces  for  transmitted  beam  intensities.  The  signals  are  very 
noisy,  but  beam  coupling  is  evident. 

One  cannot  make  sweeping  generalizations  from  looking  at  just  two  samples, 
but  obtaining  poor  results  from  two  crystals  of  very  different  origins  bodes  ill  for 
the  use  of  LiNb03  in  nonreciprocal  applications. 

3.3.6  Limitations 

Numerical  methods,  particularly  ones  involving  iteration,  are  prone  to  instability. 
Our  treatment  was  no  exception.  For  values  of  Ev  greater  than  15  kV/cm  and 
positive  fringe  velocities,  the  computer  program  would  not  converge  to  a  particular 
solution.  Instead  it  seemed  to  oscillate  between  two  solutions.  Nonlinear  materials 
often  exhibit  bistable  behaviour,  and  the  mathematical  model  may  be  pointing  to 
a  region  where  LiNb03  is  optically  bistable — two  possible  intensity  states  for  the 
same  input. 

The  presence  of  absorption  dramatically  slowed  down  convergence,  and  scores 
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Figure  16:  NR  phase  shift  vs  fringe  velocity  for  various  values  of  absorption. 


of  iterations  typically  were  required.  This  is  likely  due  to  the  poor  initial  guess 
calculated  from  the  numerical  treatment  of  Ref  [59].  But  absorption  has  been 
studied,  and  Fig  16  shows  that  a  healthy  attenuation  constant  of  2  cm-1  reduces 
the  maximum  NR  phase  shift  by  a  factor  of  2.  The  conditions  are  the  same  for  Fig 
13a  with  Ev  -  10  kV/cm. 

The  question  of  modulation  depth  is  important  here.  The  solutions  to  the  in¬ 
tensity  coupled  wave  equations  shown  in  Figs  4a,  5a  and  6a  reveal  there  are  regions 
in  the  crystal  in  which  m  is  not  much  smaller  than  1.  Indeed,  for  the  Bi12SiOJ0 
plot,  the  intensity  solutions  cross,  so  there  is  one  point  where  m  =  1.  As  suggested 
by  Ref  [8],  we  may  expect  a  reduction  in  photorefractive  coupling  for  higher  mod¬ 
ulation  depths.  Taking  this  into  account  may  lower  the  expected  magnitude  of  NR 
effects. 
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3.4  Conclusions 

This  chapter  began  by  posing  the  question:  Can  photorefractive  materials  be  used 
in  the  design  of  a  nonreciprocal  phase  shifter?  The  answer  is  yes,  but  there  are 
several  drawbacks. 

The  effect  is  small,  several  times  less  than  that  obtainable  from  magnetooptic 
NR  phase  shifters.  For  use  in  a  fibre-optic  ring  laser  gyroscope,  this  might  be 
adequate,  for  the  round-trip  optical  path  difference  between  counterpropagating 
beams  is  typically  much  less  than  a  wavelength  [22].  However  there  is  a  concurrent 
effect  which  might  be  undesirable  in  a  resonator — power  is  transferred  between  the 
beams.  Decoupling  these  effects  would  be  difficult.  The  best  control  parameter, 
fringe  velocity,  requires  delicate  and  precise  movement  of  either  a  mirror  or  the 
crystal  itself.  Of  the  three  crystal  types  studied,  the  material  with  the  largest 
theoretical  NR  behaviour,  LiNb03,  has  produced  poor  experimental  results. 

Our  research  has  shown,  however,  that  NR  transmission  effects  in  photorefrac¬ 
tive  materials  have  great  potential.  It  has  been  suggested  to  use  a  photorefractive 
crystal  as  an  optical  diode,  but  we  have  shown  that  using  fringe  velocity  one  can 
control  the  strength  of  the  beam  coupling,  reduce  it  to  negligible  coupling,  or  even 
change  the  direction  of  power  transfer. 

An  interesting  byproduct  of  this  study  is  a  numerical  method  for  treating  two- 
wave  mixing  with  spatially  variable  coupling.  The  techniques  described  should  be 
readily  adaptable  to  other  geometries. 


Chapter  4 

Unidirectional  Ring  Resonator 

One  ring  to  bring  them  all  and  in  the  darkness  bind  them. 

J.R.R.  Tolkein 
(a  fellow  Mertonian) 

Beam  coupling  in  photorefractive  materials  often  results  in  power  transfer.  As 
discussed  in  Chaptet  2,  this  light  intensity  interchange  is  described  by  a  gain  co¬ 
efficient,  2  T  sin  which  has  a  maximum  in  Bii2SiO20  with  moving  fringes  when 
$g  =  90°.  A  photorefractive  crystal  with  high  gain  can  be  used  as  the  active  medium 
in  a  resonator.  Oscillation  will  occur  as  long  as  the  photorefractive  gain  exceeds 
losses  due  to  absorption,  diffraction  and  reflection  from  optical  elements  in  the  cav¬ 
ity.  We  will  see  that  the  oscillating  beam  chooses  a  path  whose  length  satisfies  the 
phase  condition  for  resonance — namely  one  roundtrip  represents  an  accumulated 
phase  which  is  an  integer  multiple  of  27T. 

Optical  resonators  with  photorefractive  gain  are  the  subject  of  keen  interest, 
particularly  in  concert  with  phase  conjugating  mirrors  [28].  The  interest  is  un¬ 
derstandable — phase  conjugating  mirrors  retroreflect  incident  light,  reversing  the 
phase,  so  any  aberrating  effects  in  the  resonator  cavity  may  be  corrected,  and  the 
reflectivity  can  exceed  unity.  (Of  course  there  are  different  theoretical  considerations 
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in  a  system  which,  in  effect,  folds  the  optical  path.  For  example,  the  concept 
of  cavity  length  loses  its  meaning.)  Many  of  the  applications  of  photorefractive 
crystals  discussed  in  Chapter  1,  such  as  in  optical  computing,  image  processing, 
and  associative  memories,  utilise  such  optical  resonators. 

Most  of  the  work  with  these  resonators  uses  BaTi03,  and  many  configurations 
are  discussed  in  Ref  [131) .  We  were  interested  in  seeing  if  one  design  in  particular, 
the  linear  resonator,  would  work  using  Bii2SiO20  as  the  gain  medium.  (The  linear 
resonator  is  discussed  in  greater  detail  in  Chapter  5.) 

A  good  starting  place  for  the  understanding  of  resonator  operation  with  pho- 
torefractives  is  the  unidirectional  ring  resonator  (URR).  In  this  tried-and-true  ar¬ 
rangement,  a  photorefractive  crystal  is  placed  in  a  leg  of  a  resonator  formed  by 
three  or  more  mirrors  (see  Fig  1).  A  single  input  beam  incident  on  the  crystal 
generates  scatter,  some  of  which  may  be  captured  in  the  resonator  and  amplified. 
This  arrangement  has  been  explored  using  BaTi03  by  White,  et  al.  [131],  Kwong, 
et  al.  [79] ,  and  Fischer  [38].  BaTi03  has  very  high  gain  due  to  its  large  electrooptic 
effect  and  one  can  induce  oscillation  at  very  low  light  intensities  (%  50/xW/cm2  at 
633  nm),  and  may  even  achieve  bidirectional  operation.  Rajbenbach  and  Huignard 
reported  a  study  of  the  URR  using  Bi12SiO20[105].  In  a  material  like  Bi12SiO20  wave 
propagation  in  the  resonator  is  unidirectional  because  the  applied  electric  field  en¬ 
hances  the  gain  for  beams  in  the  transmission  geometry,  but  not  in  the  reflection 
geometry. 

A  URR  is  relatively  easy  to  set  up  (as  much  as  any  resonator  can  be),  oscillation 
is  self-starting  (no  need  for  a  seed  beam),  and  the  gain  can  be  quite  large.  We  felt 
the  experience  gained  in  working  with  the  URR  would  serve  us  well  later. 

We  start  this  chapter  with  a  short  review  of  the  pertinent  theory  of  two-wave 
mixing  and  derive  an  analytic  solution  to  the  intensity  coupled  wave  equations  in 
the  presence  of  absorption.  Then  we  consider  the  feedback  of  the  amplified  beam  in 
the  ring  resonator  geometry.  Next  we  discuss  experimental  results.  Although  our 
work  reproduced  many  of  the  findings  of  Rajbenbach  and  Huignard,  we  concentrate 
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on  new  observations  or  quantify  some  of  their  qualitative  results.  In  particular  we 
discuss  the  appearance  of  diverse  higher  diffracted  orders. 

The  thrust  of  this  and  the  next  two  chapters  is  experimental,  and  they  follow  a 
particular  sequence  of  experiments.  The  fruition  of  our  labours  with  the  URR  comes 
to  bear  in  Chapter  5,  in  which  we  consider  a  new  gain  mechanism  for  photorefractive 
materials  [110].  This  mechanism  is  different  from  that  discussed  in  Chapter  2  (that 
one  depended  on  the  imaginary  part  of  the  space  charge  field).  This  new  gain 
process  is  very  similar  to  one  proposed  recently  in  which  higher  diffracted  orders 
may  amplify  a  weak  mixing  wave  [6].  We  employ  this  new  mechanism  successfully  in 
a  ring  and  linear  resonator.  Finally  in  Chapter  6  we  look  more  closely  at  the  higher 
diffracted  orders  discussed  at  the  end  of  this  chapter  with  a  special  eye  towards 
applications. 


4.1  Background 

4.1.1  Optimum  frequency  detuning  and  grating  spacing 

In  two-wave  mixing  the  combination  of  small  interbeam  angles  and  a  large  applied 

electric  field  leads  to  drift  dominated  behaviour.  In  this  regime,  except  for  very  small 

interbeam  angles  or  weak  applied  fields,  Eq  E,0  Em  Et.  Our  expression  for 

the  modulated  space  charge  field  in  this  limit  is  given  by  Eqn  2.55.  The  imaginary 

part  of  E, i  has  extrema  for  b  =  ±Em/E,o-  The  maximum  value  occurs  for  6opt  = 

-Em/E,0  and  the  positive  root,  6mjn,  gives  the  minimum  value.  For  6opt 

E,o2  Eq  /1X 

E,x=jm  ^~T+  E'  £  (*) 

in  other  words,  pure  imaginary  so  =  90°.  For  6m;n 


E, i  =  —  \m  E,, 


Here  $g  =  0.  This  positive  root  is  not  discussed  in  the  literature,  but  it  gives  insight 
into  the  resonance  condition  for  6.  The  ratio  Em/ E,0  is  small,  therefore  there  is 
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little  separation  of  6opt  and  6min— the  resonance  is  sharp. 

For  fixed  £m  and  E. 0  (ie.  constant  interbeam  angle  and  applied  field),  f>opt  is 
constant.  But  6Cpt  =  ^optri,  where  <*>u;opt  is  the  optimum  detuning  frequency  for 
the  fixed  grating  spacing.  Using  the  definition  of  Tj  and  from  Chapter  2  we 


obtain 


r  t  e^sND  T  esND  r 

°pt  "  bopt  6.7 rN;  /o  -  t.KE, o  /o 


The  optimum  detuning  frequency  is  linearly  dependent  upon  the  total  intensity. 

Both  Em  and  Eq  are  proportional  to  the  grating  spacing,  A.  We  can  optimize 
the  imaginary  part  of  E,\  in  Eqn  1  for  a  particular  grating  spacing.  The  condition 
to  be  satisfied  is  E,02  =  Em  E„,  which  leads  to 


Aoot  — 


2nE,0  I  ey. 


For  the  crystal  parameters  in  Table  2.1,  Aopt  =  15/xm.  Any  grating  spacing  close 
to  this  value  will  insure  high  gain  at  the  optimum  detuning  frequency. 

We  have  considered  conditions  favourable  to  high  gain.  Now  we  need  to  consider 
losses. 


4.1.2  Coupled  intensity  equations  with  absorption 

The  results  of  Chapter  2  show  that  for  waves  mixing  at  the  optimum  detuning 
frequency,  $9  =  90°.  That  is,  there  is  maximum  intensity  coupling  (rsin$0  = 
T),  and  no  phase  coupling.  Consider  two  waves  mixing  as  in  Fig  2.5  at  optimum 
detuning.  We  make  the  assumption  that  absorption  losses  are  sufficiently  small  that 
T  is  constant  through  the  medium.  Then  our  coupled  intensity  equations  are 


dU 

dx 

dl_ 

dx 


=  2rT7TT-^- 

=  -2r  UI- 
1+  +  /_ 


(5) 

(6) 


where  a  is  the  absorption  constant. 
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4.1.3  Ring  resonator 

We  may  now  consider  two  wave  mixing  in  which  the  amplified,  oscillating  beam 
undergoes  feedback  as  shown  in  Fig  1.  Rajbenbach  and  Huignard  have  shown  that 
in  BijjSiOjo  the  oscillating  beam  can  arise  from  the  noise  spectrum  of  a  pump  beam 
[105].  In  other  words,  oscillation  is  self-starting.  The  particular  frequency  compo¬ 
nent  of  the  noise  spectrum  corresponding  to  the  optimum  frequency  detuning  will 
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achieve  the  highest  gain,  thus  the  frequency  of  the  oscillating  beam  is  automati¬ 
cally  detuned  to  the  optimum  value.  This  means  the  intensity  expression  of  Eqn  'J 
is  valid. 

At  steady-state,  the  condition  for  oscillation  is  that  both  the  amplitude  and 
phase  of  E±  return  to  their  original  values  after  one  complete  circuit.  The  amplitude 
condition  can  be  expressed  in  terms  of  the  intensity 

MO  )--RL(d)  (11) 

where  R  is  the  combined  reflectivity  of  the  resonator  for  one  round  trip.  We  ignore 
diffraction  losses  in  this  treatment.  Substituting  Eqn  9  at  x  =  d  into  the  above 
expression  vields 

740)  =  1  -  _ 

Pr°  1.(0)  1  -  Re-** 

For  oscillation  to  occur,  this  ratio  must  be  greater  than  0,  so  T  must  satisfy 

2Td  >  ad  —  In  R  (13) 

Eqn  12  gives  the  steady-state  beam  ratio  of  the  resonator. 

The  phase  condition  for  oscillation  is  that  the  phase  delay  of  the  round  trip 
must  be  an  integer  multiple  of  27T.  This  is  not  a  very  restrictive  condition — there 
are  many  components  of  the  noise  spectrum  of  the  input  pump  beam  which  can 
find  optical  paths  around  the  resonator  satisfying  the  phase  requirement. 

In  Bii2Si02o  oscillation  in  a  ring  resonator  is  unidirectional.  For  the  set-up  in 
Fig  1  we  take  the  crystal  orientation  to  give  efficient  power  transfer  for  transmission 
geometry.  A  counterpropagating  beam  would  mix  with  the  pump  beam  to  form  a 
reflection  grating  with  K  in  the  i  direction.  But  the  applied  field  (necessary  for  high 
gain  in  Bii2Si02o)  is  in  the  z  direction  and  cannot  enhance  the  reflection  coupling. 
Thus  the  gain  for  the  reflection  case  will  not  be  large  enough  to  sustain  oscillation. 
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mirror 

Figure  2:  Experimental  layout. 


4.2  Experimental  Results 

4.2.1  Layout  and  resonator  alignment 

The  layout  of  the  URR  we  studied  is  shown  in  Fig  2.  The  resonator  cavity  was 
formed  by  3  dielectric  plane  mirrors.  One  mirror  (C)  was  mounted  on  a  piezoelectric 
driven  translation  stage,  and  by  applying  a  voltage  to  the  driver,  the  geometry  of 
the  cavity  could  be  slightly  changed.  The  Bi12Si02o  crystal  sample  had  5mm  x 
5mm  front  and  rear  faces  and  a  length  of  5.33mm.  The  crystal  had  silver  painted 
electrodes  on  the  001  face,  and  the  orientation  used  was  as  in  Fig  2.5.  The  crystal 
was  mounted  in  the  holder  described  in  John  M.  Heaton’s  D.Phil.  thesis  [49] .  This 
holder  has  plunger  type  contacts  and  was  connected  to  a  Brandenburg  HV  power 
supply.  Unless  otherwise  stated,  the  applied  voltage  was  4.5  kV.  Reflection  off  a 
glass  plate  inserted  in  one  leg  of  the  resonator  provided  a  sample  of  the  oscillating 
beam.  The  strength  of  the  oscillating  signal  was  measured  using  a  photodiode.  An 
iris  aperture  allowed  only  low  order  modes  to  oscillate  in  the  resonator  when  it  was 
stopped  down  (vi  iable  down  to  slightly  less  than  1  mm  in  diameter). 
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The  pump  beam  was  the  expanded,  spatially  filtered  and  collimated  output 
of  an  argon  ion  laser.  The  wavelength  was  514  nm  and  the  light  was  vertically 
polarised.  Part  of  this  beam  was  split  off  to  provide  a  normalising  signal,  as  well  as 
a  reference  beam  for  interferometric  measurements. 

The  resonator  alignment  is  critical  to  achieve  oscillation.  A  systematic  place¬ 
ment  of  the  mirrors  using  the  following  method  worked  well.  The  resonator  was 
built  up  “backwards”,  using  an  injected  beam  from  a  He-Ne  laser  source.  The 
He-Ne  laser  was  put  in  the  position  of  the  photodiode  in  Fig  2.  The  He-Ne  beam 
was  parallel  to  the  optical  table  and  the  beam  height  was  level  with  the  centre  of 
the  crystal.  Mirrors  A,B,  and  C  were  positioned  sequentially,  forming  the  general 
arrangement  of  the  ring  resonator.  The  line  AB  was  at  an  angle  of  %  1.9°  to  the 
pump  beam  to  give  a  grating  spacing  of  15  /tm.  The  height  at  which  the  He-Ne 
beam  struck  each  mirror  was  constant.  Any  uphill  or  downhill  slope  to  the  align¬ 
ment  beam  was  detected  easily  using  a  stopped  down  iris  aperture  of  the  correct 
height.  Care  was  taken  to  insure  the  reflection  off  mirror  C  hit  the  glass  plate  at 
the  exact  spot  as  the  injection  of  the  be;  *n  from  the  He-Ne  laser.  The  iris  aperture 
was  placed  downstream  of  the  glass  plate  now.  Any  misalignment  resulted  in  the 
beam  reflected  off  mirror  C  missing  the  aperture.  Each  mirror  mount  had  vertical 
and  horizontal  tilt  controls  and  these  were  adjusted  until  any  reflected  beam  re¬ 
traced  the  path  of  the  original  alignment  beam.  When  this  was  accomplished,  the 
resonator  was  sufficiently  lined  up  to  obtain  some  sort  of  oscillation. 

Additional  fine-tuning  of  the  resonator  was  made  in  two  stages.  First  the  iris 
aperture  was  opened  to  its  maximum  diameter.  The  pump  beam  was  turned  on, 
a  voltage  applied,  and  oscillation  started.  The  oscillating  beam  consisted  of  a 
writhing,  complicated,  superimposed  mode  structure.  We  replaced  the  He-Ne  laser 
with  the  photodiode  and  mirror  C  was  adjusted  to  optimize  the  signal  strength  of 

t 

the  oscillating  beam  sample.  (We  found,  by  trial  and  error,  that  this  signal  was 
most  sensitive  to  changes  in  the  tilt  of  mirror  C,  off  which  the  oscillating  beam  had 
the  largest  angle  of  incidence.)  The  second  stage  was  to  close  the  iris  aperture  to 


! 
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its  minimum  diameter  and  then  adjust  the  mirrors  to  optimize  the  symmetry  in  the 
lower  order  modes.  For  this  tuning,  mirror  B  worked  best. 

4.2.2  Comparison  with  previously  reported  results 

Oscillation  was  self-starting  in  our  ring  resonator.  As  mentioned  above  the  aperture 
could  be  used  to  select  single  modes.  Not  a  specific  lower  order  mode,  but  one  of 
many.  A  sample  ot  these  modes,  as  well  as  an  example  of  the  structure  with  a  fully 
open  aperture,  is  given  in  Fig  3.  There  are  two  images  in  each  photograph  because 
there  were  reflections  off  both  the  front  and  rear  surfaces  of  the  glass  plate.  Here 
the  photodiode  in  Fig  2  was  replaced  by  a  lens  and  a  camera  back. 

The  main  difference  between  these  modes  and  those  reported  by  Rajbenbach  and 
Huignard  is  the  obvious  rectangular  symmetry.  Ref  i  1 05]  reported  modes  exhibiting 
circular  symmetry.  This  means  in  our  resonator  the  limiting  aperture  must  be  the 
crystal  itself. 

Mode  stability  was  fairly  good — from  a  few  seconds  to  two  minutes.  High  fre¬ 
quency  vibration,  such  as  produced  by  rapping  on  the  optical  table  did  not  initiate 
mode  transition,  but  low  frequency  disturbances,  such  as  caused  by  air  currents, 
did.  Often  a  particular  oscillating  mode  would  fade  and  there  would  be  a  delay 
of  some  tens  of  seconds  before  oscillation  would  build  up  again.  At  other  times  a 
single  stable  mode  would  give  way  to  a  waver  of  competing  modes.  The  struggle 
would  have  a  single  stable  mode  emerge  as  the  victor.  This  competition  occured 
most  frequently  between  modes  TEMoo  and  TEMn,  or  between  modes  TEMio 
and  TEM02.  Ref  [25]  describes  mode  competition  in  a  URR  using  Bi12Si02o-  In 
particular,  by  placing  lenses  in  the  resonator,  the  authors  were  able  to  emulate  a 
confocal  resonator  and  limit  the  number  of  higher  order  modes  by  decreasing  the 
Fresnel  number  of  the  system. 

We  used  low  input  intensities,  less  than  15  mW/cm2.  This  is  an  order  of  mag¬ 
nitude  less  than  that  used  by  Rajbenbach  and  Huignard.  One  reason  we  chose  to 
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Figure  4:  Measurement  of  detuning  frequency. 

work  in  this  range  was  necessity.  Our  laser  output  was  less  stable  at  high  power 
settings.  A  second  reason  was  that  it  was  more  convenient  when  performing  inter¬ 
ferometric  time  studies.  The  frequencies  we  measured  were  on  the  order  of  hertz  to 
tens  of  hertz  instead  of  the  hundreds  of  hertz  characteristic  of  high  light  intensities. 

The  interferometric  arrangement  is  shown  in  Fig  4.  The  reflection  of  the  oscil¬ 
lating  beam  off  the  front  face  (%  20%)  is  combined  with  part  of  the  input  beam  in  a 
Mach-Zehnder  scheme.  Any  frequency  difference  between  the  two  beams  results  in 
an  interference  pattern  of  fringes  moving  at  the  detuning  frequency.  This  frequency 
can  be  measured  by  analysis  of  an  oscilloscope  trace. 

We  had  a  number  of  questions  about  the  detuning  frequency  we  wanted  to 
answer. 

1.  Does  the  detuning  frequency  depend  on  mode? 

2.  Is  the  behaviour  linear  with  input  intensity,  as  suggested  by  Eqn  3? 

3.  How  will  changing  the  length  of  the  cavity  affect  th»  d«*tunir.g  frequency? 
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Figure  6:  Change  in  cavity  length  due  to  displacement  of  mirror  C.  For  our  case: 
AB  =  65.2  cm,  BC  —  33.0  cm,  C A  =  56.5  cm. 

Fig  5  shows  the  frequency  difference,  6 u  =  6u/2ir,  for  various  modes  as  a  func¬ 
tion  of  input  pump  intensity.  Each  data  point  is  a  statistical  average  of  up  to  ten 
measurements  taken  for  a  particular  mode  and  a  given  intensity.  Error  bars  are  not 
shown  in  this  figure,  but  typical  error  bars  for  these  interferometric  measurements 
appear  in  Fig  7a.  The  scatter  is  somewhat  broad,  but  there  is  no  discernible  trend 
difference  between  the  modes.  Within  experimental  error  every  mode  seems  to  have 
roughly  the  same  detuning  frequency  for  the  same  intensity.  The  dependence  of  8v 
on  the  input  intensity  appears  linear  in  this  region.  The  straight  line  is  based  on 
a  linear  regression  fit  to  the  function  6v  —  cIC:  where  c  is  the  slope.  We  may  take 
Io  =  /_(0)  because  the  oscillating  beam  intensity  is  very  weak  compared  to  the 
pump  intensity.  The  best  fi*  slope  is  0  23  rr,2/Ws.  This  is  ovei  eight  times  larger 
than  the  slope  predicted  by  Eqn  3  using  the  crystal  parameters  of  Table  2.1.  Two 
factors  may  account  for  this  discrepancy.  The  table  parameter  values  were  chosen 
as  representative  of  the  literature.  There  may  be  significant  variation  among  crystal 
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samples.  In  particular,  this  Bi , 2Si O 30  sample  produced  photocurrent  densities  (pro¬ 
portional  to  sNp)  nearly  four  times  larger  than  those  for  other  BiuSiO^o  samples  for 
a  given  input  intensity  [31].  The  second  factor  concerns  the  size  of  the  applied  field. 
Some  authors  report  that  the  applied  electric  field  in  Bi!2SiO:>o  is  often  less  than 
that  expected  from  the  applied  voltage  jie.  E,q  •-  V/d).  For  example.  Jonathan, 
et  al.,  found  their  applied  field  was  half  as  large  as  expected  :64l.  This  decrease  in 
field  strength  may  be  due  to  shadowing  effects  at  the  electrodes,  or  a  loss  due  to 
poor  ohmic  contact  of  the  electrode  to  the  crystal.  A  combination  of  higher  values 
for  sNp  and  a  smaller  E, 0  would  result  in  a  larger  slope  in  Eqn  3. 

In  BaTiOa  one  can  change  the  detuning  frequency  by  displacing  one  resonator 
mirror  slightly  [79],  This  changes  the  cavity  length  and  the  difference  in  detuning 
is  necessary  to  satisfy  the  phase  condition  for  oscillation.  The  detuning  frequency 
can  take  on  positive  or  negative  values,  and  is  periodic  with  the  mirror  shift. 

We  would  not  expect  a  change  in  cavity  length  to  produce  such  a  large  range 
of  detuning  frequencies  in  Bil2SiOjo  because  the  gain  is  very  strong  only  for  values 
close  to  5u/opt.  It  is  likely  that  oscillation  would  cease  before  a  large  change  in 
detuning  took  place. 

Fig  6  shows  the  change  in  geometry  when  mirror  C  is  displaced  a  distance  61 
along  the  bisector  of  lACB.  Both  sides  AC'  and  BC'  are  shortened  by  a  length 
81  cos  The  total  change  in  cavity  length  is 

AL  =  -261  cos  \9C  (H) 

For  our  case  9c  =  89.5°.  The  piezoelectric  translator  had  an  elongation  of  1  pm 
/ 1 00  V.  Thus  applying  10  V  gave  SI  =  0.1pm.  This  resulted  in  AL  =  0.2SA,  where 
A  =  514  nm.  This  change  in  length  was  confirmed  by  injecting  a  beam  into  the 
resonator  with  no  applied  field  to  the  crystal  (no  beam  coupling)  and  setting  up  an 
interferometer  as  in  Fig  4.  An  applied  voltage  of  40  V  to  the  piezoelectric  translator 
shifted  the  interference  pattern  by  almost  exactly  one  fringe. 

Fig  7a  shows  the  results  of  mirror  displacement  on  the  detuning  frequency  of 
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mode  TEMi0  for  an  incident  intensity  of  9.7  m\V  cm2.  The  error  bars  represent 
one  standard  deviation  of  the  experimental  values.  Fig  7b  includes  various  modes. 
Again  the  different  modes  all  seem  to  be  at  roughly  the  same  detuning  frequency 
and  this  is  independent  of  mirror  displacement. 

The  ordinate  axis  in  Figs  7a  and  b  has  quotation  marks  around  the  word 
“change'’.  Did  we  really  change  the  cavity  length  by  displacing  mirror  C?  As 
soon  as  the  voltage  to  the  piezoelectric  translator  was  modified,  the  oscillation  in 
the  URR  immediately  stopped.  After  several  seconds  to  a  few  minutes,  the  oscil¬ 
lating  beam  would  reappear — equally  intense  and  at  the  same  detuning  frequency. 
It  is  likely  that  after  the  original  oscillation  was  disturbed  by  shifting  mirror  C, 
the  new  oscillating  beam  travels  along  a  slightly  different  path  which  satisfies  the 
twin  conditions  of  optimum  detuning  and  proper  cavity  length.  We  still  obtained 
oscillation  when  the  path  choices  were  restricted  by  putting  in  one  or  two  small 
apertures  in  the  system.  In  this  regard  oscillation  in  the  URR  is  very  forgiving. 

One  benefit  of  using  low  incident  power  was  the  discovery  of  a  threshold  for 
oscillation.  No  oscillation  occured  for  pump  intensities  less  than  0.8  mW/cm2.  This 
is  seen  in  Fig  8.  Here  the  ratio,  (3r ,  of  the  intensity  of  the  oscillating  beam  ( I+(d )) 
to  the  pump  intensity  (/_(0))  is  plotted  as  a  function  of  input  intensity.  According 
to  Eqns  9  and  10,  (3r  =  R,0j R  and  thus  depends  on  T,R,  and  a,  all  of  which 
are  independent  of  intensity.  We  should  mention  that  Rajbenbach  has  reported 
a  region  of  reduced  gain  in  Bi12SiO20,  lor  nearly  degenerate  two-wave  mixing  at 
similar  intensities  [106].  It  is  possible  that  at  low  light  intensities  not  all  acceptor 
sites  are  filled,  ie.  is  reduced,  which  leads  to  reduced  coupling. 

Threshold  devices  are  increasingly  important  in  many  optical  applications  such 
as  image  processing  [70]  and  associative  memories  [112].  A  URR  using  Bi12SiO20 
could  be  used  as  a  threshold  detector,  however  the  mode  structure  of  the  oscillating 
beam  would  be  a  limiting  factor  on  the  spatial  resolution  of  the  output. 

The  highest  beam  ratio  we  measured  was  0.08.  This  is  lower  than  the  steady- 
state  beam  ratio  of  0.6  reported  in  Ref  [107].  However,  Rajbenbach  and  Huignard 
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Figure  7:  Detuning  frequency,  6 v,  vs  “change”  in  cavity  length,  that  is,  a  mirror  dis¬ 
placement  designed  to  give  this  change  in  length,  (a)  TEMi0.  Error  bars  represent 
one  standard  deviation. 
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Figure  7:  Detuning  frequency,  6v,  vs  “change”  in  cavity  length,  that  is,  a  mirror 
displacement  designed  to  give  this  change  in  length,  (b)  Various  modes. 
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Figure  8:  Ratio  of  oscillating  intensity  to  input  intensity,  (3r,  vs  input  intensity.  No 
oscillation  occurred  for  an  input  less  than  0.8  mW/cm2. 
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used  a  ionger  wavelength  (568  nm)  to  reduce  the  absorption  loss,  and  had  a  longer 
interaction  length  (1  cm).  We  can  use  our  value  and  Eqn  10  to  estimate  the  gain 
coefficient  in  our  system.  For  this  crystal  a  =  1.1  cm’1,  d  =  0.53  cm.  and  R  =  0.65. 
(To  calculate  R:  Mirrors  .4  and  B  had  reflectivities  of  0.95,  C  had  a  reflectivity  of 
0.98,  the  transmission  of  the  glass  plate  was  0.92.  and  there  was  a  20%  reflection 
loss  off  the  front  face  of  the  crystal.)  In  terms  of  3r 

r=-^n(R(e-°'{l,R3r)-3r))  (15) 

For  our  experimental  values.  F  =  1.05  cm-1.  This  is  lower  (by  a  factor  of  ten) 
than  that  calculated  from  the  treatment  of  Secs  2.1.8  and  2.3.2.  but  those  did  not 
take  into  account  optical  activity  or  higher  diffracted  orders.  Further,  as  we  have 
discussed,  E, 0  may  be  reduced. 

4.3  Higher  Diffracted  Orders 

4.3.1  Emergence  of  higher  diffracted  orders 

In  the  derivation  of  the  coupled  wave  equations  in  Chapter  2  we  neglected  coupling 
to  higher  diffracted  orders  (see  Fig  2.6).  We  argued  there  that  this  simplification  was 
valid  if  the  interbeam  angle  was  large.  However,  to  get  the  high  gain  necessary  for 
oscillation  in  the  URR  using  Bii2SiO20  the  interbeam  angle  must  be  small.  Higher 
diffracted  orders  (HDOs)  were  visible  whenever  the  oscillating  beam  was  present. 
In  general,  besides  the  pump  beam  (0  order)  and  the  oscillating  beam  (-1  order), 
the  4-1,  -2,  and  -3  orders  were  seen  (see  Fig  9a). 

Au  and  Solymar  have  described  the  emergence  of  higher  diffracted  orders  and 
have  shown  that  significant  power  transfer  to  these  orders  is  possible  [7].  They  have 
also  identified  an  off-Bragg  coupling  mechanism  by  which  power  transfer  between 
the  principal  mixing  waves  may  occur,  even  when  the  phase  mismatch  angle, 
equals  0  (6].  We  will  not  go  into  the  full  theory  in  this  thesis,  but  an  off-Bragg 
coupling  mechanism  similar  to  that  presented  in  Ref  [6]  is  described  in  Chapter  5. 
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Figure  9:  Two  species  of  higher  diffracted  orders  in  the  URR.  (a)  HDOs.  (b) 
Companion  orders. 
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The  kev  feature  here  is  that  conditions  for  HDOs  improve  as  the  interbeam  angle 
becomes  smaller. 

4.3.2  Companion  orders 

The  reflection  of  the  oscillating  beam  off  the  glass  plate  revealed  an  interesting 
phenomenon.  The  oscillating  beam  had  an  entourage.  A  number  of  other  beams 
(which  we  shall  cali  companion  orders)  were  visible.  These  represented  a  different 
“species"  of  higher  diffracted  orders.  They  appeared  on  either  side  of  the  oscillating 
beam  and  were  evenly  spaced  along  a  line  some  35°  out  of  the  piane  of  the  resonator. 

Closer  inspection  revealed  that  companion  orders  also  appeared  around  the 
pump  beam,  the  HDOs,  and,  surprisingly,  the  reflection  of  the  oscillating  beam 
off  the  front  crystal  face.  (We  shall  call  the  pump  and  oscillating  beams  and  the 
HDOs  due  to  their  two-wave  mixing  the  principal  beams ,  in  relation  to  the  compan¬ 
ion  orders.)  Several  views  of  these  companion  orders  are  seen  in  Fig  10. 

Several  observations  may  be  made  here: 

•  The  angular  separation  of  companion  orders  was  much  smaller  than  that  for 
the  HDOs.  (We  will  measure  this  separation  in  the  next  subsection.) 

•  There  seemed  to  be  a  focusing  effect  in  the  companion  orders — as  one  moved 
away  from  the  principal  beam,  on  one  side  the  orders  became  increasingly 
focused  while  on  the  other  side  the  companion  orders  were  increasingly  defo- 
cused. 

•  The  companion  orders  shared  the  same  mode  structure  as  the  principal  beam. 

•  The  spacing  and  orientation  of  the  companion  orders  were  not  affected  by  a 
decrease  in  mterbeam  angle  between  the  pump  and  oscillating  beam  to  1.44°. 

•  The  number  of  companion  orders  around  each  principal  beam  diminished  as 
the  order  of  the  HDO  got  larger  (i.e.  as  one  moved  away  from  the  pump 
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Figure  10:  Companion  Orders  in  the  Ring  Resonator,  (a)-(c)  Oscillating  beam, 
(d)  Pump  beam  (one  side  only),  (e)  Higher  diffracted  -2  order,  (f)  Reflection  of 
oscillating  beam  off  front  face.  In  all  but  (d),  the  principal  beam  is  the  second  from 
the  bottom. 
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beam).  The  pump  beam  usually  had  eight  or  nine  companion  orders,  the 
oscillating  beam  five  to  seven.  IIDOs  -2  and  —1  three  or  four  and  HDO  5 
two  or  three.  It  was  not  unusual  to  find  that  a  single  input  to  the  crystal, 
the  pump  beam,  produced  over  30  separate  beams.  This  was  a  very  prolific 
crystal! 

How  do  the  companion  orders  arise?  It  should  be  emphasized  that  because 
of  the  iris  aperture  and  the  angular  separation  of  the  companion  orders,  only  the 
oscillating  beam  could  circumnavigate  the  resonator  and  be  reinforced.  Only  the 
oscillating  and  pump  beam  entered  the  front  face  of  the  crystal.  The  companion 
orders  were  not  due  to  any  reflection  external  to  the  crystal,  they  had  to  be  produced 
inside  the  crystal. 

One  possibility  we  considered  was  that  the  companion  orders  were  some  type  of 
preferential,  amplified  scatter  from  the  noise  spectrum  of  the  principal  beams.  The 
question  we  asked  was,  is  the  electrooptic  effect  in  Bii2Si02o  enhanced  for  a  grating 
vector  out  of  the  plane  of  the  resonator  at  35°?  The  answer,  as  demonstrated  in 
Appendix  B,  is  no. 

A  telling  clue  to  the  origin  of  the  companion  orders  is  seen  in  Fig  11.  An 
expanded,  incident  beam  with  no  voltage  applied  to  the  crystal  produced  a  trans¬ 
mitted  beam  as  well  as  a  second  beam  at  an  angle  of  35°.  It  seemed  our  crystal 
sample  had  a  slight  wedge  angle. 

4.3.3  Measurement  of  the  wedge  angle 

This  particular  crystal  of  Bii2Si02o  had  come  back  recently  from  being  recut  and 
polished.  It  was  part  of  a  crystal  of  original  dimensions  5mm  x  5mm  x  10mm.  After 
a  long  and  protracted  bout  with  an  internal  fracture,  the  first  crystal  split  in  two 
(see  Appendix  A  describing  damage  in  Bii2SiO20  crystals).  The  larger  piece  was 
saved  and  became  the  crystal  we  used  for  these  experiments.  It  is  very  likely  that 
during  the  cutting  or  polishing  the  wedge  was  introduced. 
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Figure  11:  Transmission  of  a  uniform  beam  through  the  crystal  sample.  The  pres¬ 
ence  of  a  multiple  reflection  suggests  a  wedge  angle. 


Figure  12:  Multiple  reflections  in  a  crystal  with  a  wedge  angle. 
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Figure  13:  Interferometric  measurement  of  wedge  in  BiizSiOjo  crystal. 

The  angular  separation  of  the  companion  orders  in  air.  A 9.  was  measured  as  the 
ratio  of  spatial  separation  over  a  given  distance  behind  the  crystal.  This  was  found 
to  be  0.19°  for  all  the  companion  orders  observed. 

As  shown  in  Fig  12,  the  wedge  angle,  9W ,  is  related  to  this  angular  spacing  by 

A 9  =  9.t  -90  =  2 n9w  (16) 

where  n  is  the  index  of  refraction.  For  our  measurement,  9W  =  0.036°.  Notice  that 
the  angular  separation  of  multiple  reflections  with  the  input  beam  is  a  multiple  of 
2 9W.  In  other  words,  the  higher  multiple  reflections  are  at  the  same  angular  spacing 
as  higher  diffracted  orders  produced  by  two- wave  mixing  of  the  input  beam  and  the 
first  co-directional  multiple  reflection. 

The  presence  of  a  wedge  in  the  crystal  was  confirmed  using  a  Specac  interferom¬ 
eter.  This  is  shown  in  Fig  13,  a  sketch  of  the  interference  pattern.  A  fringe  analysis 
was  difficult  because  both  the  wedge  angle  and  crystal  face  were  small.  The  orien¬ 
tation  of  the  fringes  were  at  approximately  35°  to  the  001  crystal  axis.  The  fringes 
were  bent,  suggesting  one  or  both  surfaces  were  curved.  This  may  explain  the  fo¬ 
cusing  effects  we  saw.  The  two  interference  fringes  is  too  small  a  number  to  give 
a  value  for  the  wedge  angle  of  0.036°,  but  reflection  off  the  curved  surfaces  might 
emulate  a  larger  wedge  angle. 
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4.3.4  Nature  of  the  companion  orders 

Are  the  companion  orders  merely  the  result  of  amplification  of  multiple  reflections? 
This  cannot  be  the  complete  answer  because  we  see  companion  orders  on  either  side 
of  the  principal  beam.  Multiple  reflections  exist  only  on  one  side.  In  the  same  wav 
the  interaction  of  the  pump  and  oscillating  beams  produce  HDOs,  the  two-wave 
mixing  between  any  principal  beam  and  its  first  codirectional  multiple  reflection 
gives  rise  to  higher  diffracted  “companion”  orders.  Additional  weak  multiple  re¬ 
flections  at  the  proper  angular  displacements  may  well  facilitate  the  emergence  of 
higher  diffracted  companion  orders.  Because  the  wedge  angle  was  smaller  than  the 
interbeam  angle  between  the  pump  beam  and  the  oscillating  beam  (roughly  ten 
times  smaller),  a  greater  number  of  companion  orders  were  present  than  HDOs. 
Moving  out  from  the  pump  beam,  the  intensity  of  the  principal  beams  decreased, 
and  thus  fewer  companion  orders  were  seen. 

Once  the  HDOs  and  companion  orders  are  established,  an  intricate  network  of 
beam  coupling  builds  up.  For  example,  each  separate  beam  couples  to  the  pump 
beam.  This  is  evident  because  when  the  oscillating  beam  is  suddenly  interrupted, 
the  HDOs  and  their  companion  orders  decrease  in  brightness,  then  gradually  fade 
away  completely.  Clearly  the  pump  beam  is  reading/erasing  a  number  of  gratings. 
A  complete  description,  in  which  each  beam  is  coupled  to  every  other  beam,  would 
be  complicated  indeed. 

4.4  Conclusions 

Our  goal  in  these  experiments  was  to  gain  expertise  with  a  resonator  using  Bi^SiOzo 
as  an  active  medium.  In  the  course  of  the  experiments  we  have  put  numbers  to 
many  qualitative  findings  previously  reported,  and  have  found  reasonable  agreement 
with  theory,  if  one  allows  for  variation  in  crystal  parameters  and  applied  fields. 
We  have  demonstrated  the  crucial  importance  of  the  optimum  detuning  frequency 
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to  oscillator  operation.  We  discovered  that  oscillation  in  a  URR  resonator  using 
Bi12Si02o  has  a  threshold  for  incident  intensity.  For  our  crystal  sample  no  oscillation 
occurred  for  intensities  less  than  0.8  mW  'em’.  And  we  have  observed  and  explained 
the  presence  of  various  higher  diffracted  orders,  some  of  which  are  due  to  a  slight 
wedge  angle  in  the  crystal  sample  we  used.  In  Chapter  6  we  will  discuss  potential 
applications  for  a  crystal  with  a  wedge  angle. 


Chapter  5 

Wave  Amplification  with  Two 
Pump  Beams 


So  we  grew  together, 

Like  to  a  double  cherry,  seeming  parted, 

Bui  yet  an  union  in  partition; 

Two  lovely  berrtes  moulded  on  one  stem. 

A  Midsummer  Night’s  Dream 
Act  III,  Scene  II 

Thus  far  we  have  looked  at  wave  mixing  in  photorefractive  materials  with  two 
input  beams.  A  characteristic  of  the  wave  interaction  is  power  transfer  in  one 
direction,  eg.  I+  increases  at  the  expense  of  but  not  vice  versa.  The  coupling 
term  which  describes  the  power  transfer  depends  on  the  imaginary  part  of  the 
modulated  space  charge  field. 

But  there  are  other  means  of  redistributing  energy  if  a  greater  number  of  waves 
interact  and  the  interbeam  angles  are  small.  A  recent  article  by  Ringhofer  and 
Solymar  suggested  a  new  mechanism  for  wave  amplification  in  this  case  [110].  In 
this  new  scheme,  four  waves  interact  as  shown  in  Fig  1.  Power  is  transferred  inwards 
from  the  pump  beams.  Unlike  our  traditional  two-wave  mixing  mechanism  which 
is  necessarily  on-Bragg  (the  beams  which  write  the  grating,  simultaneously  read  it 
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pr  xtal 


Figure  1:  Four  wave  mixing.  Beams  1  and  2  are  strong  pump  beams. 


directly),  this  new  mechanism  depends  on  an  indirect,  off-Bragg  coupling.  This  is 
similar  to  the  process  by  which  beams  can  be  amplified  via  a  higher  diffracted  order 
[6]. 


In  this  chapter  we  review  the  theory  of  Ref  [110]  (which  gives  only  the  final 
expression  for  the  gain  coefficient).  We  will  concentrate  on  qualitative  features  of 
the  theory.  One  important  criterion  is  that,  for  high  gain,  the  four  mixing  waves 
must  be  frequency  shifted  relative  to  one  another.  A  common  method  of  achieving 
this  small  frequency  shift  experimentally  is  to  introduce  a  Doppler  shift  by  reflection 
off  a  slowly  moving  mirror.  The  mirror  is  typically  mounted  on  a  piezoelectric 
translator  driven  by  a  sawtooth  voltage  signal  [109].  This  works  fine  for  one  beam, 
but  coordinating  the  motions  of  two  or  three  mirrors  to  get  multiple  frequency 
shifts  is  difficult  and  requires  more  equipment  than  we  had  available  to  us.  One 
way  out  is  to  use  a  ring  resonator.  The  oscillating  beam  in  the  ring  resonator  is 
automatically  frequency  shifted  to  give  optimum  gain.  In  this  arrangement,  one  of 
the  internal  beams  (3  or  4  in  Fig  1)  is  reinforced  in  a  resonator,  and  one  needs  only  a 
single  moving  mirror  to  detune  a  second  pump  beam.  The  remainder  of  the  chapter 
relates  experimental  work  with  this  resonator  and  a  similar  arrangement — a  linear 
resonator. 


Our  goals  in  the  experiments  were  to  verify  the  existence  of  the  new  gain  mech¬ 
anism,  to  test  certain  qualitative  features  of  the  theory,  and  to  establish  some 
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benchmark  by  which  the  new  mechanism  could  be  compared  to  the  traditional  gain 
mechanism. 

5.1  Theory 

Consider  the  four  waves  mixing  in  Fig  i.  The  interference  pattern  produced  by  the 
overlapping  of  these  waves  is  merely  a  superposition  of  those  patterns  produced  by 
each  pair  <~f  waves.  The  total  intensity  is 

1(f)  -  \E\2  =  i  £  f 

(=i 

=  +  l^m  (1) 

1=1  m=l 

where  =  ui|  -  are  the  detuning  frequencies,  and  Ktm  =  -  fcm  are  the 

resulting  grating  vectors.  We  assume  nearly  degenerate  mixing. 

We  can  follow  the  general  treatment  of  Chapter  2  and  derive  a  modulated  space 
charge  field  for  each  grating.  As  in  that  chapter,  we  neglect  higher  diffracted  or¬ 
ders  and  assume  small  modulation.  We  can  reduce  the  complexity  somewhat  by 
assuming  beams  1  and  2  are  very  much  stronger  than  beams  3  and  4.  Thus  we 
may  neglect  coupling  between  beams  3  and  4.  In  Chapter  2  we  found  that  power 
transfer  was  in  one  direction  only.  Let  us  take  this  direction  in  a  clockwise  sense 
in  reciprocal  space  (Fig  2a)  so  power  is  transferred  1  — ♦  3  — ♦  4  — ►  2.  We  let  our 
subscripts  reflect  this  direction.  The  modulated  dielectric  constant  is 

e,  =  er0  + erJr  + 

\  io  io 

4-  Euit  •)  +  Eu„  0  4. 

Io  Io 

+  E„„  1)  +  c  c_  )  (2) 

Io 

where  Iim  —  E*Em  and 


CHAPTER  5:  WAVE  AMPLIFICATION  WITH  TWO  PUMP  BEAMS 


no 


Figure  2:  Nearly  degenerate  four-wave  mixing,  (a)  Propagation  vectors  in  reciprocal 
space.  K  =  Jcx  -  fcj.  (b)  Off-Bragg  parameter,  V>- 

reflects  the  fact  that  the  interbeara  angles  are  small  and  we  are  in  the  drift  domi¬ 
nated  regime.  Here  b,m  -  6u>imTd,  £m(m  =  jrN^/fxKIm,  and  E,lm  =  eN^/t,Ktm. 


5.1.1  On-Bragg  and  off-Bragg  coupling 

The  writing  beams  read  the  various  gratings  formed  in  the  photorefractive  medium. 
Finding  the  coupling  between  the  beams  is  rather  like  doing  a  jigsaw  puzzle.  We 
need  to  find  some  combination  of  a  writing  beam  and  a  grating  vector  which  will 
exactly  reproduce  any  one  of  the  writing  beams,  ie.  duplicate  the  propagation  vector 
and  frequency.  If  the  pieces  fit,  we  have  vector  closure,  and  the  on-Bragg  condition 
of  Fig  2.6  is  satisfied.  For  example,  beam  1  combines  with  grating  vector  K 13  and 
part  of  it  is  diffracted  into  beam  3. 

Sometimes  the  jigsaw  pieces  do  not  meet  exactly,  and  we  need  a  hammer  to  force 
the  fit.  The  hammer  is  the  off-Bragg  vector,  rj>,  shown  in  Fig  2b.  Vector  closure  is 
given  by 

£3  =  fci  +  fcj  —  £4  +  t? 


(4) 
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Thus,  for  example,  beam  1  combines  with  grating  vector  A"42  o.nd  v  to  reproduce 
beam  3.  This  is  indirect,  or  off-Bragg,  coupling.  .Vote,  for  this  to  work  in  general, 
we  must  match  the  frequency  detunings.  Su)!3  —  8u > 42  —  C«5wi2  and  5u;!4  =  8ui 23  = 
(1  -  <l)6u )2,  where  (  is  the  scaling  parameter  shown  in  Fig  2a.  This  condition 
suggests  itself,  regardless,  because  of  the  symmetry  of  the  mixing.  This  also  ensures 
that  all  gratings  rmve  at  the  same  velocity. 

In  the  undepleted  pump  approximation,  where  /3,/4  <C  I\,I2,  the  coupled  dif¬ 
ferential  equations  for  E3  and  E4  become 

djE'x  o  a 

=  -jyo(Ev„I ,  +  E„„M2)  E3  -  jf(E^r  -  £«„)  ExE2E\t~^  (5) 

dEt  a  a 

=  - jfo(E E<  -  j±  (Eult-  4-  Eu„)  E,E2E‘3e~^  (6) 

where  q  =  ~k0n3r.  The  first  term  on  the  right  hand  side  of  the  above  equations 
represents  on-Bragg  coupling,  the  second  term  off-Bragg. 


5.1.2  Analytic  treatment  of  the  coupled  equations 

To  make  Eqns  5  and  6  more  tenable,  we  introduce  a  change  of  variables 
B3  =  E3  exp  (jy),  =  El  exp  ( ) 


Then  our  coupled  equations  become 


dB3 

■  ( 

—  F  -L 

dx 

}q\i+P 

ro 

-jq{E„„' 

+  Eu  „ ) 

dB< 

.  (  Pro 

F  *  _L_ 

dx 

0 

+ 

r-H 

+  iq{EWii' 

+  ) 

-E'  -  B 


dx  -\l+(3ro  l+/3ro 

*  '  Pro 

where  0r<>  =  h/h-  These  equations  are  of  the  form 


1  +  fir 


l+Pro 

rri-E...  -  T)  b,  + 


B3  —  —j{a\\B3  -f  &12B4) 
B\  —  —j{anB3  -f  ojjB4) 


(7) 

(8) 

(9) 

(10) 

(11) 
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where 
an  —  9 


r _ 1 _  p  . 

1  +  3ro  ""  1  -r  Jro 


aiz  —  '/ 


V  3ro  /  r-  . 

a21  -  -9"  7~  (£*33 

1  -r  Pro 


V  3ro 

1  -  3re 

3 , 


(£-„‘  -  £ 


'•*'4  2  ) 


0,22  —  —'l 


1  -  Jr 


.  F  •  _ 


_ = — E 

1  ?  w° 
i  Jro 


The  solutions  are 

E3{x)  =  B3  exp(-j^-)  =  £3(0)  exp(-j(T -r  |)x)  (12) 

E4(x)  =  54  exp  (j )  =  £4(0)  exp(j(T: -  ^  )*)  (13) 

where 

T  =  2 (^n  —  a22)  ~  ;iy  (an  —  a22)2  -  4aj2a2i  (14) 


5.1.3  The  gain  coefficient 

The  coupling  term  T  is  complex.  Let  T  =  Tr  +  jT,.  Then 

£»(*)  =  E,(0)exp(Ti*-j(Tr  +  ^)*)  (15) 

E4(x)  =  £l4(0)exp(Tix+i(TF-lV')*)  (16) 

The  intensities,  I3  and  /4,  will  exhibit  simple  exponential  gain,  expressed  by  a  gain 
coefficient,  2T<,  which  is  the  new  mechanism  analogue  of  2T  sin$9.  An  interesting 
feature  is  the  phase.  If  \ip  <C  Tr,  then  beams  3  and  4  are  phase  conjugates. 

Let  us  look  at  T  more  closely.  We  can  simplify  Eqn  14  somewhat  by  noting 
=  Ewa,  and  E„it  =  Eu „.  So  ai2  =  — aji-  Our  general  expression  for  T  is 


*  ?’{i  +J${E-  ~  "  2^E-»  +  E-"> + *' 

(17) 

So  our  gain  term  for  intensity  is 

2T<  =  Re 

2 q*(E„n  +  Eux;)  -i>>-  j^(E„„  +  £„„*)’ 

(18) 
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The  first  term  on  the  right  hand  side  of  the  equation  is  due  to  on- Bragg  coupling 
between  the  internal  beam  and  the  two  pump  beams.  Power  is  transferred  from 
pump  2  to  the  internal  beam,  but  lost  to  pump  1.  The  second  term  is  due  to 
off-Bragg  coupling. 

It  is  worthwhile  to  look  at  two  special  cases. 

I.  Three-wave  mixing.  (J  =  ~. 

Here  both  internal  waves  are  coincident,  so  we  have  in  effect  three- wave  mixing.  In 
this  case,  E ’WM  =  Eqn  18  becomes 


2T,  -  -2gJ— ^  Im(  £«>.)+  \  4^  Ref  EUJ  -  -  4q'-  (R  e(£„14))2 

1  4-  dro  \  (  *■  f  dTO ) 

(19) 

Here  it  is  obvious  that  on-Bragg  coupling  depends  on  the  imaginary  part  of  the 
modulated  space  charge  field,  and  the  off-Bragg  coupling  on  the  real  part.  Two  other 
features  are  important  here.  The  sign  of  g,  which  depends  on  rejf,  can  be  positive 
or  negative  by  virtue  of  crystal  orientation.  This  is  a  consideration  in  optimizing 
2Ti.  For  example,  to  get  large  gain  when  q  is  positive,  we  must  minimise  Im(£Ull) 
and  balance  the  two  terms  containing  Ref EUii)  simultaneously.  The  second  feature 
is  the  pump  ratio,  >3ro.  For  positive  q  and  Im( Euit  ),  a  beam  ratio  greater  than  one 
is  necessary  for  effective  on-Bragg  coupling.  If  the  beam  ratio  is  less  than  one,  then 
this  term  detracts  from  the  overall  gain. 

II.  Equal  pump  intensities.  (3ro  —  1. 

Here 


2T<  =  Re  ^2g^(E„„  +£„,/) -tf1) 


(20) 


This  is  the  expression  from  Ref  [110).  There  are  two  strategies,  mathematically,  for 
maximising  the  red  part  of  this  square  root  (we  assume  C,  ^  \).  Which  strategy  we 
choose  depends  on  the  the  size  of  the  off-Bragg  parameter,  ifi.  It  has  to  be  different 
from  zero  to  have  gdn,  but  if  it  is  too  large,  the  rp7  term  dominates  and  gain  is 
reduced  or  lost  completely.  In  our  jigsaw  puzzle  andogy,  when  the  hammer  is  too 
large  it  smashes  the  pieces  rather  than  making  them  fit.  Therefore  we  expect  the 
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gain  to  he  effective  for  small  interbeam  angles.  However,  for  very  large  values  of 
V  one  can  achieve  high  gain  if  the  real  part  of  ( Eu^  —  * )  is  minimised  while 

its  imaginary  part  is  maximised  (and  has  a  positive  sign).  For  this  to  be  effective, 
though,  (  must  be  different  from  1  In  the  case  of  three- wave  mixing  then 

t 

2Ti  =  y Aqtt>  -  v2  (21) 

The  best  test  for  the  presence  of  the  new  gain  mechanism  is  to  ha/e  equal  pump 
intensities.  There  is  no  on-Bragg  coupling  in  this  case. 

5.1.4  Optimization  of  the  gain 

Let  us  look  at  the  problem  of  achieving  high  gain  for  a  particular  case.  We  will 
restrict  ourselves  to  the  combination  of  three- wave  mixing  and  equal  pump  inten¬ 
sities.  The  gain  is  optimized  for  the  maximum  value  of  q  R e(  For  the  case 

Eqit  »  E. o  we  have 

E.o  fl  +  bn 

Re(E^)  =  - (22) 

(1  +  6l<zC)  +b“ 

which  has  extrema  for 


Fig  3  shows  2T,  as  a  function  of  pump  detuning  frequency  (bvu)  for  various 
interpump  angles.  Fig  3a  is  for  positive  q,  and  Fig  3b  is  for  negative  q.  For  all  the 
theoretical  plots  in  this  chapter  the  crystal  parameter  values  are  those  which  give 
a  good  qualitative  fit  to  the  experimental  data — that  is,  providing  the  correct  peak 
frequency  detuning  for  the  symmetric  case  of  three-wave  mixing.  The  differences 
between  these  and  those  in  Table  2.1  are  —  2  x  10JJ  m~3,  and  the  ratio  fj,/yr  = 
1.9  x  10,J  (V-m)-1.  We  use  an  applied  field  of  8  kV/cm  and  a  total  pump  intensity 
of  15  mW/cm2. 

It  is  clear  there  is  one  best  overall  interpump  angle.  For  these  parameters,  it  is 
approximately  3.2°  for  both  positive  and  negative  q.  This  angle  gave  us  the  best 
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experimental  gain.  An  important  point  is  the  magnitude  of  the  gain  coefficient.  At 
optimum  detuning  and  grating  spacing,  the  gain  is  37  cm  ‘.  In  the  case  of  nearly 
degenerate  two-wave  mixing,  we  found  at  small  angles.  ■E'u.opti  —  E If  we  compare 
mixing  between  pump  1  and  the  inside  beam,  then  the  same  for  the  same  crystal 
parameters  the  largest  expected  gain  coefficient  is  35.5  cm'1.  (Note:  These  values 
are  theoretical  values.  In  practice,  because  of  loss,  this  gain  coefficient  is  much 
less- — a  typical  optimum  value  in  Bi^SiOjo  is  12  cm'1.)  This  new  gain  mechanism 
seems  to  be  slightly  more  efficient  than  the  traditional  mechanism. 

5.1.5  Variation  in  C 

If  £  1  then  we  have  four- wave  mixing.  Fig  4  shows  the  gain  coefficient  as  a  function 

of  £  for  various  pump  detuning  frequencies  for  two  extremes  of  in-air  angles  between 
the  pump  beams.  Fig  4a  is  for  an  interpump  angle  of  1.6°  and  positive  q,  and  Fig  4b 
for  negative  q.  For  a  given  value  of  £  there  is  an  optimum  detuning  frequency. 
Figs  4c  and  4d  show  an  interpump  angle  of  10°.  This  gain  is  only  possible  when 
the  imaginary  part  of  (EUiJ  +  /£«,/)  is  maximised  and  (  much  different  than 
The  symmetry  and  interrelatedness  of  the  internal  beams  is  very  important  in  these 
plots.  Both  beams  experience  the  same  gain. 

At  this  point  we  should  mention  the  appearance  of  spatial  subharmonics  be¬ 
tween  mixing  beams  in  BijjSiC^o-  This  has  been  reported  by  Mallick,  et  al.  [89). 
The  authors  observed  K/2,K/2 ,  and  K/ 4  subharmonics  (where  K  is  the  grating 
vector  formed  by  two  nearly-degenerate  mixing  waves).  These  are  very  distinct 
beams  which  may  attain  an  intensity  relative  to  the  pump  beams  of  40%.  Which 
submultiples  were  visible  depended  on  the  amount  of  frequency  detuning.  Ref  [89] 
suggested  these  beams  were  diffracted  orders  which  overlapped  and  were  mutually 
reinforced. 

The  spatial  subharmonics  reported  correspond  to  £  =  1/2, 1/3,  and  1/4  for  our 
four- wave  mixing.  Fig  4  shows  that  high  gain  is  possible  for  these,  and  other,  values 
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Figure  5:  Gain  coefficient,  2T^  vs  pump  beam  ratio  for  positive  and  negative  q. 

of  (.  Perhaps  these  intrapump  beams  arise  from  noise  amplified  by  a  combination 
of  the  gain  mechanism  described  in  this  chapter,  and  a  reinforcing  mechanism. 
The  phenomenon  of  subharmonics  does  not  seem  to  be  a  general  characteristic  of 
BiijSi02o — we  have  looked  for  it,  but  have  never  observed  distinct  subharmonics  in 
any  of  our  crystal  samples. 

5.1.6  Variation  in  pump  beam  ratio 

The  relative  strength  of  the  pump  beams  is  important  to  the  gain  term.  If  the 
ratio  is  much  different  than  1,  the  on-Bragg  coupling  becomes  important  and  may 
dominate  the  gain.  In  any  case,  pumps  of  unequal  intensity  produce  a  term  inside 
the  square  root  of  Eqn  19  which  detracts  from  the  off-Bragg  coupling. 

The  effect  of  changing  0TO  on  the  gain  coefficient  is  seen  in  Fig  5.  Here  (  = 
the  interpump  angle  is  3.2°  and  for  positive  g,  hv\i  is  taken  as  -10.5  Hz,  while  for 
negative  q,  6 i/\j  is  -12.5  Hz.  These  frequencies  maximise  the  gain  coefficient  for 
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Figure  6:  Unidirectional  ring  resonator  with  two  pump  beams. 

5.2  Unidirectional  Ring  Resonator:  Experimental 
Results 

5.2.1  Experimental  layout 

We  set  about  to  prove  experimentally  the  existence  of  the  new  gain  mechanism.  The 
obvious  way  of  doing  the  experiment  is  to  have  3  or  4  input  beams  and  measure  the 
power  transfer  into  the  signal  beams.  But  the  difficulty  is  then  to  adjust  carefully 
the  frequency  detuning  of  each  beam  needed  to  maximise  the  interaction.  We  opted 
for  a  different  experimental  set-up  (Fig  6)  in  which  there  are  only  two  input  beams 
(the  two  pump  beams)  and  then  one  of  the  signal  beams  may  emerge  from  the 
noise.  In  order  to  help  the  emergence  of  the  signal  beam  we  created  a  feedback  loop 
by  constructing  a  three-mirror  resonator.  This  is  not  an  ideal  test  of  the  theory 
because  the  interaction  will  now  depend  on  additional  phenomena  like  the  scattering 
of  the  input  pump  beams  and  on  new  parameters  like  the  resonant  frequency  of  the 
cavity,  and  in  our  undepleted  pump  approximation  there  would  be  no  limit  to  the 
exponential  growth  of  the  signal  beam.  However,  this  beam  can  only  emerge  in  the 
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resonator  if  there  is  an  amplification  mechanism.  Therefore  the  present  experiments 
can  prove  the  existence  of  the  mechanism,  but  they  are  not  suitable  for  checking 
quantitatively  the  predictions  of  the  theory. 

The  photorefractive  material  used  was  the  Bi12Si02o  sample  described  in  Chap¬ 
ter  4.  For  all  runs  a  4  kV  d.c.  voltage  was  applied  along  the  001  crystal  axis. 
Vertically  polarized  light  of  wavelength  514  nm  was  incident  on  the  110  face.  One 
pump  beam  (which  we  shall  call  driving  pump)  was  Doppler  shifted  using  a  mirror 
mounted  on  a  piezoelectric  translator.  The  input  signal  to  the  translator  was  a 
triangular  wave.  During  half  the  period  the  frequency  shift  was  favourable  for  os¬ 
cillation  and  amplification,  and  during  the  other  half  conditions  were  unfavourable. 
The  ring  resonator  cavity  was  aligned  using  the  main  pump  beam.  An  iris  aperture 
in  one  of  the  cavity  legs  allowed  only  low  order  modes  to  oscillate.  The  oscillating 
beam  was  sampled  from  reflection  off  a  glass  plate. 

5.2.2  Verification  of  new  gain  mechanism. 

First  the  set-up  was  aligned  so  that  the  angle  between  the  main  pump  beam  and 
the  oscillating  beam,  9U,  was  the  same  as  the  angle  between  the  driving  pump 
beam  and  the  oscillating  beam,  do-  Then  we  switched  on  the  main  beam  on  its  own 
and,  as  expected,  obtained  oscillation  due  to  the  traditional  gain  mechanism,  that 
is  the  oscillating  beam  adjusted  its  own  frequency  detuning  so  as  to  maximise  the 
imaginary  part  of  the  space  charge  field.  Next  we  switched  on  the  driving  beam  of 
the  same  intensity  as  the  main  beam.  The  oscillation  stopped,  as  predicted,  because 
the  two  pump  beams  worked  against  each  other.  The  oscillating  beam,  however, 
reappeared  if  we  slightly  detuned  the  driving  beam.  Such  detuning  would  not  have 
helped  the  traditional  gain  mechanism  (recall  there  is  no  gain  via  the  imaginary  part 
of  the  space  charge  field  when  there  are  two  pump  beams  of  equal  intensity)  but 
could  maximise  the  new  gain  mechanism.  Hence  the  reappearance  of  the  oscillating 
beam  is  a  proof  for  the  existence  of  the  new  gain  mechanism. 
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Figure  7:  Gain  vs  ratio  of  pump  beam  intensities.  Gain  is  defined  as  the  photodiode 
signal  with  both  pumps  over  the  photodiode  signal  with  only  the  main  pump. 

5.2.3  Comparison  with  traditional  gain  mechanism 

The  relative  intensity  of  the  oscillating  beam  as  a  function  of  pump  beam  ratio 
is  plotted  in  Fig.  7.  The  interpump  angle  in  air  was  3.2°  with  the  oscillating 
beam  bisecting  this  angle.  The  intensity  ratio  {Imam/ Idrive)  was  varied  using  a 
neutral  density  filter  wheel.  The  photodiode  signal  was  taken  at  optimum  driving 
pump  detunirg  frequency  and  this  value  was  normalised  to  the  average  signal  for 
oscillation  with  no  driving  pump.  The  driving  pump  had  no  noticeable  effect  for 
ratios  >  50:1 — the  oscillating  signed  was  indistinguishable  from  that  due  to  the 
main  pump  alone.  For  ratios  <  1:10  oscillation  was  surpressed.  Presumably  for 
very  weak  main  pumps  the  new  gain  mechanism  is  not  efficient — the  direction  of 
energy  exchange  from  the  driving  pump  for  the  traditional  gain  mechanism  is  wrong 
for  our  geometry.  The  optimum  beam  ratio  was  unity,  giving  an  increase  in  signal 
intensity  by  about  a  factor  of  four.  This  overall  behaviour  agrees  qualitatively  with 
Fig  4a  (for  the  case  of  negative  q ),  although  the  peak  there  did  not  occur  for  0Te  =  1. 
Our  crystal  was  set  up  so  the  energy  transfer  was  2  — »  4  — »  3  — »  1,  corresponding 
to  negative  q.  We  also  see  that  the  new  gain  mechanism  is  more  efficient  than  the 
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traditional  one.  In  a  system  '.vith  feedback,  such  as  this  resonator,  a  small  increase 
in  2T,  can  be  magnified  into  a  large  increase  in  gain. 

Each  data  point  represents  the  average  of  the  photodiode  signal  for  10-30  cycles 
of  the  triangular  input  to  the  piezoelectric  translator.  The  error  bars  represent  one 
standard  deviation  for  the  data.  The  likely  reason  for  the  noisiness  of  the  data  is 
the  rather  short  interaction  length  of  the  crystal,  only  5.3  mm,  thus  giving  only  a 
limited  amount  of  amplification.  The  overall  plot  and  trends,  however,  were  easily 
repeatable. 


5.2.4  Optimization  of  gain 

The  dependence  of  the  oscillating  signal  on  the  driving  pump  beam  detuning  is 
shown  for  a  pumpbeam  ratio  of  1:2  (driving  pump  intensity  =  10  mW/cmz).  Fig  8a 
shows  the  symmetric  case  of  0M  =:  Op  =  1.6°,  with  error  bars.  The  next  two  plots 
(without  error  bars)  show  the  effect  of  keeping  9m  =  1.6°  and  varying  Op,  and  for  9p 
=  1.6°  and  varying  0m-  The  trends  are  quite  clear.  In  all  cases  there  is  an  optimum 
detuning  and  the  signal  declines  as  Op  or  0m  increases. 

5.2.5  Search  for  subharmonics 

Let  us  discuss  some  further  experiments  related  to  those  of  Mallick  et  al.  who  also 
had  two  pump  beams  incident  on  a  Bi^SiOjo  crystal  and  observed  spontaneously 
emerging  output  beams  at  subharmonic  positions.  We  asked  the  following  question: 
Is  the  signal  beam  in  the  oscillator  higher  at  a  subharmonic  position  that  in  its 
vicinity,  or  in  other  words,  are  subharmonic  positions  in  any  way  distinguished? 

To  answer  this  question  we  set  0M  =  1-6°  and  Op  =  3.2°,  that  is  when  the 
oscillating  beam  trisects  the  inter-pump  angle  and  we  chose  unity  beam  ratio.  This 
corresponds  to  the  position  of  the  third  subharmonic.  By  varying  either  0m  or  &D 
or  both  we  did  not  observe  any  preference  for  the  subharmonic  position.  This  seems 
to  indicate  that,  as  predicted  by  the  theory  of  Ref  [110],  the  subharmonic  positions 
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are  by  no  means  preferred  and  gain  is  continuously  present  for  the  range  of  s 
and  Op's  measured,  not  only  in  discrete  positions. 

There  is  no  doubt  that  subharmonics  were  observed  repeatedly  and  reproduciblv 
by  Mallick  et  al.  Why  is  it  then  that  neither  the  theory  of  Ref  ;  1 1 0 j  nor  our  exper¬ 
iments  can  find  any  preference  for  subharmonics?  It  is  possible  the  theory  fails  to 
find  the  subharmonics  because  of  the  approximations  used,  namely  small  modula¬ 
tion  and  no  higher  diffracted  orders.  Our  experiments  might  miss  the  subharmonics 
because  for  some  reason  it  is  the  presence  of  the  resonator  that  ensures  the  con¬ 
tinuously  available  gain.  Another  possibility  is  that  the  subharmonics  found  by 
Mallick,  et  al.  are  due  to  some  properties  of  their  experimental  set-up,  e.g.  those 
of  the  piezoelectric  drive. 

We  should  mention  an  effect  observed  when  in  the  range  of  pump  frequency 
detuning  correspoding  to  high  gain.  If  the  resonator  was  not  aligned,  or  a  screen 
was  placed  in  one  of  the  cavity  legs,  no  oscillation  occurred.  However,  a  large  patch 
of  optical  noise  appeared  in  between  the  pump  beams  when  the  driving  mirror 
movement  direction  was  favourable  for  gain.  This  region  could  be  fairly  intense, 
roughly  20  times  less  than  that  for  the  oscillating  beam.  There  were  even  “higher 
diffracted  orders’’  produced —  less  intense  regions  of  noise  on  either  side  of  the 
pumps.  It  is  clear  from  Fig  4a  that  gain  is  high,  for  small  angles,  for  a  range  of 
around  j.  Noise  in  this  region  undergoes  high  amplification. 

5.2.6  The  second  signal  beam 

The  appearance  of  the  second  signed  beam  is  an  important  verification  of  the  theory. 
The  new  gain  mechanism  depends  on  four- wave  mixing.  We  cannot  have  amplifica¬ 
tion  of  one  signal  beam  in  the  absence  of  the  other.  There  is  a  crucial  symmetry  in 
these  beams.  In  one  experiment,  we  varied  the  interpump  angle  from  3.4°  to  5.8°, 
keeping  the  angle  9M  constant  as  1.7°.  The  angle  between  the  second  signal  beam 
and  the  driving  pump  was  measured  at  1.7°  in  every  instance. 
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Figure  9:  Simultaneous  oscilloscope  traces  of  signal  strength  of  beams  3  and  4. 
Bottom  trace  is  oscillating  beam. 

Both  signal  beams  could  be  sampled  using  a  photodioae.  Fig  9  shows  a  typical 
simultaneous  trace  for  the  emergence  of  the  two  beams.  The  second  signal  beam 
was  weaker  (by  a  factor  of  5  to  10),  but  the  close  relationship  between  the  beams  is 
unmistakable.  If  the  oscillating  beam  was  suddenly  interrupted,  both  oscilloscope 
traces  immediately  dropped  to  background  level. 


5.3  Linear  Resonator  with  Two  Pump  Beams 

A  significant  finding  of  the  last  section  was  that  for  the  ring  resonator  the  gain  in 
the  system  with  two  pumps  was  larger  that  using  a  single  pump  (i.e.  exploiting  the 
traditional  photorefractive  gain  mechanism).  This  gave  us  encouragement  to  try  the 
two  pumps  with  a  linear  resonator  arrangement  with  Biia£  'O2o-  Using  just  a  single 
pump  we  have  been  unable  to  obtain  oscillation  in  Bi^SiOzo  with  this  geometry. 

The  linear  resonator  using  photorefractive  crystals  as  the  gain  medium  has  been 
investigated  theoretically  [24,65]  and  experimentally  using  BaTi03[131].  In  these 
studies,  a  cavity  is  formed  by  two  mirrors  (see  Fig  10),  and  a  single  pump  is  sufficient 
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Figure  10:  Linear  Resonator:  Experimental  layout. 

for  oscillation  to  begin.  The  oscillating  beams  arise  from  the  noise  spectrum  of  the 
input  pump  beam.  Those  components  which  satisfy  cavity  resonance  conditions 
and  which  have  the  appropriate  frequency  detuning  can  be  amplified.  If  these 
counterpropagating  beams  are  sufficiently  intense,  a  phase  conjugate  of  the  input 
pump  beam  is  produced. 

That  oscillation  in  the  linear  resonator  is  more  difficult  to  start  and  sustain  is 
readily  seen  in  Fig  10.  Amplification  is  effective  only  when  the  oscillating  beam 
is  codirectional  with  the  pump  beams.  The  oscillating  beam  must  traverse  the 
crystal  twice  to  get  back  to  a  direction  where  gain  is  possible.  We  have  twice  the 
absorption  loss.  The  situation  is  reminiscent  of  the  ferrymen  who  hauled  supplies 
and  convicts  out  to  the  prison  on  Alcatraz  island  in  San  Francisco  Bay — the  trip  was 
only  profitable  on  the  way  out.  This  configuration  using  Bii2SiO20  was  described 
by  Rajbenbach  and  Huignard  [105],  but  it  is  unclear  what  success,  if  any,  they  had 
in  achieving  oscillation. 

Everything  about  the  layout  was  the  same  as  for  the  ring  resonator,  except  a 
d.c.  voltage  of  4.75  kV  was  used.  Both  pumps  had  intensities  of  4  mW/cm2.  The 
k  interpump  angle  was  3.2°  in  air  and  the  oscillating  beam  bisected  this  angle.  As 

before,  the  driving  pump  was  detuned  using  a  mirror  mounted  on  a  piezoelectric 
translator.  We  used  a  triangular  ramp  signal  to  drive  the  mirror.  This  effectively 
changes  the  sign  of  the  frequency  detuning  halfway  through  each  cycle  and  for  one 
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Figure  11:  (a)  Typical  data  for  Linear  Resonator.  Top  trace  is  ramp  signal  to 
piezoelectric  mirror,  (b)  Oscillating  beam  intensity  vs  pump  detuning  frequency. 


CHAPTER  5:  WAVE  AMPLIFICATION  WITH  TWO  PUMP  BEAMS 


131 


direction  of  mirror  motion  we  expect  favourable  gain,  but  none  for  the  opposite 
direction.  This  effect  is  readily  seen  in  Fig  11a  which  depicts  some  sample  data. 
Here  6uj  =  87  s_1.  There  is  a  slight  delay  as  the  ramp  voltage  increases  before 
oscillation  begins.  This  might  be  explained,  in  part,  by  an  initial  slow  amplification 
from  the  noise  spectrum.  More  likely  it  is  due  to  some  chirp  in  the  detuning  of  the 
driving  pump  beam  due  to  inertia  as  the  mirror  changes  direction  of  motion. 

It  was  essential  for  both  pump  beams  to  be  present  to  maintain  oscillation.  There 
was  no  oscillation  for  either  pump  alone.  Thus  it  appears  that,  as  for  the  case  of 
the  ring  resonator,  this  new  gain  mechanism  is  more  efficient  than  the  traditional 
one. 

The  oscillating  beam  was  sampled  using  the  reflection  (20%)  off  the  front  face 
of  the  crystal.  Fig  lib  shows  the  variation  of  the  oscillating  beam  intensity  as  a 
function  of  the  frequency  detuning  of  the  driving  pump  beam.  There  is  clearly  an 
optimum  value.  Unfortunately  the  intensity  of  the  oscillating  beam  was  not  high 
enough  to  observe  phase  conjugate  reflection. 

5.4  Conclusions 

In  this  chapter  we  have  investigated  a  new  mechanism  for  wave  amplification  in 
photorefractive  materials.  This  mechanism  uses  two  pump  beams  and  is  very  effec¬ 
tive  in  Bij2Si02o  at  small  interbeam  angles  where  the  off-Bragg  parameter  is  small. 
Therefore  it  seems  to  be  a  useful  complement  to  nearly-degenerate  two- wave  mixing 
described  in  Ref  [109].  We  have  used  this  mechanism  successfully  in  two  different 
resonator  configurations. 


Chapter  6 

Applications  of  Wave- Mixing  in  a 
Crystal  Sample  with  a  Wedge 
Angle 

There  ij  nothing  insignificant. 

Samuel  Taylor  Coleridge 

The  photorefractive  effect,  when  first  reported,  was  nearly  dismissed  as  a  mere 
nuisance.  The  opening  paragraph  of  the  paper  by  Ashkin,  et  al.  [5]  addressed  the 
optically-induced  refractive  index  inhomogeneity  in  LiNb03  and  LiTa03: 

The  effect,  although  interesting  in  its  own  right,  is  highly  detrimental 
to  the  optics  of  nonlinear  devices  based  on  these  crystals. 

And  the  closing  remarks  proclaim  (italics  mine): 

Work  on  this  mechanism  and  its  possible  elimination  is  continuing. 

Chen  recognized  the  potential  of  the  change  in  refractive  index  as  a  means  of  storing 
volume  holograms  [20]  and  the  rest,  as  described  in  Chapter  1,  is  history. 
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The  moral  here,  of  course,  is  that  in  the  construction  of  technology,  sometimes 
stumbling  blocks  make  the  best  building  material. 

In  our  experiments  with  the  unidirectional  ring  resonator  described  in  Chapter 
4,  we  discovered  a  crystal  sample  of  Bi12Si02o  with  a  wedge  angle.  The  multiple 
reflections  created  by  the  wedge  and  their  resulting  higher  diffracted  orders  intro¬ 
duced  new  channels  through  which  incident  power  was  transferred.  The  result  was 
a  loss  in  efficiency  of  the  resonator.  It  would  be  simple  to  denounce  the  wedge  angle 
as  an  unfortunate  mishap  and  have  the  sample  repolished.  But  before  taking  that 
step  we  should  ask:  Is  there  any  use  for  a  crystal  with  a  wedge  angle? 

One  application  is  immediately  obvious  and  was  mentioned  in  Chapter  4.  One 
can  produce  many  output  beams  with  a  single  input  beam.  Since  the  wedge  angle 
is  fixed,  the  results  are  highly  reproducible— something  not  always  guaranteed  with 
photorefractive  crystals. 

This  chapter  addresses  other  ways  of  utilizing  our  “unfortunate  mishap”.  First 
we  look  at  some  interesting  transient  phenomena  when  the  applied  voltage  is  first 
turned  on.  This  provides  some  insight  into  how  power  is  transferred  among  the 
higher  diffracted  orders.  Finally  we  use  a  He-Ne  633  nm  beam  to  probe  the  grating 
formed  by  argon  ion  laser  light  of  514  nm.  The  probe  beam  was  coincident  with  the 
writing  beam.  We  observed  a  field/intensity  dependent  scattering  that  has  potential 
application  in  laser  beam  steering. 

6.1  Transient  Effects  in  Higher  Diffracted  Orders 

6.1.1  Transient  energy  transfer  between  mixing  beams 

There  is  a  certain  inertia  during  grating  formation  in  photorefractive  materials.  It 
takes  time  for  the  space  charge  field  to  build  up  to  modulate  the  refractive  index. 
The  build-up  period  is  characterised  by  transient  effects  and  these  can  be  quite 
spectacular.  For  example,  in  drift-dominated  two-wave  mixing,  the  power  transfer 
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is  oscillatory  and  during  the  undulations  of  intensity,  a  weak  signal  beam  can  nearly 
deplete  a  strong  pump  beam  ,45l. 

The  time  dependence  of  the  modulated  space  charge  field  produced  by  two 
mixing  waves  is  described  by  the  well-known  differential  equation  ( 1 1 6 i 

BE, ,  „  ,  m 

—Qt —  gE.  i  =  fvn  (1) 


where  m  is  the  modulation  ratio  of  the  the  mixing  waves  and 


+ 


Eqn  1  reduces  (as  it  must)  to  Eqn  2.19  in  the  steady-state.  Because  the  coefficients 


are  complex  in  this  first  order  differential  equation,  the  time  behaviour  of  Eti  can 
be  oscillatory. 

In  short,  transient  effects  are  described  by  Eqn  1  and  the  coupled  wave  equa¬ 
tions.  John  M.  Heaton’s  D.Phil.  thesis  provides  many  theoretical  and  experimental 
details  of  two-wave  mixing  transient  behaviour  in  Bii2Si02o  [49].  In  particular  he 
developed  both  numerical  and  analytic  solutions  to  the  coupled  time  and  space 
differential  equations.  Both  treatments,  though,  are  computationally  complicated. 
For  example,  the  analytic  solution  works  only  for  the  undepleted  pump  approxima¬ 
tion  and  requires  complex  integration  and  the  evaluation  of  Bessel  functions  with 
complex  arguments.  Appendix  C  reveals  a  new,  simple,  analytic  formula  we  have 
developed  for  the  signal  beam  complex  amplitude  which  agrees  remarkably  well  to 
the  numerical  solution  of  Heaton  and  Solymar  [48]  for  a  wide  range  of  important 
and  practical  cases.  The  savings  in  computer  time  is  marked.  A  typical  calcula¬ 
tion  of  220  time  points  using  either  method  of  Heaton  takes  20  minutes  on  an  IBM 
personal  computer,  while  the  computation  of  the  same  number  of  points  using  the 
simple  formula  requires  less  than  5  seconds. 

Transient  effects  have  been  studied  experimentally  in  two-wave  mixing  in  LiNb03 
[75],  Bii2GeO20  [45],  and  Bi12SiO20  [49,64].  Observed  effects  are 
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•  Power  is  transferred  from  the  strong  oeam  to  the  weak  beam. 

•  For  small  interbeam  angles  and  small  applied  fields  the  transfer  is  oscillatory. 

•  For  strong  applied  fields  the  gtating  writing  times  are  longer. 

•  For  smaller  grating  spacings  (larger  interbeam  angles)  the  grating  writing 
times  are  longer. 

Our  study  is  the  first  experimental  look  at  transient  power  transfer  among  higher 
diffracted  orders. 


6.1.2  Videotape  studies 

The  business  of  applying  a  voltage  to  a  crystal  may  seem  rather  mundane,  but  there 
are  definite  considerations.  Bi12SiO20  is  a  very  good  insulator  in  the  dark  (resistivity 
=  5  x  1013fi  cm  [2]).  One  must  be  careful  of  having  a  high  voltage  across  an 
unilluminated  crystal.  The  danger  is  an  in-air  spark  between  the  electrodes  which 
might  damage  the  crystal  (see  Appendix  A).  The  crystal  must  be  flooded  with  light 
to  reduce  the  resistivity  several  orders  of  magnitude  before  the  voltage  is  applied. 

The  turn-on  characteristics  of  our  HV  power  supplies  were  not  step-functions. 
There  was  no  instantaneous  jump  from  0  V  to  the  desired  potential  difference. 
Instead  the  voltage  (as  measured  from  the  current  across  a  fixed  10  Gfl  resistor — a 
value  close  to  the  resistance  of  our  illuminated  crystal)  had  rise  time  constants  on 
the  order  of  a  few  milliseconds,  and  generally  experienced  some  type  of  overshoot 
above  the  final  steady-state  voltage.  For  example,  if  the  power  supply  was  turned 
to  1  kV,  the  voltage  peaked  at  2.7  kV  (an  overshoot  of  1.7  kV)  before  dropping  back 
to  1  kV.  For  a  desired  voltage  of  4.0  kV  the  rise  time  was  longer  and  the  overshoot 
was  1.0  kV  (underdamped  behaviour). 

Modelling  a  time-varying  applied  field  like  this  is  difficult,  and,  for  the  most 
part,  unnecessary.  Studies  of  transient  behaviour  in  photorefractive  crystals  employ 
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choppers  or  shutters  for  one  of  the  mixing  beams  with  the  applied  voltage  at  its 
steady-state  value  i 49] . 

But  the  simple  fact  remains:  at  some  instant,  one  has  to  turn  on  the  voltage 
with  the  crystal  illuminated.  Not  so  simple  is  that  when  this  occurs,  a  pulse  of  light 
energy  escapes  from  the  crystal,  away  from  the  region  of  transmitted  light. 

This  effect  is  observed  to  some  extent,  in  all  of  our  Bo^SiC^o  crystal  samples. 
The  pulse  duration  is  short — less  than  one  second  for  intensities  on  the  order  of 
1  mW/cm2.  This  is  too  fast  for  ordinary  camera  film  to  pick  up,  but  not  for  a 
TV  camera.  We  observed  scores  of  transient  sequences  and  recorded  on  videotape 
several  of  them  for  varying  crystal  sample,  incident  intensity,  and  magnitude  and 
polarity  of  applied  voltage. 

Figs  la-f  show  a  transient  sequence  for  one  crystal.  This  was  a  sample  from 
Sumitomo  Electrics.  Its  front  face  was  a  1  cm  square  and  it  measured  3mm  thick. 
A  voltage  of  5  kV  was  applied  along  the  001  axis  and  the  incident  light  was  vertically 
polarised,  of  wavelength  514  nm,  and  intensity  0.3  mW/cm2.  Fig  la  starts  with  the 
transmitted  light  with  no  applied  field.  Scatter  is  seen  predominantly  on  sides  of  the 
crystal  image.  The  top  and  bottom  scatter  regions  are  due  to  an  incomplete  fit  of  the 
crystal  in  its  holder — there  were  gaps  through  which  light  could  pass.  It  is  likely  the 
scatter  on  the  left  results  from  a  reflection  off  an  internal  surface,  although  we  took 
efforts  to  ensure  the  incident  light  was  normal  to  the  front  surface.  The  subsequent 
photographs  show  the  emergence  of  the  light  pulse  on  the  left  side  at  intervals  of  20 
ms.  The  pulse  manifests  itself  as  somewhat  staccato  horizontal  streamers  that  reach 
out  from  the  main  transmitted  beam.  There  are  two  counterpropagating  pulses — 
one  on  either  side  of  the  crystal,  but  the  one  corresponding  to  the  primary  direction 
of  energy  transfer  is  the  most  prominent. 

We  can  explain  this  pulse  in  terms  of  amplification  of  scatter  noise.  As  soon 
as  a  voltage  is  applied  to  the  crystal,  power  can  be  transferred  to  scattered  beams 
making  small  angles  with  the  incident  light.  These  beams  are  amplified,  generating 
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Figure  1:  Transient  effects  in  Sumitomo  BiuSiOjo  crystal  sample,  (a)  No  applied 
field,  (b)  20  ms  after  application  of  field,  (c)  40  ms.  (d)  60  ms.  (e)  80  ms.  (f)  100 


ms. 
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Figure  2:  Transient  effects  in  Crystal  Technology  BiuSiOjo  crystal  sample,  (a)  No 
applied  field,  (b)  Field  turn  on.  (c)  20  ms  after  application  of  field,  (d)  60  ms.  (e) 
80  ms.  (f)  100  ms. 


CHAPTER  6:  CRYSTAL  SAMPLE  WITH  A  WEDGE  ANGLE 


139 


\ 


(g) 


00 


0) 


(j) 


Figure  2:  Transient  effects  in  Crystal  Technology  Bi]2Si02o  crystal  sample,  (g)  120 
ms.  (h)  160  ms.  (i)  200  ms.  (j)  400  ms. 
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more  scatter  which  can  be  amplified,  creating  additional  scatter  and  so  on.  Power 
is  passed  on  like  the  baton  in  a  relay  race.  This  transfer  is  very  effective  in  the 
transient  phase  because  of  the  presence  of  moving  fringes  [49]. 

Through  our  studies  we  can  make  these  observations.  Keeping  everything  else 
constant: 

•  Pulse  duration  increases  with  decreasing  incident  illumination. 

•  Pulse  duration  and  spatial  extent  decreases  with  decreasing  applied  voltage. 

•  The  pulse  is  seen  predominantly  on  the  side  corresponding  to  the  direction  of 
normal  power  transfer. 

•  The  pulse  reoccurs  when  the  voltage  is  turned  off. 

Figs  2a-j  show  the  Crystal  Technology  Bi^SiOzo  sample  with  the  wedge  angle 
(0.036°).  Here  the  bright  pulse  of  light  moves  horizontally  to  the  right  (the  correct 
direction  for  normal  power  transfer  for  this  orientation),  and  we  see  the  higher 
diffracted  orders  emerge  on  either  side  in  sequence.  This  is  a  very  pretty  thing  to 
watch.  The  intensity  of  the  orders  undulates  as  another  pulse  moves  out  along  the 
35°  line. 

Clearly  this  phenomenon  gives  us  rare  insight  into  energy  transfer  among  higher 
diffracted  orders.  It  is  worth  looking  at  this  in  greater  detail. 

6.1.3  Photodiode  measurements 

A  closer  inspection  of  the  transient  behaviour  was  made  using  the  experimental 
layout  in  Fig  3.  The  path  length  from  crystal  to  photodetector  was  approximately 
3  meters.  This  length  was  necessary  to  provide  adequate  spatial  separation  of  the 
orders  to  steer  them  to  the  appropriate  detector.  A  storage  oscilloscope  monitored 
two  orders  at  a  time.  We  numbered  the  orders  as  in  Fig  4.9a,  that  is  the  main 
incident  beam  was  the  0  order,  the  first  multiple  reflection  (on  top)  -1,  the  second 
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Figure  3:  Study  of  transients  in  higher  diffracted  orders.  Experimental  layout. 

-2,  etc.  Those  orders  below  the  0  order  wave  were  numbered  4-1,  4-2,  and  so  on. 
The  -1  order  was  always  monitored  and  this  provided  a  reference.  An  attempt  to 
measure  the  main  beam  revealed  no  discernible  change  in  intensity  during  transient 
effects,  ie.  an  undepleted  pump.  Repeatibility  was  excellent  in  these  runs. 

Fig  4  shows  some  typical  data.  This  is  for  an  applied  voltage  of  4  kV  and 
incident  intensity  of  1.4  mW/cm2  for  beams  -1  and  4-1.  We  see  that  both  beams 
gain  in  intensity  after  the  field  is  applied,  but  beam  4-1  gains  more  slowly  and  peaks 
while  beam  -1  is  losing  strength.  Both  beams  show  ringing  effects,  but  in  beam  4-1 
the  oscillations  are  more  damped. 

Fig  5  shows  oscilloscope  traces  for  beams  4-1,  4-2,  -1  and  -2  (the  only  beams 
measurable)  for  the  same  intensity  of  0.45  mW/cm2  and  an  applied  field  of  4  kV. 
This  is  a  composite  of  3  plots  in  which  the  -1  traces  were  matched.  The  dashed 
line  marks  the  instant  the  voltage  is  first  applied.  The  important  feature  here  is 
the  confirmation  of  what  was  observed  in  the  video  runs.  The  intensity  peaks  move 
out  sequentially  from  the  main  beam.  Also  higher  diffracted  orders  seem  to  grow 
in  pairs — orders  4-1  and  -2  peak  at  roughly  the  same  time. 

We  can  explain  the  relative  delay  in  power  transfer  to  the  outlying  orders  in  two 
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Figure  4:  Simultaneous  oscilloscope  traces  for  photodiode  signals  of  beams  -1  and 
+  1.  Incident  intensity  is  1,4  mW/cm2.  Applied  voltage  is  4  kV.  Dashed  line  repre¬ 
sents  the  instant  voltage  is  first  applied. 

ways.  First,  efficient  power  transfer  depends  on  a  strong  modulated  space  charge 
field.  The  coupling  of  an  order  to  its  adjacent  orders  is  expressed  in  the  modulation 
ratio,  m,  in  Eqn  1.  Until  the  adjacent  order  amplitude  grows  large,  the  driving 
term,  /m,  is  weak.  Thus  the  growth  of  one  order  couples  it  more  strongly  tc  the 
next  order  and  so  on.  This  is  why  optical  intensity  seems  to  be  handed  from  order 
to  order.  The  second  effect  is  due  to  coupling  with  the  main  beam.  The  larger  the 
interbeam  angle,  the  longer  the  grating  writing  time.  So  receiving  power  directly 
from  the  main  beam  takes  more  time  for  distant  orders. 

We  can  take  the  time  from  first  applying  the  voltage  to  the  first  peak  in  intensity 
as  a  repeatible,  representative  time  value  for  the  photorefractive  transient  effect. 
Because  we  are  using  low  light  intensities  the  transient  behaviour  has  time  constants 
on  the  order  of  tenths  of  seconds  or  longer.  The  time  constants  for  the  rise  of  the 
voltage  delivered  by  the  power  supply  was  on  the  order  of  milliseconds — so  we  may 
neglect  transient  circuit  effects.  Fig  6a  shows  the  reciprocal  of  this  time  value 
plotted  as  a  function  of  incident  intensity  for  the  -1  beam.  The  fit  is  remarkably 
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Figure  5:  Composite  of  oscilloscope  traces  for  photodiode  signals  of  beams  -1,-2, 
+  1  and  -f2.  Incident  intensity  is  1.45  mW/'cm2.  Applied  voltage  is  4  kV.  Dashed 
line  represents  the  instant  voltage  is  first  applied. 
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Figure  6:  Reciprocal  of  time  from  the  instant  voltage  is  first  applied  to  the  first 
intensity  peak  vs  incident  intensity.  Applied  voltage  is  4  kV.  (a)  -1  order,  (b)  +1 
order. 
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Figure  6:  Reciprocal  of  time  from  the  instant  voltage  is  first  applied  to  the  first 
intensity  peak  vs  incident  intensity.  Applied  voltage  is  4  kV.  (t)  -2  order. 

linear.  This  suggests  the  time  value  is  in  some  way  proportioned  to  the  dielectric 
relaxation  time  constant,  rd,  which  we  know  is  inversely  proportional  to  intensity. 
Figs  6b  and  c  show  similar  plots  for  orders  +1  and  -2.  Here  the  slopes  of  the  best 
linear  fits  are  very  close  (6.6  and  7.3  mW/cm2-s  respectively),  which  again  suggest 
the  emergence  of  these  orders  as  a  pair. 

Finally  we  need  to  ask  what  is  the  dependence  of  the  transient  behaviour  on 
applied  field.  In  two- wave  mixing  an  increase  in  applied  field  tended  to  slow  down 
the  transient  effects.  This  is  seen  in  higher  orders  as  well,  as  evidenced  in  Fig  7. 

Fig  8  shows  the  effect  of  applied  voltage  on  the  reciprocal  time  from  the  instant 
the  voltage  is  first  applied  to  the  first  intensity  peak.  Again  we  see  a  slowing  of 
transient  behaviour  for  large  applied  electric  fields. 

Our  wedge  angle  has  provided  a  simple  means  of  producing  highly  repeatible 
transient  effects  in  higher  diffracted  orders. 


Figure  7:  Simultaneous  oscilloscope  traces  for  photodiode  signals  of  be*“s 
-2  for  various  applied  voltages.  Incident  intensity  is  2.6  mW/cm  .  (a) 

2.5  kV.  (c)  0.5  kV. 
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Figure  8:  Reciprocal  of  time  from  the  instant  voltage  is  first  applied  to  first  intensity 
peak  vs  applied  voltage.  Incident  intensity  is  2.6  mW/cm2. 

6.2  Application  to  Laser  Beam  Steering 

6.2.1  Probe  of  the  grating 

One  of  the  exciting  prospects  of  wave  mixing  using  a  crystal  with  a  wedge  angle  is 
the  possibility  of  writing  holograms  with  very  large  grating  spacings.  This  extreme 
drift-dominated  regime  has  never  been  explored  and  might  have  some  interesting 
features.  On  an  optical  table  interbeam  angles  are  limited  by  the  geometry  and 
size  of  optical  elements  such  as  prisms  and  mirrors.  The  smallest  in-air  angle  we 
could  achieve  was  %  1.2°.  The  wedge  gave  an  effective  in-air  angle  between  beams 
of  0.19°  which  results  in  a  grating  spacing  of  150  fim. 

Our  crystal  with  a  wedge  effectively  produces  very  small  grating  vectors.  We 
decided  to  write  a  grating  using  a  514  nm  beam  and  probe  collinearly  with  a  633 
nm  beam.  This  probe  beam  will  not  satisfy  the  Bragg  condition  exactly,  however 
the  angles  involved  are  sufficiently  small  that  efficient  diffraction  may  be  expected. 
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Figure  9:  Probing  the  grating.  The  wavelengths  of  the  writing  (A,)  and  reading 
(A,)  beams  are  different. 

Fig  9  shows  that  the  angle  of  the  diffracted  probe  beam  is 


Here  the  subscript  r  refers  to  the  reading  beam  and  g  to  the  grating  writing  beam. 
For  Bi12SiOJ0  the  indices  of  refraction  at  514  and  633  nm  are  2.62  and  2.52  re¬ 
spectively.  So  we  expect  9r  =  1.28  0,.  We  should  see  diffracted  633  nm  beams  at 
different  angular  spacings  than  the  514  nm  orders. 

The  idea  of  a  collinear  probe  beam  which  satisfies  the  Bragg  condition  was  first 
investigated  by  Huignard  and  Ledu  [55].  This  is  discussed  also  in  the  article  by 
Mallick,  et  al.,  on  spatial  subharmonics  [89].  The  K / 2  grating  vector  (for  their  case 
a  grating  spacing  of  80  fim)  produced  by  light  of  wavelength  514  nm  was  read  by 
light  of  wavelength  1.06  /zm  and  part  of  the  reading  beam  was  diffracted.  This  is 
significant  because  of  the  importance  of  the  near-infrared  in  optical  communica¬ 
tions.  Having  the  writing  (514  nm)  and  reading  (1.06  jxm)  beams  collinear  might 
be  necessary  if  the  beams  are  restricted  to  the  same  optical  path,  say  along  an 
optical  fibre.  The  photorefractive  crystal  could  act  as  a  switch  to  diffract  the  signal 
beam,  depending  on  the  presence  of  the  grating  writing  beam  and  other  control 
parameters.  Our  scheme  might  be  called  collinear  near-Bragg  diffraction. 


Diff  ractea 
reac  beam 
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Figure  10:  Use  of  a  633  nm  probe  beam  to  read  a  grating  formed  by  514  nm  light. 
Probe  and  input  writing  beam  are  collinear. 

Fig  10  shows  the  experimental  layout  for  this  test.  An  expanded,  collimated, 
uniform  intensity  beam  of  wavelength  514  nm  was  at  normal  incidence  on  the  crys¬ 
tal.  A  beam  splitter  allowed  the  injection  of  a  collinear,  unexpanded  beam  from  a 
He-Ne  laser.  (An  attenuating  filter  to  used  to  reduce  the  intensity  of  the  633  nm 
in  some  experimental  runs,  but  the  results  were  the  same  as  for  an  unfiltered  probe 
beam.)  A  narrow  band  filter  was  used  behind  the  crystal  to  remove  the  green  light 
for  many  of  the  photographs. 

We  did  not  observe  collinear  near-Bragg  diffraction.  The  results  were  much  more 
surprising  and  suggest  a  novel  means  of  laser  beam  steering. 

6.2.2  Control  of  direction  of  scattered  beams 

Figs  11  and  12  show  the  results.  Fig  11a  depicts  the  main  transmitted  beam  and 
its  first  multiple  reflection  for  no  applied  field.  Fig  lib  is  the  multiple  reflection 
pattern  of  the  He-Ne  beam  with  no  voltage  applied  to  the  crystal.  Fig  11c  shows 
both  beams.  When  a  voltage  was  applied  both  beams  produced  higher  diffracted 
orders.  However,  the  red  orders  overlaid  the  green.  There  was  no  discernible  reading 
of  the  514  nm  grating  by  the  633  nm  light.  The  red  higher  diffracted  orders  seem 


Figure  12:  Apparent  rotation  of  line  of  diffracted  orders  for  applied  voltage  of  5 
kV.  Change  in  writing  beam  intensity:  (a)  Change  in  writing  beam  intensity:  0.2 
mW/cm3.  (b)  2  mW/cm2.  (c)  4  mW/cm1.  (d)  5  mW/cm3.  (e)  9  mW/cm3.  (f)  11 
mW/cm3.  (g)  15  mW/cm3.  (h)  20  mW/cm  2.  (i)  30  mW/cm2. 
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to  arise  from  the  same  mixing  mechanism  as  the  green  ones — via  multiple  reflection 
due  to  the  wedge  angle.  This  is  somewhat  unexpected  because  BiijSiCKo  is  relatively 
insensitive  at  red  wavelengths.  For  example,  the  quantum  efficiency  (a  measure  of 
carrier  production  per  absorbed  photon)  is  an  order  of  magnitude  less  at  633  nm 
than  at  514  nm  431.  However  photorefractive  gratings  at  633  nm  in  Bi^SiOio  have 
been  measured  by  Powell  and  Petts  103!.  They  found  the  sensitivity  in  the  red 
could  be  increased  by  a  factor  of  16  if  the  crystal  was  heated  to  300°  C. 

Figs  12a-i  show  the  effect  of  increasing  the  intensity  of  the  green  light,  with  a 
constant  applied  voltage  of  5  kV.  The  green  light  is  filtered  out  for  these  figures.  The 
line  of  diffracted  orders  appears  to  rotate.  Actually  the  bottom  diffracted  orders 
disappear  and  the  top  orders  move  horizontally  at  spatial  rates  which  preserve  the 
line.  This  is  a  slow  process,  taking  from  10  seconds  to  over  a  minute  to  reach  steady- 
state.  The  effect  can  be  observed  as  well  by  keeping  the  writing  beam  intensity 
constant  and  by  varying  voltage.  If  either  the  intensity  is  suddenly  turned  down  or 
the  applied  voltage  removed,  the  diffracted  beams  slowly  drift  back  to  their  original 
positions  (as  in  Fig  lib). 

We  can  easily  exploit  this  effect  for  laser  beam  steering.  We  know  the  higher 
diffracted  orders  are  “copies”  of  the  input  beam  (in  the  URR  they  had  the  same 
mode  structure  as  the  oscillating  beam).  By  adjusting  either  the  applied  voltage 
or  the  intensity  of  the  writing  beam,  we  can  direct  the  a  copy  of  the  633  signal 
beam.  Figs  13a  and  b  show  the  apparent  rotation  of  the  line  of  diffracted  orders  as 
a  function  of  writing  intensity  and  applied  voltage.  The  angle  of  rotation  appears 
to  be  linear  with  intensity  and  quadratic  with  appliea  voltage. 

How  can  we  explain  the  movement  of  the  diffracted  orders?  It  is  unlikely  that  it 
is  due  to  a  piezoelectric  deformation  of  the  crystal.  Piezoelectricity  is  not  intensity 
dependent.  There  is  a  slight  fracture  in  the  crystal.  Could  the  presence  of  an  applied 
field  and/or  high  intensity  light  cause  expansion  of  the  fissure?  It  is  possible,  but 
close  inspection  of  the  diffracted  orders  using  just  514  nm  light  suggests  this  cannot 
be  the  mechanism  responsible.  When  a  voltage  is  applied,  increasing  the  incident 
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intensity  also  shifts  these  diffracted  orders.  The  outlying  orders  tend  to  smear  out 
into  regions  of  muted  noise  and  only  the  -1  order  persists.  This  tends  to  move 
horizontally,  but  as  Figs  2g-i  reveal  there  is  still  a  vestige  of  the  original  order  left 
behind  (seen  in  the  photograph  as  a  crescent  region).  The  responsible  mechanism 
must  depend  on  both  intensity  and  applied  electric  field.  It  is  characterised  by  time 
constants  much  longer  than  those  observed  for  the  transient  effects. 

What  is  particularly  interesting  about  this  new  phenomenon  is  the  integrity  of 
the  spatially  shifted  orders.  As  discussed  above,  the  -1  order  using  514  nm  light 
leaves  a  ghost  in  the  original  location  (the  position  corresponding  to  no  applied 
voltage).  But  with  633  nm  light  the  shift  seems  to  be  complete — the  entire  beam 
moves.  Fig  12i  shows  the  -1  order  is  stretched  out  somewhat  for  very  intense  incident 
beams.  It  is  interesting  to  note  that  the  -2  order  moves  twice  as  far  as  the  -1  order, 
thus  preserving  the  line  of  diffracted  orders. 

We  will  discuss  two  possible  mechanisms.  The  first  is  a  thermal  effect.  If  a 
temperature  gradient  or  a  region  of  inhomogeneity  exists  in  the  crystal,  the  refrac¬ 
tive  index  may  vary  spatially.  This  could  divert  the  beams,  and  the  heating  might 
well  depend  on  incident  intensity  and  applied  voltage.  Such  an  effect  would  tend  to 
preserve  the  line  of  orders.  Because  the  angles  involved  are  small,  all  the  beams  we 
observe  travel  through  the  same  part  of  the  crystal,  therefore  the  effect  of  a  nonuni¬ 
form  index  will  be  felt  by  all.  The  main  beam  passes  through  the  crystal  once, 
the  first  multiple  reflection  makes  two  additional  transits,  and  the  second  multiple 
reflection  two  trips  beyond  that.  Thus  we  would  expect  an  increased  deviation  for 
the  higher  orders.  If  this  is  the  mechanism  then  this  effect  will  differ  from  crystal 
sample  to  crystal  sample,  and  perhaps  even  within  a  single  crystal  sample  if  the 
probe  beam  passes  through  a  different  site. 

The  second  possible  mechanism  is  that  of  noise  gratings.  The  growth  and  effi¬ 
ciency  of  noise  gratings  depends  on  intensity  and  strength  of  applied  electric  field, 
and  they  emerge  on  time  scales  similar  to  those  describing  the  movement  of  the 
diffracted  orders.  The  effect  we  see  may  be  considered  similar  to  that  which  causes 
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the  horizontal  pulse  of  light  energy  at  field  turn-on  discussed  in  Sec  6.1.2.  That 
is,  scattered  light  is  amplified  in  a  preferred  direction.  This  effect  is  well-known 
in  BaTiOs  and  is  called  fanning  [341.  although  it  has  never  been  reported  in  the 
literature  on  Bi12SiO20.  In  BaTiOs  beams  interacting  in  the  crystal  can  suffer  sig¬ 
nificant  bending,  particularly  in  the  case  of  a  self-pumped  phase  conjugate  mirror 
using  internal  reflections  [27,28). 

Another  remarkable  feature  is  that  the  movement  is  in  the  opposite  direction 
to  that  expected  for  normal  power  transfer  for  vertically  polarised  light.  Foote  and 
Hall  have  demonstrated  that  the  direction  of  power  transfer  depends  on  the  polar¬ 
isation  of  the  mixing  waves  in  Bii2Si02o[39].  In  the  presence  of  an  applied  electric 
field  Bi12SiO20  is  birefringent,  and  for  light  incident  upon  the  110  face,  two  orthogo¬ 
nally  polarised  states  exist  in  the  medium.  These  two  states  correspond  to  different 
directions  for  power  transfer.  However,  birefringence  cannot  be  the  complete  solu¬ 
tion  to  the  wayward  spatial  shift  of  the  orders.  We  would  expect  to  see  a  splitting 
of  the  diffracted  orders,  not  a  unidirectional  movement.  Further,  measurements  of 
the  polarisation  states  of  the  diffracted  orders  show  that,  as  expected,  they  are  el- 
liptically  polarised,  primarily  in  a  direction  commensurate  with  incident  vertically 
polarised  light.  Finally,  only  the  diffracted  orders  move,  the  main  transmitted  0 
order  does  not  move. 

This  is  a  very  interesting  phenomenon  and  one  which  merits  further  research. 


6.3  Conclusions 

Our  crystal  sample  with  a  wedge  angle  has  provided  a  wealth  of  novel,  interesting 
effects  to  observe.  In  this  chapter  we  concentrated  on  two  effects  with  particular 
applications.  The  first  concerned  transient  power  transfer  as  a  voltage  was  applied  to 
the  crystal.  A  first-ever  study  of  transient  effects  in  higher  diffracted  orders  revealed 
that  light  energy  is  passed  on  from  order  to  order  sequentially.  The  orders  seem 
to  emerge  in  pairs.  Finally  we  observed  a  new  intensity/electric  field  dependent 
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scattering  in  BiuSiOio  reminiscent  of  beam  bending  in  BaTi03.  This  last  effect 
allows  for  angular  control  of  one  or  more  of  the  higher  diffracted  orders. 


Chapt  er  7 

Conclusions  and  Suggestions  for 
Future  Work 


Let  us  not  go  over  the  old  ground,  let  us  rather  prepare  for  what  is  to  come. 

Cicero 


What  can  I  know  ?  What  ought  I  to  do?  What  may  I  hope  ? 


Immanuel  Kant 


A  multitude  of  interesting  and  important  optical  effects  arise  when  beams  of  light 
mix  and  interact  inside  a  photorefractive  crystal.  Through  the  first  six  chapters  and 
three  appendices  of  this  thesis,  we  have  investigated  a  wide  variety  of  some  of  these 
phenomena.  Our  main  interest  has  been  what  photorefractive  materials  do  best. 
They  provide  a  medium  for  the  coupling  of  beams.  Like  some  agent  or  broker, 
they  arrange  for  trades  in  energy  or  differential  rates  in  the  acquisition  of  phase. 
Photorefractive  crystals  can  achieve  all  this  in  real  time  and  in  a  controllable  fashion. 
This  remarkable  ability  stems  from  the  nonlinear  nature  of  the  polarisation  of  the 
medium  in  the  presence  of  a  light  wave. 

This  last  chapter  consists  of  three  parts.  The  first  is  a  brief  survey  of  some 
of  the  accomplishments  described  in  this  thesis.  The  second  is  a  guide  for  further 
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studies  in  photorefractives.  either  here  at  Oxford,  or  elsewhere.  Finally  we  sav  a 
few  encouraging  \  ords  about  the  future  of  photorefractive  crystals. 


7.1  Accomplishments 

Ignoring,  for  the  moment,  the  sage  advice  of  Cicero  at  the  beginning  of  this  chapter, 
let  us  make  a  quick  review  of  the  key  points  of  this  thesis.  This  will  help  define  the 
directions  of  future  work. 

We  began  with  a  description  of  the  photorefractive  effect  and  the  place  of  pho¬ 
torefractive  crystals  in  the  long,  proud  history  of  crystals.  This  led  to  a  discussion 
of  the  evolution  of  our  current  understanding  and  applications  of  photorefractivity. 

The  nonlinear  effects  observed  in  these  crystals  are  numerous,  but  they  are 
basically  variations  on  a  theme.  The  theme  is  two- wave  mixing,  or  a  superposition 
thereof.  Chapter  2  discussed  this  fundamental  interaction,  particularly  in  terms  of 
differences  in  the  relative  sizes  of  the  characteristic  fields,  and  the  optimising  of 
power  transfer  between  the  mixing  beams. 

The  preceding  provided  a  necessary  background  for  the  understanding  of  novel 
phenomena  and  techniques  discussed  in  the  remaining  thesis. 

Two-wave  mixing  for  counterpropagating  beams  was  the  subject  of  Chapter  3. 
We  demonstrated  that  for  the  important  cases  involving  moving  intensity  fringes  or 
in  the  presence  of  an  applied  voltage,  the  degree  of  beam  coupling  varies  spatially. 
All  previous  treatments  of  this  problem  assumed  a  constant  coupling  term.  A  com¬ 
parison  of  the  solutions  of  the  coupled  wave  equations  for  constant  and  variable 
coupling  revealed  several  cases  in  which  the  results  differ  significantly.  Our  numer¬ 
ical  method  was  employed  next  to  make  the  first  systematic  theoretical  assessment 
of  the  control  of  nonreciprocal  transmission  and  phase  shift  in  LiNbC>3,  a  represen¬ 
tative  photorefractive  crystal.  For  both  NR  effects,  the  velocity  of  the  light  intensity 
fringes  offers  the  greatest  range  for  control.  This  is  particularly  attractive  because 
the  photocurrent  in  the  circuit  containing  the  crystal  was  shown  to  depend  on  the 
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degree  of  wave  interaction.  The  photocurrent  has  its  minimum  values  when  beam 
coupling  is  strongest.  This  should  provide  a  useful,  nonintrusive  means  to  monitor 
NR  transmission.  A  comparison  of  the  theoretical  results  for  different  crystal  types 
suggests  both  LiNb03  and  BaTi03  should  show  strong  NR  effects,  however  our 
LiNb03  crystal  samples  proved  too  erratic  in  coupling  to  allow  us  to  measure  any 
effects.  Our  final  conclusion  about  using  photorefractive  crystals  for  NR  purposes 
was  that  they  look  promising  for  NR  transmission  (eg.  as  an  optical  diode),  but 
not  for  NR  phase  shift.  Moving  intensity  fringes  and  an  incident  beam  ratio  much 
different  from  1  are  required  to  obtain  a  suitably  large  NR  phase  shift  ( s;  tr/2),  a 
value  readily  realizable  using  NR  magnetooptical  devices. 

In  the  next  chapter  we  began  a  study  of  oscillation  in  resonators  with  pho¬ 
torefractive  gain.  We  looked  at  the  unidirectional  ring  resonator  and  made  several 
new  observations.  We  found  a  threshold  for  incident  light  intensity  below  which 
no  oscillation  occurred  in  Bi12Si020.  For  our  sample  this  was  0.8  mW/cm2.  We 
measured  the  frequency  detuning  of  the  oscillating  beam  (relative  to  the  incident 
light)  and  found  no  measurable  difference  between  the  various  resonator  modes.  In 
BaTi03  when  one  slightly  repositions  one  of  the  mirrors  forming  the  ring  resonator, 
the  resulting  change  in  resonator  length  is  accompanied  by  a  change  in  frequency 
detuning  of  the  oscillating  beam.  We  discovered  this  was  not  the  case  for  Bii2SiO20. 
Here  the  detuning,  which  is  a  function  of  light  intensity,  applied  electric  field,  and 
interbeam  angle,  must  be  close  to  an  optimum  value  to  achieve  the  large  gain  neces¬ 
sary  for  oscillation.  A  surprising  new  result  was  a  preponderance  of  higher  diffracted 
orders  for  one  crystal  sample.  Many  of  these  additional  beams  were  due  to  a  slight 
wedge  angle  in  the  crystal.  Different  diffracted  orders  (which  we  dubbed  compan¬ 
ion  orders)  arose  from  two-wave  mixing  between  a  beam  in  the  crystal  and  its  first 
codirectional  multiple  reflection. 

Chapter  5  is  significant  because  it  provides  an  experimental  verification  of  a 
new  gain  mechanism  for  photorefractive  materials.  The  new  amplification  scheme 
involves  off-Bragg  coupling  and  results  in  a  power  transfer  from  two  pump  beams 
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inward  to  two  symmetric  signal  beams.  The  gain  is  readily  controllable  by  chang¬ 
ing  the  frequency  detuning  of  the  pump  beams.  At  its  optimum  value,  the  new 
gain  mechanism  in  Bi12Si02o  gives  slightly  larger  gain  than  the  traditional  two- 
wave  mixing  mechanism.  This  was  confirmed  experimentally  in  a  unidirectional 
ring  resonator  where  the  small  increase  in  the  gain  coefficient  using  two  pump 
beams  resulted  in  a  steady-state  oscillating  beam  intensity  four  times  larger  than 
for  a  single  pump.  The  increased  effectiveness  of  resonator  operation  permitted  the 
first  reported  quantitative  measurement  of  oscillation  in  a  linear  resonator  using 
Bi12SiO20.  Finally  we  searched  for  evidence  of  spatial  subharmonics  between  the 
two  pump  beams,  but  found  no  preferential  position  for  the  signal  beams. 

We  then  returned  to  a  study  of  the  crystal  sample  with  the  wedge  angle.  The 
minuteness  of  the  wedge  angle  (0.036°  as  measured  by  the  angular  spacing  of  adja¬ 
cent  orders)  allows  us  to  investigate  the  case  of  very  large  grating  spacings — much 
larger  than  has  been  looked  at  before.  Transient  effects  in  this  region  give  keen 
insight  into  the  transfer  of  power  between  diffracted  orders.  Using  a  single  uniform 
input  light  beam,  the  transient  phase  of  grating  formation  shows  power  is  passed 
sequentially  from  order  to  order.  In  other  words,  one  order  must  be  amplified  before 
it  can  effectively  pass  on  power  to  an  adjacent  order.  We  have  identified  a  transient 
pulse  of  light  which  occurs  when  a  voltage  is  first  applied  to  an  illuminated  Bi12SiO20 
crystal.  All  our  crystal  samples  exhibited  this  pulse  to  some  extent.  As  an  aside, 
we  developed  a  new,  simple,  analytic  formula  which  describes  transient  behaviour 
in  two- wave  mixing  and  which  results  in  a  remarkable  savings  in  computation  time 
over  other  numerical  and  analytic  treatments.  Finally  we  demonstrated  a  novel 
method  of  laser  beam  steering  in  which  the  positions  of  the  higher  diffracted  orders 
of  a  collinear  probe  beam  (of  wavelength  different  than  the  writing  beam)  can  be 
shifted  by  changing  either  the  applied  voltage  or  the  intensity  of  the  writing  beam. 
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7.2  Future  Work 

There  are  many  directions  the  work  described  in  this  thesis  can  go.  We  will  consider 
these  directions  by  subject. 

1.  Nonreciprocal  Effects.  The  theory  we  considered  in  Chapter  3  assumed  small 
modulation  depth  for  the  intensity  fringes  (ie.  small  m).  There  are  some 
cases  of  interest  where  this  approximation  may  not  be  valid,  for  example  when 
the  incident  beams  have  roughly  the  same  intensity  and  strong  absorption  is 
present.  It  would  be  worthwhile  to  do  an  analysis,  such  as  was  done  recently 
by  Au  and  Solymar  [8],  for  mixing  beams  in  the  reflection  geometry  and  large 
modulation  depths. 

Experiments  using  LiNbC>3  as  the  medium  were  unsuccessful,  however  NR 
transmission  effects  have  been  reported  in  BaTiC>3  and  KNb03.  One  of  these 
materials  could  be  used  to  verify  our  theory,  particularly  the  aspects  involving 
control  of  these  effects. 

2.  New  Gain  Mechanism.  As  mentioned  in  Chapter  5,  the  best  way  to  evaluate 
this  new  mechanism  is  to  use  three  or  four  beams,  each  of  which  can  be 
frequency  shifted  relative  to  the  others.  Besides  experimental  confirmation  of 
the  theoretical  results,  these  questions  should  be  looked  at: 

•  What  is  the  relation  between  the  signal  beams?  Specifically,  are  they 
phase  conjugates  or  even  near  phase  conjugates?  On  a  related  theme: 
what  type  of  cross-talk  exists  between  the  signal  beams?  Is  spatial  and/or 
temporal  information  shared  between  them? 

•  What  effect  will  incompatible  frequency  detuning  have  on  the  gain?  By 
incompatible  we  mean  that  the  frequency  differences  between  the  beams 
give  different  speeds  for  the  moving  fringe  patterns  inside  the  medium. 
Our  theory  depended  on  a  scaled  frequency  shift  for  the  beams  to  give 
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all  the  fringe  patterns  the  same  velocity.  However,  if  gain  is  possible  or 
even  enhanced  for  different  combinations  of  frequency  shift,  this  would 
be  significant,  and  might  explain  why  we  observed  oscillation  in  the  ring 
resonator  for  such  a  wide  range  of  frequency  detuning  of  the  driving 
pump  beam. 

•  What  is  the  intensity  profile  of  the  region  of  internal  optical  noise  between 
two  frequency  shifted  pump  beams?  How  does  this  optical  noise  relate  to 
the  spatial  subharmonics  observed  by  Mallick,  et.  al.  [89]?  We  described 
this  phenomenon  in  Section  5.2.5.  The  optical  noise,  which  looks  rather 
like  a  swarm  of  irritated  fireflies,  seemed  concentrated  halfway  between 
the  two  beams,  but  extended  well  out  of  the  plane  of  the  pump  beams.  As 
we  changed  the  frequency  of  the  driving  pump  beam,  the  profile  seemed 
to  stay  the  same,  but  either  dimmed  or  brightened  overall.  We  would 
expect  to  see  the  noise  follow  the  pattern  suggested  by  Fig  5.4,  namely, 
for  a  specific  frequency  detuning  value,  the  most  intense  peak  would  be 
in  the  centre,  but  a  change  in  detuning  should  result  in  a  symmetric, 
bimodal  intensity  distribution  corresponding  to  a  new  optimum  value 
of  the  scaling  parameter,  (Fig  5.2a).  Our  group  recently  ordered  a 
vidicam  and  digitizer,  and  this  equipment  would  be  ideal  for  looking  at 
the  intensity  profiles. 

•  What  would  be  the  result  if  both  signal  beams  were  reinforced  by  (sep¬ 
arate)  unidirectional  ring  resonators?  We  made  a  brief  attempt  at  this 
set-up,  with  minor  success.  The  difficulty  is  insuring  that  one  is  reinforc¬ 
ing  the  second  signal  beam,  and  not  just  an  arbitrary  noise  beam  close 
to  the  symmetric  position.  If  the  two  URRs  are  spot  on,  one  might  be 
able  to  achieve  significantly  enhanced  gain  in  both. 


3.  Linear  Resonator.  The  intensity  of  the  oscillating  beam  in  our  linear  res¬ 
onator  described  in  Chapter  5  was  too  weak  to  generate  a  phase  conjugate 
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beam.  We  recently  had  the  surfaces  of  a  1  cm  cube  of  BitiSiOio  antirefiection 
coated.  With  reduced  reflection  loss  and  a  longer  interaction  length,  amplifi¬ 
cation  may  be  high  enough  to  produce  phase  conjugation.  But  this  raises  an 
interesting  question:  what  type  of  phase  conjugation,  if  any,  is  produced  in  a 
linear  resonator  fed  by  two  input  beams?  With  a  single  input,  the  counter- 
propagating,  oscillating  beams  become  the  “pump"  beams,  and,  if  they  are 
sufficiently  intense,  a  phase  conjugate  of  the  input  beam  emerges.  With  two 
input  beams  do  we  get  two  phase  conjugates?  WhaL  if  the  oscillating  beam 
does  not  bisect  the  angle  between  the  two  input  beams?  If  each  input  gen¬ 
erates  a  phase  conjugate,  and  if  there  is  a  second  signal  beam  between  the 
two  inputs,  then  we  could  have  the  case  of  seven-wave  mixing — a  complicated 
situation  indeed! 

4.  New  Crystal  Cut.  In  Appendix  B  we  suggested  a  new  crystal  cut  which  would 
result  in  an  enhanced  electrooptic  effect  in  BiuSiC^o  for  wave  mixing  in  the 
absence  of  optical  activity.  Optical  activity  will  tend  to  reduce  the  net  pho- 
torefractive  gain,  but  a  question  worth  asking  is,  will  the  overall  gain  for  a 
BiuSiOjo  crystal  cut  with  a  111  face  be  higher  than  for  a  001  face  cut?  The 
difficulty  here  would  be  a  fair  test —  different  crystal  samples  usually  have 
different  crystal  parameters.  We  have  a  rare  opportunity  to  make  a  definitive 
test.  In  Appendix  A  we  described  crystal  damage  in  a  1  cm  cube  of  Bii2SiO20- 
The  undamaged  part  of  the  crystal  should  be  large  enough  to  cut  into  two 
samples  of  equal  length.  Both  samples  should  have  identical  parameters.  The 
original  orientation  and  that  proposed  in  Appendix  B  use  polished  1 10  and 
llO  faces,  so  the  new  crystal  cut  will  merely  involve  a  rotation  of  54.7°  about 
the  normal  axis  through  the  polished  faces. 

5.  Wedge  Angle.  Once  the  test  on  the  new  crystal  cut  is  finished,  a  slight  wedge 
angle  can  be  introduced  into  the  standard  cut  crystal.  It  would  be  interesting 
to  compare  wedge  effects  with  the  sample  we  used  in  Chapter  6.  In  particular 
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we  need  to  see  if  the  spatial  shift  of  the  higher  diffracted  orders  of  a  collinear 
probe  beam  is  reproducible  in  a  different  crystal  sample. 


7.3  Conclusions 

Research  in  photorefractives  is  at  a  very  exciting  stage.  The  field  is  still  young  at 
age  20,  due  in  part  to  a  fallow,  contemplative  phase  in  the  1970s.  But  the  result  was 
a  coherent  charge  transport /field  equation  theory  which  seems  to  explain,  at  least 
qualitatively,  the  effects  seen  in  hundreds  of  experiments.  Now  new  experimental 
results  come  fast  and  furious,  as  if  straining  at  the  leash  of  the  theory,  testing  its 
limits. 

The  next  stage  will  involve  systematic  exploitation  of  the  effects  being  demon¬ 
strated  in  the  laboratory  today.  Many  devices  use  electrooptic  crystals  such  as 
BitjSiOjo  in  spatial  light  modulators  and  LiNbC>3  in  optical  waveguides,  but,  as 
of  this  writing,  no  commercial  products  exist  which  utilize  photorefractive  wave 
mixing.  Still,  the  promise  is  there,  and,  if  the  quality  of  materials  becomes  more 
uniform,  photorefractive  crystals  will  find  their  place  in  optical  interconnects,  res¬ 
onators,  and  computers  of  the  future. 

It  it  our  sincere  hope  that  the  work  represented  in  this  thesis  furthers  the  tran¬ 
sition  to  this  next  stage. 


Appendix  A 

Damage  in  Bi12SiO20  Crystal 
Samples 


Chapters  4,  5,  and  6  chronicle  work  done  with  a  repolished  piece  of  Bii2Si02o  from  a 
damaged  crystal.  The  original  sample  was  5  mm  x  5  mm  x  10  mm.  The  damage  was 
first  observed  as  a  small  fracture  extending  about  1  mm  into  the  crystal,  roughly 
3  ram  from  one  end  on  the  long  edge  which  served  as  the  cathode.  Over  the 
course  of  several  weeks  of  experimental  work  (two-wave  mixing  and  photocurrent 
measurements)  the  fracture  grew  in  extent  and  depth  until  the  crystal  split  in  two. 

At  the  time  we  felt  the  damage  was  initiated  by  an  air  spark  which  passed  over 
the  front  face  of  the  crystal  when  the  flooding  illumination  was  interrupted  and  a  10 
kV/cm  electric  field  was  applied.  Extreme  care  was  taken  with  new  crystal  samples 
to  prevent  a  repeat  occurence.  We  had  limited  success  with  this  as  evidenced  by  Figs 
la  and  lb.  These  photographs  show  a  second  damaged  Bi12Si02o  crystal  sample. 
This  was  a  1  cm  cube  from  Sumitomo,  Ltd.  Fig  la  depicts  the  damaged  crystal 
and  a  new  replacement  crystal  for  comparison.  Fig  lb  shows  details  of  the  internal 
damage.  As  with  the  previous  crystal  the  cracks  originated  on  the  cathode  surface. 

It  is  very  easy  to  chip  the  edges  of  a  Bil2SiO20  crystal.  The  best  way  to  minimize 
edge  damage  is  to  limit  handling  of  the  crystal.  The  mere  act  of  setting  the  crystal 
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Figure  1:  Damage  in  BiuSiOjo  crystal  samples,  (a)  New  Sumitomo  sample  (left) 
and  damaged  Sumitomo  sample  (right),  (b)  Interior  view  of  damaged  sample.  Top 
of  picture  was  crystal  cathode. 
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in  a  holder  with  spring  plunger  electrical  contacts  invariably  causes  some  minor 
chipping.  But  this  does  not  seem  to  be  the  source  of  the  fractures.  The  original 
cracks  on  both  damaged  crystals  started  on  the  sides  of  the  crystals  some  distance 
from  the  edges. 

The  key  to  the  damage  seems  to  be  the  cathode.  Bryksin.  et  al.,  reported  in 
1982  and  1983  the  presence  of  strong  electric  fields  near  the  cathode  in  a  uniformly 
illuminated  sample  of  Bi12Si02o[15,14j.  Very  close  (<  1  mm)  to  the  surface  the 
field  strength  could  be  in  excess  of  100  kV/cm.  Oberschmid  made  a  very  thorough 
study  of  photocurrent  in  Bii2Si02o  and  Bii2GeO20[96].  He  found  that  electrons  were 
injected  into  the  crystal  at  localized  sites  on  the  cathode.  The  current  from  these 
sites  spread  out  along  conical  paths.  The  presence  of  both  positive  and  negative 
charges  at  the  cathode  creates  instabilities  and  Oberschmid  observed  strong  local 
fluctuations  in  the  electric  field.  Anderson  and  Luke  reported  in  1987  beam  induced 
damage  in  Bi13Si02o  in  a  PRIZ  configuration  [4].  The  PRIZ  is  a  spatial  light 
modulator  in  which  the  front  and  real  crystal  faces  are  polished  and  have  transparent 
electrodes.  The  authors  found  fractures  (originating  at  the  cathode)  extending 
about  65  pm  into  the  crystal.  These  arose  using  relatively  low  illuminations  (100 
pW/cm2)  and  weak  applied  fields  (2  kV/cm).  It  appears  damage  manifests  itself 
in  the  form  of  small  pits,  2-3  pm  in  diameter  on  the  surface  of  the  cathode.  Striae 
may  emanate  from  these  pits  and  may  run  between  several  adjacent  pits.  These 
running  striae  can  grow  into  fractures.  It  is  very  likely  that  the  pits  serve  as  the 
local  electron  injection  sites  mentioned  by  Oberschmid. 

What  can  be  done  to  minimize,  or  eliminate,  crystal  damage?  We  used  painted 
silver  electrodes  on  unpolished  crystal  faces.  Surely  the  problem  of  pits  or  rough 
areas  is  worse  on  an  unpolished  surface.  It  might  be  beneficial  to  polish  all  the 
surfaces  of  a  Bij2Si02o  crystal  before  applying  electrodes.  Once  damage  has  been 
identified,  reduced  applied  electric  fields  and  low  illuminations  may  prolong  the 
useful  life  of  the  crystal.  At  least  we  may  state  confidently  the  reverse  is  true — strong 
fields  and  high  intensities  seem  to  aggravate  the  condition.  In  the  Sumitomo  sample 
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we  were  trying  to  obtain  oscillation  in  a  linear  resonator  configuration  (see  Chap  6). 
We  knew  a  fracture  was  present  in  the  crystal  and  we  limited  the  applied  field  to  7 
kV/cm  and  the  intensity  to  10  mW/cm*.  The  crack  extended  approximately  2  mm 
into  the  crystal  and  its  growth  was  slow.  We  increased  the  field  to  10  kV  /cm  and 
the  illumination  to  30  mW/cm2.  The  crystal  suddenly  lit  up  with  bright  internal 
sparks.  We  immediately  turned  off  the  applied  field.  The  fracture  was  now  over 
5  mm  long  and  had  developed  branches.  After  this  episode  the  damage  increased 
visually  each  time  the  crystal  was  used. 

The  internal  sparking  is  a  tell-tale  sign  of  crystal  damage.  In  our  optical  set-ups 
there  is  generally  quite  a  bit  of  scattered  light  falling  on  the  walls  of  the  laboratory 
(the  crystals  send  out  a  lot  of  scatter).  The  first  evidence  of  damage  in  the  Sumitomo 
crystal  was  a  semi-periodic  pulsation  in  a  region  of  scattered  light  when  an  electric 
field  was  applied  to  the  crystal.  When  the  field  was  increased,  the  frequency  of 
the  pulses  increased.  Close  inspection  of  the  crystal  revealed  a  small  area  near  the 
cathode  growing  alternately  bright  then  dark.  This  was  found  later  to  be  the  site 
of  a  small  fracture. 

I  included  our  experiences  with  damage  in  Bii2SiO20  in  a  talk  on  crystal  diag¬ 
nostics  given  at  Cambridge  in  September,  1987.  The  response  of  the  audience  was 
rather  interesting.  My  confession  of  destroying  the  crystals  prompted  several  others 
to  bare  their  souls.  One  individual  exclaimed,  almost  with  pride,  that  his  group  had 
ruined  two  Bil2Si02o  crystals — one  having  literally  exploded  when  an  a.c.  electric 
field  was  applied  in  the  kilohertz  range. 

At  a  previous  scientific  meeting  we  heard  the  woeful  tale  of  a  research  student 
and  a  crystal  of  BaTiOa.  This  crystal  has  a  Curie  temperature  around  5°C.  Care 
must  be  taken  to  avoid  abrupt  temperature  changes  near  this  value.  The  well- 
meaning  student  was  cleaning  the  surfaces  of  the  crystal  using  methanol.  When  the 
methanol  evaporated  the  temperature  change  was  sufficient  to  shatter  the  crystal. 

Relating  such  experiences  is  invaluable  to  research.  We  often  learn  more  from 
mistakes  than  from  successes.  It  would  be  useful  to  include  a  short  forum,  perhaps 
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at  the  end  of  one  or  more  sesssions.  at  conferences  to  pass  on  warnings  and  pitfalls. 
Barring  this,  experimenters  should  seek  out  fellow  workers  to  discuss  things  to  be 
wary  of. 

Equally  as  important  as  understanding  the  uses  of  photorefractive  materials  is 
to  be  familiar  with  the  abuses  of  photorefractive  materials. 


Appendix  B 

General  Treatment  of  the 
Electrooptic  Effect  in  Bi12Si02Q 


The  appearance  of  the  companion  orders  in  the  ring  resonator,  described  in  Chapter 
4,  prompted  the  following  question:  Is  the  electrooptic  effect  in  Bii;SiO20  enhanced 
for  an  electric  field  in  some  direction  other  than  along  the  001  axis?  If  the  answer 
is  yes  then  the  companion  orders  might  be  due  to  some  favourable  combination  of 
grating  spacing  and  orientation  formed  by  the  principal  beam  and  noise. 

Consider  a  general  grating  formed  by  the  mixing  of  two  waves  in  BiuSiCho-  The 
modulated  space  charge  field  is  given  by  Eqn  2.20.  We  should  mention  that  in  the 
derivation  of  E, we  use  the  projection  of  the  applied  electric  field  in  the  direction 
of  the  grating  vector,  K,  as  E, o.  The  magnitude  of  E,x  will  determine  the  index 
modulation  via  the  electrooptic  effect.  E,i  has  the  following  projections  onto  the 
principal  axes  of  the  crystal  with  no  applied  field 

E,  i  =  Eti  +  Eyii  +  Ezz  (1) 

In  the  presence  of  E,x ,  the  index  ellipsoid  is  transformed  to 

x2  y1  z 2  ... 

"7  i - :  i - -  +  2 rEtxy  4-  2rEyxz  +  2 rExyz  —  1  (2) 

n*  n *  n* 

where  r  =  ra3. 
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Figure  1:  Two-wave  mixing  with  grating  vector  not  coincident  with  applied  field. 

We  want  to  find  the  new  principal  coordinate  system  of  the  crystal.  Rewrite  the 
above  as  the  quadratic  form 

Q  =  xrCx  (3) 

^3  ctE  brE 

where  xT  —  (x,y,z)  and  C  -  CrE  \  arE  with  a  =  EXI'E.  b  =  Ey/E.  c  = 

brE  arE  -4 

n * 

Ei,  E.  and  E  —  Efl  .  To  find  the  new  principal  coordinate  system  we  diagonalize 
Q  and  to  do  that  we  need  the  eigenvectors  of  C. 

The  characteristic  equation  for  C  is 

(1-A)  -2r3E3abc-r2E2{a2  -b7 -c2)[~2-\)  =0  (4) 

This  expression  can  be  simplified  somewhat  if  we  assume  eigenvalues  of  the  form 
-  e ,rE  for  i  =  1,2.3. 

e3  -  (a3  -  b2  ~  c3)e  -  2afcc  *  0  (5) 

Eqn  5  is  a  cubic  equation  in  e  whose  solution  in  terms  of  the  coefficients  is  well 
known,  although  a  bit  complicated  '  104 ! .  The  roots  to  this  equation  give  the  eigen¬ 
values  which  represent  the  lengths  of  the  axes  of  the  transformed  index  ellipsoid. 
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The  eigenvectors.  satisfy 

(C  -  AtI)u>,  =  0  (6) 

The  transformation  from  the  old  primary  axes  system  to  the  new  one  is  made 
by 

x'  =  WTx  (7) 

where  T  means  transpose  and  IT  is  the  matrix  whose  columns  are  the  eigenvectors 
u>, .  In  other  words,  the  eigenvectors  are  the  directions  of  the  new  principal  axes. 
To  find  the  net  electrooptic  effect  in  the  new  axes  system,  transform  the  incident 
polarisation  vector  to  primed  coordinates  and  sum  its  projections  onto  the  new 
principal  axes  of  the  index  ellipsoid. 

Our  case  of  interest  is  shown  in  Fig  1-  The  two  mixing  waves  are  at  near  normal 
incidence  to  the  front  surface,  so  the  grating  vector.  A.  is  parallel  to  this  surface 
and  makes  an  angle  ft  with  the  z  (001 )  axis.  Then  a  =  sin  ft  .  v  2,  b  =  sin  ft  .  v  2.  and 
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Figure  3:  New  principal  axes  system  idark  axes  i  for  electric  field  at  an  angle  of 
54.7°  to  the  original  001  (light)  axis. 

c  =  cos  ft.  The  characteristic  equation  is 

e3  -  e  +  sin2  f 2  cos  ft  —  0  (8) 

Fig  2  shows  the  roots  of  Eqn  8  for  ft  between  0  and  180°.  To  find  the  maximum 
absolute  value  for  a  normalised  eigenvalue  look  at  the  general  solution  to  Eqn  9. 

C3  +  djE  -r  CL2  =  0  (9) 

Following  Ref  1041,  we  introduce  the  quantities  Q  =  —  oj/3,  R  =  a2/ 2.  and  9  = 
arccos(./2/v  Q3).  Then  our  three  roots  are 


=  -2y/Qco*y- 

Ej  =  -2y/(?COS^ 


«3  —  — 2y 'Q 


cos 


6 
3 

9  4-  2tt 
3 

9  4-  47T 


(10) 

(ID 

(12) 


For  our  case  ax  =  -1  and  a2  —  sin2  ft  cos  ft.  If  we  optimise,  sav,  with  respect  to 
ft  then  we  find  we  must  satisfy 

sin2  ft  =  ^  (13) 

O 

or  ft  =  .54.7°  with  a  corresponding  eigenvalue  of  -1.155.  The  other  two  eigenvalues 
are  e2  =  0.577  and  e3  =  0.577.  The  eigenvectors  are 
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The  new  principal  axes  are  as  shown  in  Fig  3.  The  first  thing  we  shouid  notice  is 
that  we  can  have  an  enhanced  electrooptic  effect.  If  light  is  polarised  along  the  c 
axis,  then  the  normalised  electrooptic  effect  is  1.155.  or  over  15%  larger  than  for  the 
case  shown  in  Fig  2.4.  We  shouid  recall  that  we  also  have  a  reduced  applied  field 
component  in  this  case  (  E,0  cos  fi).  To  take  full  advantage  of  this  enhancement,  one 
should  have  electrodes  normal  to  this  optimum  axis,  namely  the  111  or  111  face. 

The  111  cut  for  Bi^SiOio  is  used  for  devices  such  as  spatial  light  modulators, 
as  is  the  cut  shown  in  Fig  2.4  (53,84.4  .  Marrakchi.  et  al.,  claim  that  the  crystal 
cut  in  Fig  2.4  is  the  optimum  one  for  the  electrooptic  effect  in  BinSiOjo  for  wave 
mixing  [901.  Optical  activity  will  have  a  more  significant  effect  on  the  electrooptic 
effect  for  crystals  cut  with  a  111  face.  The  normalised  changes  in  refractive  indicies 
for  the  two  principal  axes  for  light  polarised  in  a  plane  parallel  to  the  110  face  are 
1.15  and  -0.577,  where  for  the  crystal  in  Fig  2.4  these  values  are  1  and  0. 

However,  the  enhanced  electrooptic  effect  cannot  explain  our  companion  orders, 
because  there  the  angle  the  grating  vector  made  with  the  z  axts  was  approximately 
35°.  This  is  not  the  optimum  angle,  and  the  oriental’''"  of  the  new  principal  axes 
for  vertically  polarised  light  and  this  angle  of  grating  vector  gives  a  net  normalised 
electooptic  effect  of  -0.885,  a  magnitude  less  than  that  for  the  case  of  ft  -  0°.  This 
supports  the  claim  that  companion  orders  are  due  to  a  wedge  in  the  crystal  sample. 


Appendix  C 

A  Simple  Analytic  Formula  for 
Transient  Two- Wave  Mixing 


In  this  appendix  we  derive  an  analytic  formula  for  the  complex  optical  amplitude 
of  a  wave  undergoing  transient  amplification  in  two-wave  mixing.  We  show  this 
method  compares  remarkably  well  to  a  previously  reported  numerical  method  [48]. 
This  is  important  because  our  method  leads  to  more  efficient  computation. 

In  the  undepleted  pump  approximation  for  two- wave  mixing  the  space  and  time 
differential  equations  describing  transient  behaviour  are 


dS+ 

dx 

dE.  t 
dt 


-jT'E.r 

-gE.  1  -r  fSl 


(1) 

(2) 


Here  S+  —  E+/  v  I0  where  I0  is  the  pump  beam  intensity,  P  =  ^fc0n3refT>  an(^  9  and 
/  are  given  in  Chapter  6. 

Our  situation  is  this:  a  photorefractive  crystal  is  uniformly  illuminated  by  a 
single  beam.  An  externally  imposed  d.c.  electric  field  is  present  in  the  material.  At 
time  t  =  0  a  second  mixing  beam  is  introduced.  The  numerical  method  of  Heaton 
and  Solymar  [48]  to  solve  Eqns  1  and  2  for  this  situation  is  as  follows.  At  t  =  0,  5* 
and  E,i  are  constant— 5^.  =  5+0  and  Etl  =  0  for  all  x.  For  a  short  time  interval 
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-If.  £_  rnav  be  considered  unchanged.  the  space  '.marge  Held  grows  because  <u 
the  presence  of  a  modulated  intensity  pattern.  £,u  At  1  is  calculated  :rom  Eqn  - 
and  this  is  inserted  into  Eqn  i.  The  space  equation  is  integrated  yielding  '.aiues  tor 
S.(x.  At).  This  can  be  inserted  into  Eqn  2  for  the  next  time  step  and  the  process 
see-saws  between  time  and  space  solutions  until  the  required  number  of  time  steps 
are  taken.  For  the  case  of  small  interbeam  angles  and  large  applied  Helds  I  drift 
dominated  case)  the  solutions  are  highly  oscillatory  and  many  steps  are  required  to 
prevent  instabilities.  Even  with  very  small  time  steps  (on  the  order  of  the  dielectric 
time  constant)  our  experience  has  shown  the  intensity  solution  may  go  negative 
(a  unreal  physical  situation l  for  strong  coupling  and  large  interaction  lengths.  In 
other  words,  this  numerical  treatment  is  prone  to  give  solutions  which  are  more 
oscillatory  than  other  methods  [49). 

Our  method  is  similar  in  philosophy,  we  attack  one  differential  equation,  solve 
it,  and  use  the  solution  in  the  other  equation.  Basically  ours  is  an  iterative  method. 
We  start  with  Eqn  2.  We  assume  5*  is  constant  and  find 

E.i(z,t)  =  E.1(x.Q)e-3t  *  -e"**)  (3) 

9 

but  £",1(1.0)  =  0.  We  put  this  solution  into  the  space  equation 

^  =  - j  T'  C'  S+.  (4) 

dx 

C  =  —  (1  —  e~3t)  (5) 

9 

This  has  a  simple  exponential  solution.  In  theory  we  could  continue  the  iterative 
process  going  back  and  forth  between  differential  equations,  but  further  iterations 
will  require  numerical  solutions.  Let  us  stop  here  and  see  how  good  our  answer  is. 
In  terms  of  intensity  and  phase  we  have 


/+(i ,t)  =  /,.0  exp  (2  T'  Im(C)  x) 

<Mz,0  =  r'Re(C)x 


(6) 

(7) 
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Figs  la-d  show  a  comparison  between  the  analytic  and  numerical  methods  for 
interbeam  angles  in  Bi12SiO20  of  10°, 4°,  and  0.5°.  The  crystal  parameters  are  taken 
from  Table  2.1  except  that  in  Fig  lc,  n  =  10-s  m2/ Vs,  and  for  the  first  three  plots 
r  =  10~12  m/V.  In  the  last  plot  r  is  taken  as  3.4  x  10-12  m/V.  The  crystal  length  is 
1  cm,  I0  =  100  W/m2,  the  wavelength  is  514  nm,  the  crystal  temperature  is  300  K, 
and  there  is  no  frequency  detuning  between  the  beams.  For  the  interbeam  angle  of 
10°  the  applied  field  is  10®  V/m  and  for  the  other  plots,  2  x  106  V/m. 

The  agreement  between  the  two  methods  is  very  good,  regardless  of  the  inter¬ 
beam  angle,  for  the  case  of  weak  coupling  (small  value  of  r  or  small  interaction 
lengths).  For  strong  coupling  the  two  methods  still  agree  for  large  angles  (>  8°) 
but  at  smaller  angles  the  oscillations  for  the  analytic  method  are  more  damped. 

The  great  savings  is  in  computation  time.  Calculating  and  storing  220  time 
points  using  the  analytic  formula  took  less  than  5  seconds  on  an  IBM  personal 
computer.  The  numerical  method  needed  20  minutes  to  compute  the  same  number 
of  points. 

Finally  we  show  that  the  analytic  expression  Eqn  4  is  the  same  if  we  start  the 
iteration  with  Eqn  1  and  take  one  more  iterative  step.  This  fact  substantiates  the 
observation  that  Eqn  4  is  a  reasonable  approximation  to  the  transient  behaviour 
calculated  by  a  more  thorough  treatment. 

The  solution  to  Eqn  1  assuming  E,\  constant  is 


S+(x,t)  =  5+0exp(-;T'E#i*x) 

If  we  put  this  into  Eqn  2  and  separate  the  variables  we  have 

rE*  (-•*)  dEtl 

exp (-BE,\)  -  AE,i 


fS+o  t  —  f 
Jo 


where  A  —  g/fS+0  and  B  =  —jT'x. 


(8) 

(9) 


This  integral  cannot  be  solved  analytically  unless  we  can  expand  the  exponential 
term  in  a  series  and  ignore  most  higher  orders.  For  a  crystal  of  length  1  cm, 
an  applied  field  of  10®  V/m,  assuming  E,x  ~  mEt0  where  m  is  the  modulation 
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ratio,  and  taking  the  crystal  parameters  of  Table  2.1  we  find  \BEal\  %  4m.  Thus 
BE, i|  <C  1  for  sufficiently  small  modulation  ratios  (which  we  assume  is  the  case  in 
the  undepleted  pump  approximation).  A  similar  check  of  crystal  parameter  values 
shows  that  |£|/|A|  <C  1. 

Expanding  the  exponential  to  the  first  order  term  and  integrating  gives 

(1  -  exp(-/S+0(A  4-  B)t ))  (10) 

We  can  neglect  B  with  respect  to  A  and  we  reproduce  Eqn  3,  and  thus  the  rest  of 
the  derivation  of  Eqn  4  follows. 
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